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Shock Waves in Air at Very High Pressures * 
By.Ds KR. DAVIES 


Communicated by W. G. Penney ; MS. received 5 September 1947, and in 
amended form 18 February 1948 ; read 9 January 1948 


ABSTRACT. ‘The equations of the shock wave in air at N.T.P. are solved numerically for 
pressures on the high pressure side up to 1000 atmospheres. The solution enables all the 
physical entities, such as temperature, wave velocity, particle velocity, air density, etc., on 
the shock wave front to be expressed numerically as a function of the wave pressure. 'To 
solve these equations it is necessary to know the internal energy (Z) and volume (wv) of one 
gramme of air over certain regions of the two-dimensional range 1<p<1000 atmospheres 
and 273 <T<16000° k. Calculations are carried out to assess the numbers of the various 
types of molecules, atoms and ions present at any p and J. ‘The only E& values needed at 
high p are those for which T is also high, and the simple gas laws may therefore be assumed 
for the volume determinations. The resultant energy is obtained by summing the internal 
energies of the different groups of atoms and molecules present. Statistical mechanics 
furnishes equations whose solution fixes the composition. ‘This set of equations is difficult 
to solve to any degree of accuracy if oxides of nitrogen are taken into account ; of these only 
nitric oxide, NO, ever exceeds a concentration of 1% by weight, and its maximum concen- 
tration is less than 5%. As a first approximation, the shock wave equations are worked 
out on the assumption that the species present are Nz, O., N and O. ‘These calculations 
make full allowance for all quantum states apart from ionic states. They are refined later 
by ee a consideration the presence of NO and argon (1:3°% by weight) and ionization 
possibilities. 


Salse WIN TROMB HOC IDIOMS] 


URING the early part of the recent war the problem of working out theor- 
etically the distribution of pressure and velocity behind a spherical shock 
wave in water was resolved by solving the hydrodynamical equations 

for the motion of water behind the waves by numerical methods. ‘The similar 
problem of propagation of shock waves in azr, however, proved to be much more 
difficult, as the equation of state of air at high temperatures presented considerable 
difficulties. The first step in the solution consisted of working out the relationship 
between shock wave velocity, mass velocity, temperature, pressure, etc., on the 
shock wave itself, assuming a sudden discontinuity at the wave front. It is 
this problem which is described in the paper. 

A system of equations was derived by Lewis and von Elbe (1938) relating 
the various physical entities connected with a shock wave front. ‘The notation 
used in the paper for these physical entities is as follows : shock wave velocity 
U, mass velocity u, temperature 7, volume of unit mas3 of air v, internal 
energy per unit mass of air /, and the air pressure p. These shock wave equations 
may be solved if we have a complete knowledge of the internal energy of air in 
the neighbourhood of the shock wave conditions. ‘They completely determine 


* This paper forms part of a thesis presented to the University of London for the Ph.D. degree. 
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all the variables U, u, T, p, v, etc., af one of the set is known, in terms of the known 
conditions ahead of the wave front. It is assumed in this paper that the air 
on the low pressure side is at N.T.p. The independent variable chosen is p, 
the pressure on the high pressure side. Provided that only differences of pressure 
amounting to not more than a few atmospheres on the two sides of the wave 
are considered, simple calculations lead to reliable values of U, u, and T, etc., in 
terms of p, as worked out by Lewis and von Elbe. 

Solution becomes difficult for p values above five atmospheres due to the 
variation of the specific heats of air from the values obtained at ordinary tempera- 
tures and pressures, as very high temperatures are generated on the wave front 
at high pressures. For pressures up to about thirty atmospheres the only 
complication to the simple theory arises out of the necessity to allow for the 
vibrational energy of air molecules. This paper is concerned with the pressure 
range from thirty to one thousand atmospheres, where, due to several complica- 
tions, an elaborate calculation becomes necessary. Dissociation of the air 
molecules must be taken into account: molecules and atoms appear on electronic 
energy levels other than their ground states, and allowances must be made, too, 
for the changing composition of the air as the pressure and temperature conditions 
change. At pressures greater than one thousand atmospheres ionization com- 
plicates matters further, although dissociation has ceased to be troublesome. 
These calculations are not treated in this paper. 

It becomes necessary as a first step to evaluate the internal energy F and the 
volume wv of one gramme of air over most of the range 1<p<1000 atmospheres 
and 273<7T<16000°xK. At this first stage in the solution, standard methods 
of statistical mechanics furnish equations whose solution fixes the composition 
of the air. If great accuracy is required, these equations provide numerical 
difficulties. In addition to the formation of nitrogen and oxygen atoms, argon 
and oxides of nitrogen have to be dealt with. However, only nitric oxide, NO, 
and argon, A, ever appear in concentrations greater than 1°%4 by weight. Once 
the composition has been fixed, the volume follows immediately from the law 
p=nkt; and the internal energy is the sum of the internal energies of the different 
groups of atoms and molecules present. 

As a first approximation only, nitrogen and oxygen molecules and atoms are 
considered, and E and wv are evaluated over the range specified above. These 
functions enable the shock wave equations to be solved and thus specify approxi- 
mately the neighbourhood within which the wave conditions lie. These calcula- 
tions are later refined by taking into account (a) formation of NO, (b) formation 
of argon and (c) a more recent evaluation by A. G. Gaydon and W. G. Penney 
(1945) of the dissociation energy of nitrogen. This later work is greatly simplified, 
as the results of the first approximation considerably narrow the two-dimensional 
range of T and p, within which the numerical work has to be carried out. 


§2, INTERNAL ENERGY PER GRAMME FOR THE MAIN 
CONSTI EUENDSSO ReAiRe Olea ax O00 sk 
Air at normal temperatures and pressures consists primarily of oxygen and 
nitrogen molecules with a small proportion of inert gases, mainly argon. The 
percentages by weight are : nitrogen 75-5, oxygen 23-2, argon1-3. Ata tempera- 
ture of about 2000°x. the oxygen commences to dissociate, and near 3000°K. 
the nitrogen constituent behaves similarly. At the 4000°. stage a considerable 
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proportion of the nitrogen atoms combine with oxygen atoms to form oxides 
of nitrogen. Of these only nitric oxide ever exceeds the figure of 1°, by weight. 
‘The components considered in this calculation are nitrogen and oxygen atoms 
and molecules, nitric oxide and argon. The proportion by weight of these 
gases, In one gramme of air, over the temperature range 4000-16000° k. varies 
considerably with temperature and pressure conditions. It becomes necessary, 
therefore, at the outset, to evaluate the energy per gramme of these components 
over the temperature range and later to assess the degrees of dissociation at 
each temperature. 

Statistical thermodynamic theory leads to simple formulae for the transla- 
tional, rotational, vibrational and electronic energies of the various components. 
‘These are: (i) the internal energy per molecule, or atom, associated with the three- 
‘dimensional degrees of freedom 


EBERT Oe, a te ih se (1) 


where R=Boltzmann’s constant; (ii) the internal energy per molecule of the 
two rotational degrees of freedom of a diatomic molecule, 


Eel ee See (2) 
(iu) the internal energy of vibration of a molecule, per molecule, 
He=nviiexp (Ry RG) tle) 9" es. (3) 


where v is the fundamental vibration frequency and h is Planck’s constant; 
(iv) the internal energy per molecule or atom, associated with the electronic 
motion 


E,= {Sg Wyexp(—W,/RT)}(Zgexp(—W,/RT)}, (4) 


where g;, is the statistical weight and W; the energy of the electronic levelz. ‘The 
ground state of an atom is taken to be D/2, where D is the energy of dissociation 
of the ground state of the molecule formed from two of the atoms. 

These formulae, and later ones covering the equilibrium constants which 
specify the degrees of dissociation, are based on certain necessary approximations. 
It is assumed that complete factorization of the partition functions, §3, upon 
which the energy functions and dissociation formulae are based, is exact, hence 
that the total energy is the sum of the separate energies. ‘This is not strictly 
true. The spacing of the vibrational levels is not constant, and the moment of 
inertia of the molecule changes with increasing vibration and from one electronic 
level to another. However, the errors introduced by these approximations are 
not large, especially in view of the order of accuracy required in the solution 


of the shock wave equations. 
Table 1 gives the spectroscopic data used in the energy and dissociation 


formulae. 


able 1 
D (ev.) v (cm?) Be (CONS) 
O, 5-082 1580 1-440 
No 9-764 2360 1:998 
NO 6:490 1900 1-684 


The relation between B, and J, the moment of inertia, is B,=h/(87*Je), 


where c is the velocity of light (3 x 10!°cm/sec.). 
8-2 
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Other spectroscopic data that we require may be summarized. The ground 
state of O is 8,; there is another state 1A, at 0-96ev., and a second excited state 
1y, at 1-62ev. The ground state of O is ?P; there are excited *D at 1-97 ev. and 
18 at4-5ev. The ground state of N, is #2 and the ground state of Nis*S. The 
excited states of N, and N are too high to be considered as they have very small 
Boltzmann factors. These ~X and 1S states of O, and O have little effect, but 
they have been included. The ground state of NO is *II and the first excited 
state 2%. High excited levels have no significant influence on the energy 
functions. 

The statistical weights of the excited states may be written down once the 
electronic levels are known. For example, statistical weight of the ground 
state of O is 9, because the level has spin degeneracy 3 (triplet nature) and orbital 
degeneracy 3 (P states). The multiplet structure has been neglected because 
the overall separation of the ?P multiplet is somewhat less than RT at temperatures. 
producing appreciable dissociation of the oxygen molecules. 

Statistical weights used for the ground state, first and second excited levels 
respectively are : for O,—3, 5, 1, for O—9, 5, 1 and for NO—2 and 2 for the ground. 
state and first excited level. 


Table 2. Internal energy of air constituents 


Total internal 


Substance IP (CO Te) peice Beane ae \t energy 
gm. 2 cal/gm.mol. (caligesy 
Oxygen 2,000 12,060 142 381 
molecules 6,000 39,600 4,975 1,393 
10,000 67,330 8,800 2,379 
Oxygen 4,000 11,910 86 750 
atoms 8,000 ; 23,830 1,454 1,580 
12,000 35,750 3,620 2,460: 
16,000 47,660 5722 Se! 
Nitrogen 2,000 11,450 0) 409 
molecules 6,000 38,710 0 1,382 
10,000 66,300 0) 2,368 
14,000 94,124 0 3,362. 
Nitrogen 2,000 5,960 0 426. 
atoms 6,000 17,870 () 1277 
10,000 29,790 0) DADS, 
14,000 41,707 0 2,978. 
Nitric 2,000 11,794 0) 393 
oxide 6,000 39,150 4 1,305. 
10,000 66,930 2 pail 2,238 
14,000 94,696 1,366 3,202 
Argon 2,000 5,957 0 149 
6,000 17,874 0) 448 
10,000 29,790 0) 746 
14,000 41,707 0) 1,193 


Note: Intermediate values for this and ensuing tables may be obtained on application | 
to the author. 


The formulae used for the electronic energies in cal/gm. atom or molecule : 
for oxygen molecules, 


Ey = 23-07 x 109(4-8 e~™ + 1-62e-)/(3 45-467). (5) 
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for oxygen atoms, 


Ey= 23-07 x 109(9-835 ee" 44-5 e-™)/(945e-M +e) >, (6) 
for nitric oxide molecules, 
Hy 23-07 x 109(5-453:e-%) (tema 9) ee... (7) 


where x, = P,/RT, x,=P,/RT, y,=Q,/RT, y,=O,/RT, z,=R,/RT, and P,, O,, 
R, are the energies of the first excited levels and P,, Q, the energies of the second 
excited levels, all in appropriate units. The relationship 23-07 cal/gm. atom 
= 1 ev. is used. 

In table 2 the energies per gramme of each constituent of air are given as a 
function of the temperature. 


§3. ESTIMATION OF THE EQUILIBRIUM CONSTANTS OF 
NITROGEN, OXYGEN AND NITRIC OXIDE UP TO 16000° K. 

It is possible with the aid of the general laws of thermodynamics to predict 
the form of the functional dependence of the degree of dissociation of a gas on 
the temperature and pressure conditions. ‘These laws, however, leave unspecified 
certain constants. A knowledge of the energy levels which the particles of the 
gas can occupy, as furnished by standard spectroscopic methods, enables these 
constants to be evaluated; the degree of dissociation can thus be computed 
numerically from thermodynamics in conjunction with spectroscopic data. 

The chemical reactions to be considered are :— 


©, O70 No — NEEN = NSIOE NO ye esaraee (8) 

Many more, whose effect is small and have been neglected, may occur, such as 

N,+O=N,0; N+O,=NO,; O=Ot +e. 

Consider the simplest type of dissociation equilibrium, AB= A +B, in which 
gaseous compound AB dissociates into the monatomic vapours A and B. Then 
the equilibrium constant K is defined by K=p, .pp/Pap, where py and pp are 
the partial pressures of the monatomic vapours and p,y, the partial pressure 


of the original compound. 
Using this definition in the nitrogen reaction N,=2N, it is easily shown that 


C= (ISG CRSP) pe ee ee (9) 
or In K, =In[4a?p,/(1 —a?)], where a, is the ratio=(weight of nitrogen atoms in 
a given volume)/(weight of nitrogen atoms plus molecules in the same volume), 
p, is the partial pressure of nitrogen (N+N,), and K, the equilibrium constant 
of the reaction at temperature 7. 

A similar equation holds for oxygen (a). Hence, once the K value is fixed 
at a given temperature, the dissociation, given by a, etc., can be evaluated. 

Calculation of equilibrium constant as a function of temperature is simple 
once the partition functions for the components have been obtained. 

The partition function P, for O and N may both be expressed as a product 
of three linear translational partition functions and one electronic partition 
function P,=73(m)P(A). A represents the atoms O or N, m the mass of the 
atom. 

Similarly the P, functions for nitrogen and oxygen molecules may be written 
P, = T3(2m) . R(1). V(v). P(M)/2, where T represents translation, R represents 
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end-over-end rotations, V represents vibration, and P(M) represents the electronic 
part of the partition function of a molecule M. The factor 4 enters for N, and 
O,, because only one half of the rotational phase space is available to these homo- 
nuclear molecules N, and O,. It does not enter for NO. ak 
These standard expressions are used: 
T(m) =v*\/(2amkT); RU) =20v/(2rIkT); 
P(O)=9 +5 exp( —Q;/RT) + exp(—Q,/RT)+...-. 3° 
:P(O,)=3 +5 exp(—P,/RT) +exp(—P,/RT)+ .. 225: 
P(NO)=2+Zexp(—R RD) + 2 5 
P(N) =4; P(N,)=1; V(v)=1/[1 —exp(—hv/RT)]. 
For the reaction ABA +B, the equation obtained for In K is 
InK = —D/kT +3/2|n[(27M,M,kT/h2M 3] +1nhB, +1n{1 —exp(—hv/kT)} 


[1 + (gi/ge) exp (~X/RT)I[1 + (eriga)exp(—V/RT)] 


fen [1+ (@[ed®)exp(-Z/kRT)] 


where D=energy of dissociation from the lowest vibrational level; M,, Mg= 
mass of atom A, B; M,,=(M, +M,)=mass of molecule AB; go, g, =statistical 
weight of atoms or molecules in the ground state and first level (the superscripts. 
are self-explanatory); X, Y, Z=energy increase from ground state to the first 
level for atoms A, B, molecules AB respectively. 

The second and higher electronic levels have no effect. On the basis of 
wave mechanics there must for N, and O, be added, to the right hand of the 
above equation, a term In2 to allow for the symmetry of molecules AA. 

The expression used for the oxygen reaction is 
InK= —D)/RT +3/2InamokT/h? + InhB,+I1n[1 — exp (hvo/RT)] +1n 54 

+ 2In[{1+0-556exp(—Q,/RT)] —In[1+ 1-667 exp( —P,/RT)], 


(m9 = mass of oxygen atom). Similarly for nitrogen 
InK= —Dy/RT + 3/2InamykT/h? + InhB, +1n[1 —exp(—hvy/RT)] +1n 32 


(mx = mass of nitrogen atoms), and for nitric oxide 
In K = —Dyo/RT + 3/2 In 2amymoRT/h2myo +1n hB,+1n[1 —exp(—hAvyo/kT)} 
+1n18+In[1+0-556 exp(—Q,/RT)] —In[1+exp { —R,/RT}]. 
EP (13) 


The B, in this formula is c times that given in table 1, since frequencies are 
usually expressed in wave numbers as in table 1. 


Table 3. Equilibrium constants of oxygen, nitrogen and nitric oxide as 
functions of temperature 


TCS K (O) (dyne/cm?) — K (N) (dyne/cm?) K (NO) (dyne/cm?) 


2,000 0:595 0 0 

4,000 Se IO 3-390 2°754 x 104 

8,000 6-740 x 10? 7°959 x 108 8715 «108 
12,000 Den tex < 1OE° 1-158 x 10° 32623 x 10% 


16,000 SLO20b allies 1-453 x 1022 2:620'<10%% 
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§4. ESTIMATION OF INTERNAL ENERGY E(p, 7), AND 
VOLUME, v(p,T), FOR ONE GRAMME OF AIR 

In this section the effects of oxides of nitrogen and argon are neglected and 

a D value of 7-35 ev. is assumed for nitrogen. The general formulae used are 


Lie Da al.) ECT) aa i tele. (14) 
and Un(p,T)=n(p, T), where n(p, 7) is the number of any class of atoms or 


molecules 7, E,(T) is the internal energy per particle of this type of particle, EF 
the total internal energy per gramme of air, and 

PU=an(peLjRIs oy 0 See (15) 
Fortunately the only values of EF required at high pressures are those for which 
the temperature is also high. At no point does v become less than about 1/60 
of the volume at N.T.P. so that the use of equation (15) is justified. 

By means of equations (9), (11) and (12), a, and a, may be evaluated over the 
dissociation region, and hence by (14) the energies per gramme of nitrogen. and 
per gramme of oxygen as functions of p, or py, and TJ. It can be shown that 

(Gitdy) py 02091 eas) Pom) 1 esa wle perce (16) 
where p,, p. are the partial pressures of nitrogen and oxygen respectively. | Hence 
Pp; and p, may be evaluated at each temperature for the required value of total 
pressure p(=p,+p.).  Interpolating at these p, and py, values in the energy 
curves for oxygen and nitrogen, and assuming that one gramme of air always 
contains 24°, by weight of oxygen and 76° by weight of nitrogen, we obtain 
the final curves for the internal energy per gramme of air for 273 <T<16000°K. 
and total pressures equal to 15, 25, 35, 45, 65, 75, 150, 200, 400, 600, 800 and 1000 
atmospheres; thus later interpolation between two p curves is easily performed. 

The general character of these results is seen in figure 1. Once the p, 
and p, values have been determined for a given p, the volume of one gramme 
of air is given by 

tkT 
‘ (pymy/(1 + 4) + (Pomo/(1 + a9) 
Some of these values are indicated in figures 1 and 2. 
It must be borne in mind, however, that these curves are intended to give 
an indication only of the actual trend of the E and v curves. ‘They depend on a 
previous determination of D for nitrogen and neglect the formation of oxides 
of nitrogen and the presence of inert gases, etc. 


§5. SOLUTION OF THE SHOCK WAVE EQUATIONS: 
FIRST APPROXIMATION 


Consider a coordinate system moving with the wave front. ‘The mechanism 
of transfer across the wave front is stationary, and consequently the mass, 
momentum and energy of a gas passing across the wave front must be conserved. 


According to Lewis and von Elbe (1938), 
ty |V1 =Up/V_3 Py +UG/0y=Po+UZ/y; Ey +uj/2 + py, = Ey + ug/2 + pods. 


From these it follows that 


REZ (Pye Py) (Cy sa) | es ee cine (18) 
U=u4,= VV {(P2 =P) (Cie Os) se ato (19) 
U =U, —Ug— V {(P2 —Pr)(@% =U ty a lesa (20) 


where p,, v, refer to the low-pressure side and p., v, to the high-pressure side. 
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internal energy per gramme of 
air, due to the passage of a shock 
wave, as a function of the pres- 
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calories 


p in atmospheres. 


AE is in 


per gramme of air; 


of one gramme of air as a 
function of the shock wave 
pressure (in atmos.). 
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These equations may be solved numerically by a process of successive 
approximation, starting with p, and assessing a value for v,. The volume and 
energy curves then give T, and E,. If on one side the air is at N.T.P., E, is known, 
and hence AZ. By comparison with the value of AE calculated from (18), 
adjustment may be made to v. Values of AE, Ty, U and u follow directly. 
These are given in table 4. 


Table 4. Solution of the shock wave equations (neglecting oxides of nitrogen 
and argon; and assuming D for nitrogen to be 7-35 ev.) 


Pp» (atmospheres) 30 100 500 1,000 
U (metres/sec.) 1,670 3,000 6,560 9,350 
wu (metres/sec.) 1,370 2,590 5,960 8,310 
IEDC Te) 1,480 3,600 7,950 15,000 
FE, (calories per gm. of air) 287 875 4,320 8,450 
UV» (cm® per gm. of air) 140 107 70 86 


§6. SOLUTION OF THE SHOCK WAVE EQUATIONS: 
INCORPORATING EFFECT OF NO AND A 

In this solution the value of D for nitrogen is 9-764. 

It has been shown that the effect of allowing for the formation of NO is a 
‘small one, resulting at its maximum (4000-5000°k.) in an increase in the 
value of EF, of about 12 calories and a decrease of p, of about 5 atmospheres. 
Allowance for argon has a similarly small effect. ‘The increase in the nitrogen 
dissociation value has a much larger influence on the shock wave conditions. 
‘The first solution is unchanged up to about 100 atmospheres. 

The dissociation equations are easily written down. Consider the species 
‘O,, O, Nz, N, NO and argon. ‘The last component does not change, so that 
its influence is easily predicted. If v is the volume of one gramme of air at 
temperature 7, « the number of oxygen atoms, present as O, O, or in NO, in 
one gramme of air, 8 the number of nitrogen atoms, present as N, N, or in NO, 
in one gramme of air, x, y the fraction of oxygen, nitrogen, in the form of atoms, 
and z is the number of molecules of NO, the equations relating x, y, z, v and 
pare 


pe [lee eloK ook Te | Lee ee (21) 
i=l ay eeblOks2bRT, 9 2) 9 Meee (22) 
Pa Geek li Kee Sy ae cae eh WT Gt ee (23) 
p'=(RT/2v)[a(e+1)+y(B+1)], sees (24) 
| es en ee cr (25) 


where K,, K,, Kz are the equilibrium constants at temperature T for oxygen, 

nitrogen and nitric oxide respectively, p,=partial pressure of the argon and 
p’ =partial pressure of remaining constituents. 

From (21), (22) and (23), 

_ (1 =s)(1 —a2/B)K/2K 3 — 20K 32"/BKy 

% (1 —az/B)Ky/2K5+aKy2/BKy 

: _ (1 =2)(1 —a2/8)Ky/2Ky—2aKy2"/BK 

ce Ap (1 —)K,/2K,+2 
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Substitution in (24) gives a bi-quadratic for z: 
Az'+B23+ Cz2*4+D2+H=0) > 40 ees (28) 


where 
+= EG) 2) 301 

ONE 0-B-B 
£09) -18)} OIE NEG 2 
NEO) G8 FB 
-(-8-8)E-2)) 
(2) 3) J] Ge) YEE AG a) a 
ELH) 0B ge) 


1 caw eat K 
D-—— 143) (2) (14 aa 
3( B/\Ks 2p") 


T!R 
a 
ean 

ea 


For a given T and p’ the biquadratic can be solved for 2; x, y and wv follow 
from equations (26), (27), and (23). From T and v, py, follows, and hence p: 
from (25). ‘The x, y, x factors fix the composition of air and hence the component 
energies. ‘The £, valueisthenthe sum of these. ‘This method is a very laborious. 
one as it entails initially the evaluation of the coefficients A, B, etc., which in 
itself is not easy, and then the solution of the biquadratic for a large number of 
p’ points. 

An alternative procedure was developed which proved much simpler. Fixing 
the temperature 7, a value of z is chosen which, with the aid of (26), (27), (28),. 
(24) and (25), yields a p, value in the neighbourhood of the shock wave conditions. 
It is here that the first solution is of use. With a little practice, choice of a = 
value at any 7, to give an approximate p, becomes an easy matter. FE, follows. 
from x, y, z and hence AF. AE is also computed from equation (18). The 
correct ‘“‘shock wave z’’ is the one which makes E equal by the two methods. 
The correct U, u, py, UV, etc., for the given 7, then follow from the shock wave 
equations. 

At 3000°k. the nitrogen component has not commenced to dissociate 
appreciably, so that no nitric oxide has been formed at this stage. At 4000°K 
air at shock wave conditions contains about 4-4% , at 5000°K. about 4-8% a 
weight of NO. The percentage then slowly decreases with increasing tempera- 
ture until at 14000°k. the content has fallen to 0- 27%. 


The results of these computations are shown in table 5. The temperature: 
T, in this case is the independent variable. 
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. At temperatures below 4000°x. the contribution made by the energy of 
dissociation to the internal energy is small, but at 14000°x. it isas much as 60%. 
It is computed by assessing the energy liberated by the initial dissociation of 
the nitrogen and oxygen molecules and subtracting a certain amount due to the 
recombination of some of the atoms to form nitric oxide. 


Table 5. Shock wave conditions in air 


Ti (OR) 4,000 - 6,000 8,000 10,000 14,000 
x 0-1571 0-7298 0-9257 0-9678 0-9917 
y 0000092 0-0092 0-0885 0-2968 07809 
z 0-1025 0-0875 0-0430 0-0226 0-0061 


°% mass in air of 


NO 4-46 3-81 1-87 0-98 0-27 
O 3-64 16-93 21-48 22-45 23-01 
O, 17-18 4-24 0-73 0-22 0-05. 
N 0-01 0-70 6-68 22:41 58:96 
Nz 73-41 73-02 67-94 52-64 16-41 
A 1-30 insOae 1-30 1-30 1:30: 
AE (cal/gm.) 1040 2024 3090 4790 8540 
p> (atrnos.) 127 236 366 561 1003 
vp (cm$/gm.) 91:7 81-8 78-4 72:5 71:9 
U (m/sec.) 3350 4540 5645 6965 9310 
u (m/sec.) 2950 4060 5075 6320 8445 


In table 6 the shock wave variables are given with the pressure as the inde- 
pendent variable. 


Table 6 
ps (atmos.) 200 400 600 800 1,000 
U (m/sec.) 4,180 5,890 7,200 8,310 9,310 
u (m/sec.) 3,750 5,300 6,550 7,560 8,450 
T> (°K.) 5,270 8,420 10,300 12,000 14,000 
E (cal/gm.) 1,730 3,280 5,100 6,670 8,420 
VU» (cm3/gm). 82-7 - 768 71:8 69-9 71-8 


The U and u values are not far removed from those given in table 4 as the: 
first solution. The 7, figures, however, have changed a good deal, and at 
p2=500 atmospheres 7, has increased by 1480°x. The AF values are still 
in close proximity and show a practically linear dependence on p,. ‘There is. 
now a minimum value of v, (about 70cm) at p,=800 atmospheres instead of 


at 600 atmospheres. 
§7. NORMAL IMPACT OF A SHOCK WAVE ON A RIGID WALL 


The pressure experienced by a rigid wall when a shock wave strikes it is. 
many times greater than the shock wave pressure p,. The simple theory is quite 
adequate for low values of pressure, but for the higher values calculation becomes. 
difficult. The method makes use of the fact that the air behind the returning 
shock wave must be at rest. Thus, if p. and v, are the pressure and volume of 
unit mass of air before impact, p, and v3 similar quantities after impact, wu the mass. 
velocity produced immediately in the rear of the shock wave by p, and v,, then 


n= Uy Ua ( Psa Paye Pm ee NS (29) 
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Outside the “dissociation region” solution is easy. Assume a value for 
ps then equation (29) gives vs and hence 7; by the perfect-gas laws. ‘The 
internal energy of air function, E(p, T), has to be plotted as a function of tempera- 
ture for equal values of p3—considerably greater in magnitude than was required 
for py. Interpolation yields (EZ; —E,) where E, and £; are the energies of air 
at the shock wave before and after impact. If these values agree with that obtained 
from the formula | 

E, —Es=4(Ps + Pe)(Y2—U3), inne (30) 
then the py chosen is the correct one. Normally p; has to be readjusted so that 
this requirement is satisfied. 

Inside the “dis- 
sociation zone’ solu- 
tion is more difficult, 
as v.(p, 7) has to be 
plotted for equal pres- 
sures (figure 8). The 
steps are : 

(1) Assume a cer- 
tain p,; value, equa- 
tiOn = (29) “givess das 
(ii) plot 8 (6) gives 
the associated 73; (iii) plot 8 (a) gives E,; (iv) hence £; —E,(known) 
=AE; (v) compare with value given by equation (30); (vi) adjust p, to obtain 
agreement. 

In the absence of NO and argon, and with the old value of D for nitrogen, 
the following values for p3/p. have been evaluated. 


(b) 


Figure 8. 


Table 7 
p (atmos.) iS 30 50 65 100 200 300 375 500 800 1000 
P3/Deo 5:87 6:83 7:8 8:2 8-9 10:3 12:8" 13-8" 14-08 1233i0e8 


It is unlikely that correction for formation of NO, etc., will affect these values 
to any appreciable extent. 


§8. EFFECT OF IONIZATION 

The method used for investigating the possible changes in the internal energy, 
equilibrium constants and shock wave properties, due to ionization in the earth’s 
surface atmosphere due to the very high temperatures considered in this paper, 
is similar to that for evaluating the ‘“‘electronic’’ contributions. Possible ionic 
formations are accounted for by a similar addition to the partition functions 
and yield contributions as follows: 

Ionic contribution to the internal energy 

E,= {tT exp(—L/RT)}gottexp(—G/RT)};  —...... (31) 

contribution to the previous In K formula. 


oe CAN A “B B 
| [1+ (exe gr) | [1+ (oe zh) | | 
In { = SO ee fe (32) 
| [P+ Bloom) | 


where 7, is the statistical weight of atoms (or molecules) in the main ionized state 
and I, the energy difference from ground state to ionized state. 
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These equations cover the transformations 
N,=Nz +e, O,=Of +e, O=O++e, N=Nt++e. —...... (33) 


The spectroscopic data required are taken from Molecular Spectra and 
Molecular Structure of Diatomic Molecules by Herzberg, and Atomic Energy 
States by Bocher and Goudsmidt. 


Electronic potentials (ev.) State Statistical weights 
Or 13-60 as + 
N* 14-54 Be 1 
OF 12-20 2118 2 
INT 15-58 ESS 2 


These electronic potentials are very considerably higher than the excited 
levels of the un-ionized gases. Indeed, substitution of these numerical values. 
in equations (31) and (32) shows immediately that ionic formation at high tem- 
peratures has a negligibly small effect on the internal energies and K values. 
It does not, therefore, affect the shock wave properties. 

For temperatures greater than 16000°k. ionization plays an increasing part; 
the calculations to allow for this are described by Dr. Kynch in his papers in 
the ‘‘Atomic Energy” series. The problem at these very high temperatures is. 
eased by complete dissociation. 
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DESC US SLON: 


Dr. A. H. Davis. The present paper deals with the relations which exist between 
pressure, wave velocity, particle velocity, etc., whereby the observations of one quantity at a 
point—say of the velocity of the wave—makes it possible to infer the magnitude of others— 
say of the pressure, temperature cycle and so on. It does not deal with the propagation of 
shock waves through the atmosphere. | Moreover, whilst theoretical work has been con- 
ducted on the propagation of sounds of “‘ small’? amplitude through the atmosphere, 
around obstacles and so on, I believe that there is not yet any theory of the propagation of 
finite shock waves under similar conditions, and consequently an interesting theoretical field 
lies open. 

The genesis of a shock wave may be observed by exploding an explosive mixture of 
hydrogen and oxygen contained in an india-rubber balloon of about 6 cubic feet in volume. 
(The operations of filling the balloon, first partially with one gas,.and then finally with the 
other, must be conducted by remote control methods, as the balloon is liable to detonate 
when handled under some atmospheric conditions.) Observations of the blast pressure at 
distances of some 3 ft. to 10 ft. from the explosion reveal a wave which is first rather non- 
descript, but which develops into a shock wave after travelling a few feet. 


Mr. A. Warren. I should like to draw attention to some unusual features in the shape 
of the shock waves recorded near the ground with small high explosive charges. 

The measurements were made (using an experimental technique recently described by 
Grime and Sheard (Proc. Roy. Soc. A, 1946, 187, 357) up to a distance at which the calculated 
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velocity of the shock wave is within }% of the normal velocity of sound, and under these 
circumstances it was anticipated that any change in the wave form of the shock pulse would 
be insignificant. : 

It was found that up to a certain distance a wave-form characteristic of normal shock 
propagation was recorded, but that subsequently a change of waveform occurred ; a 
secondary pressure peak commenced to form behind the shock front and the two peaks 
tended to separate. Since the maximum excess pressure in the initial shock front decayed at 
an abnormally high rate, the maximum excess pressure in the secondary peak soon exceeded 
that in the shock front, and in some instances the initial shock front almost disappeared at the 
greater distances. , 

It was observed that the critical distance beyond which the secondary peak pressure 
exceeded that in the shock front appeared to be greater when the propagation losses in the 
wave front were reduced. For example, it was greater (a) in the direction of the prevailing 
wind, (b) for gauges at (say) 3 feet above ground as compared with similar measurements 
very near the ground, and (c) if the surface of the ground between the charge and gauge was 
covered with sheets of smooth material, e.g. strawboard with improved reflection character- 
istics. 

It was known from the recent series of papers by Friedlander (‘“‘ Diffraction of sound 
pulse’’, Proc. Roy. Soc. A, 1946, 186, 322) that after diffraction from obstacles of various types 
a shock wave was likely to lose sharpness and become broadened. A similar effect was 
anticipated in an unpublished mathematical treatment by Prof. Sir. G. I. Taylor dealing 
with the effects due to turbulence, and to the gradient of wind velocity near the surface of the 
ground. 

It appeared from the records already obtained that modified wave-forms very similar 
to those to be anticipated on theoretical grounds from the effects of diffraction or turbulence ~ 
were observed, but that in addition, under certain circumstances, this modified wave-form 
could co-exist with an initial shock front of characteristic shape at an appropriate time 
separation. 


Mr. A. R. Bryant. I should like to ask what the temperature would be if an explosive 
were surrounded by an envelope of argon instead of air. One would expect this to be much 
higher owing to argon having a higher value of y and owing to absence of dissociation. 
With regard to the results of the shock wave calculations could the author state how sensitive 
the quantities behind the shock front are to the assumption that thermodynamic equilibrium 
is or is not attained ? 


AvTHOR’S reply. I have read with interest the remarks of Dr. A. H. Davis and Mr. A. 
Warren. An approximate solution of the problems of hydrodynamical flow behind a shock 
wave set up in air as a result of an explosion from a spherical charge was, in fact, attempted in 
1941 by Dr. W. G. Penney and myself (unpublished). At that time the adiabatics of air at 
high temperature and pressure had not been evaluated, and y was taken to be equal to 1-403. 
This assumption was likely to introduce large numerical errors, but the results showed the 
presence of a secondary pressure peak, as mentioned by Mr. Warren. The problem was 
investigated during 1943-45 in the U.S.A., but I am not aware if any of this work has yet 
been published. 

In reply to Mr. Bryant’s first question, I agree that the temperature would be higher if 
the explosive were surrounded by argon of the same density as air in the normal atmosphere, 
the internal energy of argon per gramme at a given temperature being less than that of air ; 
there is no rotational, vibrational energy, and no dissociation to be allowed for. 

With regard to Mr. Bryant’s second question, the shock wave velocity and pressure are 
probably very little affected by the assumption of thermodynamic equilibrium. The 
actual temperature, however, will be higher than those obtained on this assumption. In 
the paper it is assumed that the energy is taken up in the translational, rotational and 
vibrational forms, distributed according to the equations of statistical mechanics. A given 
volume of air receives a definite quantity of energy from the “‘ shock wave ’’, and it is probable 
that, due to the nature of the “‘ shock ’’, rather more energy is taken up into the translational 
type and less into the rotational and vibrational types than given by the theory. Hence the 
temperature would be higher than expected on the assumption of thermodynamic equili- 
brium. 
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ABSTRACT. A plane shock wave falls normally on the interface between two ideal gases 
of molecular weights M, m and constant ratios I’, y of specific heats. Asa rule, the character 
of the reflected wave then depends only on the ratio 1M/ym. In exceptional cases, however, 


the type of reflection may change at a critical value of the incident intensity. An example of 
this is given. 


of specific heats are at rest and separated by an ideally thin and light 

partition. ‘l’he pressure is py) and temperature 7, in each gas. A plane 
shock wave S,, the steady conditions behind which are indicated by suffix 1, passes 
through the first gas and falls normally on the interface. 

A wave S, will then be transmitted through the second gas; we denote 
conditions behind S, by a sufhx 2. These conditions are defined by the familiar 
Rankine-Hugonuot relation if S, is a shock wave, or by the Riemann relation for a 
progressive wave if S,isararefaction. ‘Thus, if W represents in general the mass 
velocity behind a wave (measured in the direction of propagation of S,), and z the 
pressure ratio across it, so that 7=p,/po, 72=P2/po, we have for S, 


To perfect gases with specific volumes Vy and vp and constant ratios I, y 


. Di APG osry py) Aa Re (1) 
and for S. 
( V{A —1)poro}(2—-1)//Ar2+1), if m>1 
2 | —V{Q?= pov} ne], if mcd “ 


where A=([C'+1)/((.—-1) and A=(y+1)/(y—1). It is of course physically 
obvious, and will be shown below, that S, must be a shock wave (7, >1) and not a 
rarefaction wave (7,<1). However, for the moment we need not assume this. 
A reflected wave S, will also arise in the first gas; we denote conditions behind 
by a suffix 3, so that 73;=3/p,;. S3 may be either a shock wave (73>1).or a 
rarefaction wave (73<1), and the mass velocity behind it is 
[Wy -V(A- pV} -D/V{Aws +1}, if ms >1 


a ee re RO ONE es ote Teo Raa 3 
; 1 W, + /{(A? -1)p,V3[1 —7"4 $Y], if 73<1. ©) 


Continuity of velocity and pressure across the interface requires, however, 
that W,=W, and also p2=ps, so that 7,=737,. Making these substitutions in 
(2) and using (3) and (1), we are in a position to determine 73, and so all the pro- 
perties of the reflected and transmitted waves. The procedure is now well 
known. 

It may easily be shown that Wy, as given by (2) with 7, =737,, 18 a monotonic 
function of 73 which increases continuously through zero as 73 increases through 
1/7,. In the same way, W, is a monotonic function of 7; which decreases con- 
tinuously through W, as 73 increases through unity. It therefore appears at once, 
since W, >0, 

(a) that 7, cannot be less than 1/7; in other words, that 7, must always 

exceed 1, so that the transmitted wave is necessarily a shock wave ; 
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(b) that 7, may be either greater than or less than 1, so that the reflected wave 

may be of either type, according to circumstances. 

These properties are illustrated in figures 1 and 2 by the case of a 2-atm. 
shock wave proceeding from oxygen to argon and vice versa, both gases being 
initially at 1 atm. and 273°K. 

The condition that S, be a shock wave is that the two curves intersect at a 
point for which 7, >1; but in view of the above remarks concerning the course of 
W,(73) and W,(73), this condition can be replaced by the analytically much simpler 
one that W,(73) >W,(73) when 73=1; in other words, that 


(A —1)\Vy)(Anys1) = =) Oa ae (4) 
which may be written 3 
(ae, — DCT EIA, Ger Ls) 219 (5) 
where A, V,=1 =5 9%. 
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Figure 1. Reflection of a 2-atm. shock wave in Figure 2. Reflection of a 2-atm. shock wave in 


oxygen by argon. py=1 atm.; 7>=273°k; argon’ by oxygen. pyp=1atm.; T)=273° k. 


W,, W,=mass velocities behind transmitted 


2, W;=mass velocities behind transmitted 


and reflected waves; 73=pressure ratio and reflected waves ; 73=pressure ratio across 
across reflected wave, which is seen to be a reflected wave, which is seen to be a rare- 
shock. faction. 
Shock Shock 
33 - — Rarefaction = = Rarefaction 
Four alternatives now arise: 
(a) TAg<yd9 and (I'+1)Ap <(y+1)8o. 


Then (5) is satisfied irrespective of 7,(>1), and the reflected wave is always a 
shock wave. 

(b) TAy>y5y and (I +1)Ay>(y+ 1d. : 
Then (5) cannot be satisfied for any 7,(>1), and the reflected wave is alwavs a 
rarefaction wave. f 

(c) PAy<ydp bué (0+ 1)Ap>(y + 1)8o. 
Then (5) is satisfied for 7,<7,’, but not for 7,>7,’, where 

Bie (y—1)6,—(T—1)Ay 
+ (T+ DAG (y + So" 
Thus, for 7,<7,', 5; is a shock wave, but for 7,>7,’, S3 is a rarefaction wave. - 

(d) TAp>y59 but (1'+1)Ap<(y+1)8,. 

Then, for 7,<7y’, S; is a rarefaction wave, but for 7;>7,’, S, is a shock wave. 


Reflection of a plane shock wave at a gaseous interface E26 


Cases (a) and (b) are evidently reciprocal. Thus they show that if shock 
Waves in the first gas are always reflected as shock waves by the second, shock waves 
in the second will always be reflected as rarefaction waves by the first, and vice 
versa. ‘This is physically reasonable. Cases (c) and (d) are reciprocal in a similar 
sense, but the existence of a critical incident pressure ratio 74’, at which the 


‘reflection changes type, is more surprising. It remains to be seen whether this 


effect can be realized in real gases. The necessary condition may be written, for 
convenience, 


TMSym, CMe (744) eee (7) 
where M, m denote molecular weights. Examination of tabulated specific heats 
and molecular weights shows at once that (7) is not usually satisfied, so that 


cases (a) and (6) form the rule. However, the possibility of satisfying (7) is not 


excluded if I’ and y are sufficiently different. Suppose, for example, that 
P=1-1 and y=1-667. Then (7) requires that 1:-27m<M<1-515m. It appears 
in fact that (7) is satisfied by the following pairs of gases: neon with acetylene, 
argon with butane vapour or methyl chloride vapour, for which approximate 
constants at 1 atm. and 273°x. and corresponding values of 7,’ are given in the 
following table: 


Gas IP A rA (r-+1)A 7" 
Neon (Ne) 1-642 0-966 (g/1.) 1-586 2°553 (with Ne) 
Acetylene 1.28 1-190 1-52 Dial 1-78 
Argon (A) 1-667 1-781 2-964 4-750 (with A) 
Butane vapour 1-11 2:°599 2-88 5-48 los) 
Methyl chloride vapour V7) 2-238 2:862 5-10 1-60 


Calculations for the vapours =00. 
are rather unreliable, in view of Reflected 
the high critical temperatures and 
correspondingly large departures 
from the ideal gas laws. For 
neon and actylene, however, the 
effect is probably authentic, and 
could perhaps be detected in 
Schlieren photographs by the 
reversal of the sign of 7,—1 at 
q,/ ; the wave- and mass-velocities 
do not, of course, change sign. 
The properties of transmitted 
and reflected waves in _ this 
system up to 7,=2:'5 are shown 


Reflected 
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I 
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in figure 3. p, (atm) 
For the great majority of Figure 3. Reflection of weak shock waves in acetylene by 
j neon. =1 atm.; T)>=273° K.; pi, p3=pressures,, 
Berens ‘ pa ae (7) gales aod and W, on mass polos behind Ro ied and 
the condition fora reflected shock between transmitted and reflected waves respectively; 
is then simply TM <ym, which D,, D2, D;=absolute velocities of incident, trans- 
may also be written A.M <domi or mitted and reflected waves ; d;=velocity of reflected 


wave relative to gas behind incident wave. Velocities. 


are positive when in the direction of incident wave. 


— ? 


e (5) 
Apo] s 


ar) 
Opo/s 


where Ay, a) denote velocities of sound, and S entropy. 
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ABSTRACT. 'The rate of occurrence of transient bursts of ionization within the lower 
E region of the ionosphere has been systematically measured during the period January 1945 
to July 1946. 

It has been found that the activity varies both diurnally and seasonally and in such a 
manner as to lend support to the meteor theory of burst formation. Some measure of solar 
control of the rate of burst occurrence has also been detected ; this effect may be explained 
in terms of the normal diurnal E region changes in density of ionization. From experiments 
performed during the solar eclipse of July 1945, it is concluded that the sun does not emit 
a burst-producing radiation. 

Analysis of the observations suggests that the greater proportion of the bursts are 
created within a thin layer located at a height of 86 km. ; the distribution of the bursts 
within this layer has proved to be uniform over wide areas and no latitude effect has been 
detected. It is established that the rate of incidence of bursts which present echoing areas 


between A and A-+dA m? to a radio wave of frequency v may be expressed in the form 
CdA|v?A!, with C constant. 


si INTRODUCTION 


HE density of ionization within the E region of the ionosphere normally 

varies in a regular manner with the intensity of the ultra-violet light 

received from the sun, but a number of workers have shown the temporary 
existence within the layer of regions of much higher density. These disturbances 
appear to be of two distinct types. Appleton and Naismith (1935, 1940) describe 
one form of this sporadic ionization that has come to be known as “abnormal 
or intense E”; this type usually persists for minutes or even hours and may 
extend over large areas. A second type of E-layer disturbance was also recognized 
by Appleton, Naismith and Ingram (1937) in the form of localized bursts of 
ionization whose density greatly exceeded that characteristic of abnormal E. 
These ion clouds were observed to occur both by day and by night and yielded | 


transient echoes of only about a second duration. The reflection coefficients |] 


and heights of occurrence of these clouds were later studied in more detail by 


Appleton and Piddington (1938). It is with this second type of E-layer distur- |] 


bance, or “ionospheric burst” as it is now commonly termed, that the present i | 
paper is concerned. 

The problem of E-layer disturbances was examined by Eckersley (1940) 
from the point of view of radio communication effects and he has shown that || 
the presence of signals within the skip zone is attributable to scattering from | 


ionospheric “clouds”. ‘The rate of occurrence of the clouds was found to II 


follow a diurnal cycle, while the dependence of the intensity of the scattered ] 
radiation upon the wavelength was in agreement with a fourth power law pre- | 
viously suggested by him (1932). | 

The nature of the ionizing agent responsible for these E-layer perturbations 
is a question which has given rise to speculation, but much of the evidence now 
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points towards the validity of the meteor theory. Nagaoka (1929) was the first 
to consider the possible influence of meteors on radio transmission and suggested 
that disturbed ionospheric conditions might arise from the removal of electrons 
by a process of attachment to particles of meteoric dust. Skellett (193:25,1935) 
disagreed with this theory and, supported by experimental work performed in 
collaboration with Schafer and Goodall (1932), showed that impact ionization 
by meteors was a more likely explanation of the observed phenomena. Pierce 
(1938) also found that sudden increases in the abnormal E type of ionization 
could be explained by the passage overhead of a single large meteor, while certain 
experiments reported by Eckersley (1937) on the occurrence of scattering clouds 
were considered by Skellett (1938) to be consistent with the meteor theory of 
their formation. Further convincing evidence in support of the meteor hypo- 
thesis has recently been provided by Appleton and Naismith (1946), and by 
Hey and Stewart (1946), but that certain difficulties still remain in applying 
the theory is shown by the work of Eckersley and Farmer (1945) on the rapid 
variation of the direction of arrival and state of polarization of the waves scattered 
from the bursts. 

The aim of the present investigation has been to gather further information 
upon the characteristics of ionospheric bursts by making a systematic study 
of the time variation of burst activity. Observations were also made during 
the total solar eclipse of July 1945 to detect any influence which the sun might 
exert upon the phenomenon of burst occurrence. 


§2. EXPERIMENTAL METHODS 


The ‘“‘pulse-echo”’ technique has been used in the present work for the 
detection and location of ionospheric bursts. ‘The experiments have been 
conducted at a number of Royal Air Force radar stations each of which was 
equipped with a high power, pulsed radio transmitter and a sensitive receiver. 
‘The transmitter aerial consisted of a vertical stack of horizontal, half-wave dipoles 
and was situated within 200 metres of a receiver aerial with dipole elements 
parallel to and at the same mean height as those of the transmitter. ‘The radio- 
frequencies employed ranged from 20 to 45 Mc/sec. but the pulse recurrence 
frequency was 25 pulses/sec. at all the stations. 

The echo signal returned by the burst was recorded photographically in the 
usual (P’t) type of record. With the aerial system described above most of 
the energy is radiated at comparatively low angles of elevation; a method of 
oblique incidence has therefore been used, and the P’ of the records represents 
not the usual equivalent height of the ionospheric reflecting layer but the slant 
range of the ionospheric burst. The magnitudes of the echoes were also recorded 
photographically using a ciné camera which was operated without its shutter. 
It was found that successive sweeps on the cathode-ray display tube were regis- 
tered as separate parallel traces on the film so that the echo amplitudes could easily 
be measured from the trace deflections. This simple technique permitted 
the rates of growth and decay of particular responses to be studied. 


§3. DIURNAL AND SEASONAL VARIATIONS 
In order to investigate diurnal and seasonal changes in the rate of burst 
occurrence observations were made from two identical stations at Bawdsey, 
Suffolk, on a frequency of 22:69 Mc/sec., a pulse width of 8 psec. and a power 


9-2 
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of IMw. The hourly rate of occurrence of bursts formed within a selected 
slant range interval was measured over periods of twenty-four hours. The 
observations covered the period January 1945-July 1946, with the 24-hour 
sequence of observations taking place twice each week. 

Every precaution was taken to ensure constant performance of the equipment 
but output power variations did occur and, in order to apply appropriate correc- 
tion, the relationship between the number of bursts observed and the power 
radiated was separately investigated. This work established the result that, 
for a given system and range interval, the number of bursts observed was directly 
proportional to the square root of the power. From this result it was now 
possible to reduce the average measured power over an interval to the standard 
power of 1 Mw. and to correct the burst rate accordingly. 

The number of bursts recorded during the successive hours of a day’s run 
clearly showed the presence of a diurnal cycle of burst activity. ‘These daily 
curves possessed certain features in common but marked changes also occurred 
from day to day and curves of mean values for each month were therefore plotted 
(figures 1 and 2). The curves of figure 1 show relative monthly changes in the 
daily cycle; the absolute levels of activity are not comparable since the same 
aerial system was not used from month to month. The curves of figure 2, for 
which the same aerial system was used, are however, strictly comparable both 
for total activity and for relative monthly changes. 

hicuress le@ande7 
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Figure 1. 
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approximately 0400 and 1200 vu.t., and do not appear to vary with time of 
sunrise and sunset through the year. The relative amplitudes of the two 
maxima are, however, dependent upon the season of the year; in particular, 
the early morning peak tends to be associated with the period of the 
winter solstice while the late morning peak becomes the dominant feature of 
the curve about the time of the summer solstice. This is seen in the middle 
column of curves in figure 1 : December 1945 and June 1946 show pronounced 
early and late morning peaks respectively, while September 1945 and March 
1946 have transitional types of curves with both peaks present and of roughly 
the same amplitudes. 

The change in total volume of daily burst activity throughout the year was 
studied by taking observations on convenient days when the same aerial system 
could be used. The total number of echoes during a 24-hour period observed 
within the range interval 160-320km. was determined and the average number 
of responses per hour extracted. Figure 3 shows these average daily figures 
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plotted against the appropriate dates of observation. It is seen that a seasonal 
cycle of activity is present with its maximum and minimum occurring about the 
times of the autumn and spring equinoxes respectively. 


§4. SOLAR ECLIPSE EXPERIMENTS 


An experimental investigation of solar effects on ionospheric burst formation 
was possible during the total solar eclipse of 9th July 1945. Such solar effects. 
might be expected to arise in two ways, either directly due to the emission of a 
burst-producing radiation from the sun or indirectly from the sun’s control of 
the ionized state of the ionosphere. In both these cases the rate of occurrence 
of the bursts would be the feature most likely to change during the course of a 
total eclipse of the sun. The general plan of the experiments, therefore, was 
to compare the rate of formation of the bursts on the eclipse day with the average 
rate obtained from observations made during two control periods each of seven 
days, on either side of the eclipse day. 

"The experiments were conducted from two groups of stations in the north 
of Scotland, lying within the belt of totality for a particle radiation at the range 
of velocities to be expected. The activity throughout the eclipse period was 
measured in terms of the number of echoes observed per half hour and was 
compared with the mean activity derived from a grouping of the corresponding: 
periods of the control days. The amplitudes and durations of echoes observed 
during the eclipse were similarly compared with the averages derived from the 
control periods. This procedure was followed at all the stations of table 1. 

A careful statistical analysis of all the results obtained failed to reveal any 
significant departure of iono- 
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spheric burst characteristics 
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a source of burst-producing 
particles. However, the great 
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esting conclusions to be drawn 
on the questions of time varia- 
tion and space distribution of 
burst activity. Figure 4 
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changes occur in an apparently random mannerat the different stations. The large 
differences between the mean activity levels of the stations taken over the whole 
duration of the observations are particularly noteworthy and will be discussed 
in more detail later. 


§5. RANGE DISTRIBUTION OF BURSTS 

From the (P’t) records obtained during the eclipse experiments a histogram 
was derived for each of the stations relating the number of echoes in successive 
16km. intervals of range with the mean ranges of these intervals. It was found 
that every station possessed a characteristic range-distribution histogram. 
Similarly, each of the various aerial systems used at Bawdsey was found to 
have a characteristic histogram; figure 5 shows a typical histogram secured 
from a two-stack receiver aerial system combined with an eight-stack transmitter 
array, both aerials centred on a mean height above sea-level of 255 ft. and operated 
on a frequency of 22:69 Mc/sec. Figure 6 shows the distribution corresponding 
to a single dipole receiver aerial and a single dipole transmitter, each at a height 
of one-quarter wavelength above the ground. 

In all these experiments a sharp cut-off in range was observed at about 80km., 
in agreement with the accepted view that the bursts are formed within the E 
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layer, but it was also noticed that echoes occasionally appeared at less than this 
range. Detailed analysis has shown that the form of a station’s range-distribution 
histogram may be convincingly explained in terms of the interaction of the station’s 
characteristic radiation pattern. with a burst-producing layer situated within 
the lower E region of the ionosphere. The sharpness of certain of the histogram 
minima, e.g. that at 650km. in figure 5, suggests that the bursts are probably 
confined to a very thin layer so that, to a first approximation, we may neglect 
a height distribution of bursts and seek for a theoretical development of the 
histograms on the assumption that the bursts are distributed uniformly over a 
spherical surface. 

The number of bursts occurring within the range interval R to R+dR and 
having equivalent echoing areas between A and A+dA square metres will be 
proportional to the zone of the ionosphere intercepted between these ranges 
and also to the occurrence rate per unit area of bursts of this magnitude, say 
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ndA. If his the height of the burst-producing layer and p be the radius of the 

earth, then the area of the zone is 27(1 +//p)RdR, and we may therefore write : 
Number of bursts observed in the range 

interval R to R+dR having echoing areas ocnRdAdR,  ...... (1) 

between A and A+dA: 

i.e., total number of bursts observed in e [aa [, pat apn ee (2) 

the range interval R, to Ry: JR J AL 

In this expression A,, is the smallest echoing area which can be detected by the 

given equipment at the range R and it will be shown later (equation (9)) that 


|, ndA cl|/Ay, 
Ay, 
so that (2) may be rewritten: 
Total number of bursts observed in the ae = 

5A Rope 3 
interval R, to R,: oe R, OEE On (3) 
Now a burst of echoing area Am? formed at a range R from a station radiating 
a power of P kilowatts on a wavelength of A m. will return an echo whose 
signal/noise ratio is Z, where Z is given by 


Z cc AR24/(APnyn,)b(«); 


in which m and m, are the numbers of dipoles composing the transmitter and 
receiver aerials respectively and ¢(«) is the product of the interference and 
grating factors of the arrays, expressed as a function of the angle of elevation «. 
For a given elevation this function varies with the azimuth, due to the form of 
the radiation pattern from a half-wave dipole, and it is therefore necessary to 
work with the azimuthal average value of ¢(«). 

Since the smallest echo signal which can be distinguished from the noise 
corresponds to a constant valué of Z, it follows from (4) that the least echoing 
area, A;, which can be detected at a range R and elevation « is expressed by 


Ay* cARS@4/(Pnn,)o(), a ee (5) 
Substituting this value of A, in (3) gives 


Total number of bursts observed in the range interval R, to R, 


R; 
oc | AR-1\/(Pnyn,)b(a)dR, 
R, 
which may be expressed 


R, 
oc | O(R, h)dR, (7) 
Ry 


cee eee 


where his the height of the layer of burst formation. 

In order to develop the range distribution histogram for a particular system 
the appropriate Q(R, h) was evaluated and plotted against R for a selected value 
of h, suitable correction being applied for the effects of atmospheric refraction. 
The resulting curve was integrated over successive 16km. range intervals by 
use of Simpson’s Rule and the areas of the consecutive blocks compared with 
the experimental histogram as in figure 5, where the chain lines show the theoretical 
histogram derived for the Bawdsey system. Theoretical histograms were simi- 
larly calculated for all the stations and for a number of selected values of h; 
comparison of the theoretical and experimental patterns showed that the most 
probable value for the height of the burst-forming layer was 86+2km. The 
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question of a distribution of the bursts in height was examined by a method 
which was an elaboration of that given above, but no better measure of agreement 
with experiment was obtained than that yielded by the simple assumption of a 
thin layer of origin. It was therefore concluded that for oblique incidence 
and for radio frequencies in the range 20-45 Mc/sec. the bursts may be regarded 
as producing their scattered waves from within a very thin layer situated at a 


height of 86 km. 


§6. ECHOING AREAS OF BURSTS 


The range of echoing areas presented by ionospheric bursts was studied by 
the use of a simple arrangement of two parallel, horizontal half-wave dipoles 
each supported one quarter wavelength above the ground and operated as the 
transmitter and receiver aerials respectively. The radiation pattern of this 
system consisted of a single vertical lobe. In order to permit the estimation 
of actual echoing areas as distinct from relative values only, the system was 
standardized from observations made during test flights with an aircraft of 
known echoing area. 

(Pép)erecords of bursts 4,,, 
occurrences were secured 
with this system and the 
range distribution histo- 
gram compiled as in figure 
6. The mean rate of 
observed occurrences per 
unit area of the burst layer 
could then be calculated at 
ranges corresponding to the 
centres of the successive 
blocks of the histogram; 
assuming that the layer of 
occurrence is thin, for 
oblique incidence and at 
ranges much greater than 
86km. the volumes of 
successive zones of this 9, a 1 | 301 
layer will be proportional SOMA 
to the areas intercepted on Figure 6. Range distribution. 


the assumed surface of 
occurrence situated at a height of 86km. It is then sufficient to consider the 


superficial density of occurrence of the bursts on this equivalent surface. 

Using the constants of the standardized aerial system the values were calcu- 
lated of the smallest echoing areas that could be detected at points on the surface 
corresponding to the centres of the successive 16km. zones. Let the minimum 
detectable echoing areas be A’ and A” for two zones for which the occurrence 
rate of observed echoes per unit area are N’ and N” respectively. If we assume 
that the distribution of bursts over the ionosphere is uniform then (N’ —N’”) 
must be the rate of incidence per unit area of bursts having echoing areas between 


A and A”, 


2 


NUMBER OF ECHOES 


re 
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The number of bursts per unit area of the ionosphere occurring in unit time 
which present echoing areas between A and 4A+dA we have already defined 


as ndA, so that the number included in the interval A’ to A” will be | nd A, i.e. 


(N’ —N”)=AN= lis ndA. Now this is the area under the (n, A) curve included 


between the 4’ and A” ordinates, and since the mean ordinate over this interval, 
AN|(A’' —A”), is readily calculable, the form of the (7, A) curve itself may be 
traced by plotting the equivalent rectangles included between ordinates erected 
at abscissae corresponding to the limiting echoing areas A’, A”, etc., of the suc- 
cessive zones. In figure 7 this method of plotting has been pursued and the 
straight line indicated leads immediately to the relation 


ndA AW "dA = = i eee (8) 


It now follows that the number of bursts per unit area per unit time, Ny say, 
which present echoing areas greater than A will be given by 
oa) 


Ne { ndA ie .Nical/a/4 ee Se eee (9) 
A 


Separate confirmation of these results were obtained in the investigation on 
power effects previously mentioned. 
From a set of observations made 4 

during January 1946 ona frequency 
of 22:69 Mc/sec. the constant of 7+ 
proportionality in equation (9) 4 
proved to be 643, and the units of 
area and time 1000 square miles * *+ RSG 
and 50 hours respectively, which 4 
suggests that the total number of — 


bursts per day created in the iono- = ‘} 

sphere with echoing areas greater | 
than 10 m? is of the order 2 x 107. i tect a : J 
This incidence rate may be com- Figure 7. 


pared with Skellett’s assumption 

(1935) of 10° sporadic meteors per day; it must be observed, however, that due: 
to the comparatively low pulse recurrence frequency employed in the present 
work only echoes of duration greater than 0-3 sec. could be read from the films; 
also, no allowance has been made for a possible aspect effect such as was con- 
sidered by Pierce (1938). 

No upper limit has been observed for the echoing area which it is possible: 
for a burst to present; thus in figure 5 for example the smallest echoing area 
which could be detected at the range corresponding to the first histogram 
minimum, i.e. 650km., was of the order 10’ m2, yet the figure shows that such 
echoes did occur. 

A study of the rate of growth curves obtained from ciné observations on 
burst echoes shows that, in general, the echo grows rapidly in amplitude to a 
maximum during an interval of about 2/25sec. The subsequent decay process 
is much slower and is often closely exponential in form. Many echoes, on the 
other hand, do not show a smooth decay curve but one which oscillates in magni- 
tude. Amplitude-time curves of this type have been obtained for a large number 
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of echoes originating from within the same small range interval. ‘hese curves 
show that the duration of a burst is not simply related to the maximum echoing 
area which it presents during its lifetime ; they suggest that the processes of growth 
and decay of ionization must be subject to change from incident to incident, 
and that the manner of change in geometrical size of the ion cloud due to diffusion 
must also be variable. 


§7. RELATIONSHIP BETWEEN APPARENT RATE OF OCCURRENCE 
OF BURSTS AND RADIO-FREQUENCY OF OBSERVING STATION 


We have seen from figure 4 that the average level of activity during the eclipse 
experiments differed greatly from station to station, and it is of interest to examine 
whether this effect may be explained in terms of the known characteristics of the 
stations or whether the burst activity does indeed vary from locality to locality. 
In table 1 are set forth the average number of echoes per hour for each of the 


stations, together with the corresponding values of the integral /O(R, mak, 


for the appropriate range interval. 


Table 1 
: Frequency Echoes per hour 
Station ices Seen 20-220 fo, h) dR=b m=a/b 
Class A miles=a 
Whalehead 25°10 802 109 7:4 
Hillhead 1 22-69 607 80 7:6 
Sango 25-10 256 33 7:7 
Netherb’n 26°36 332 54 6:2 
Class B 
Hillhead 2 42-27 236 109 2:2 
Inverness 43-69 105 87 1-2 
Tannach 42-27 88 94 0-9 


The general conclusion drawn from this table is that the number of bursts. 
observed by a station diminishes with increase of the radio frequency. We 
have found previously, however, that the echoing area distribution function 
expressed by equation (8) holds good for different radio-frequencies but with a 
different constant of proportionality. In order to interpret the table therefore 
it has appeared reasonable to modify equation (8) so as to include the effect 
of frequency explicitly and to write it in the form 


pied Poa ee ee (10) 


where , is the occurrence rate of bursts per unit area which present echoing 
areas between A and A +dA to a radio wave of frequency », f(v) is some function 
of the frequency and c is a constant independent of the frequency. 

The occurrence rate of bursts which present echoing areas greater than some 
value A now becomes A~!cf(v), and (7) takes the form: 

Number of bursts in the range interval R, 
R, to R, observed to occur by a station ce fr) OCR, Od Bevee. 6 (11) 
operating on the radio-frequency v: ose 

Following this simple treatment the unknown function f(y) is clearly pro- 
portional to the value of the ratio m given in the final column of table 1, and it 
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would follow that stations operating on approximately the same radio-frequency 
should yield the same m. Examination of table 1 shows that the m-values for 
the Class A stations are substantially of the same magnitude and that f(v) for 
these stations will also be approximately constant since the frequency range 1s 
small. This correspondence suggests that the differences in station activities 
are not due to local variations of burst activity but may be adequately explained 
in terms of the physical factors of position and type of aerial, frequency and 
power radiated. 

This question of local variation of burst activity was later examined in more 
detail by means of simultaneous observations made from two stations situated 
near Aberdeen and Weymouth respectively. It was not found that even this 
latitude separation produced any systematic activity differences. 

The striking dependence of observed burst activity upon the radio-frequency 
employed is shown by comparing the results for the Class B stations with the 
figures for the A group. It is seen that doubling the radio-frequency produces 
a marked decrease in activity and a corresponding decrease in the m-values. 
The results are too rough to permit accurate deductions to be drawn and more 
work is clearly needed, but comparison of the average values of the two groups 
-of stations show that /(v) varies as 1/v? and equation (10) becomes 

ndA =p A=" dA) ae eee (12) 

It is clear that this apparent variation of burst activity with radio-frequency 
must be a consequence of the different echoing areas presented by a burst to radio 
waves of different frequencies and equation (12) leads immediately to the conclusion 
that the echoing area which a given burst presents to a wave of frequency v varies 
inversely as v®. Using this result it may now be deduced from equation (1) 
that the amplitude of the echo signal which a given burst returns to a station, 
-other things being equal, will vary as the fourth power of the wavelength of the 
radio waves employed. It is interesting to note that this result is in agreement 
with Eckersley’s long wavelength observations on scatter clouds (1932). 


§8. GENERAL CONCLUSIONS WITH REFERENCE 
TO THE METEORIC THEORY 

Reference has already been made to the development of the theory associating 
the occurrence of transient bursts of ionization in the ionosphere with the incidence 
of meteors, and it is of interest to examine whether the present experimental 
results are in agreement with this theory. 

In a comprehensive account of meteor phenomena given by Hoffmeister 
(1937) a distinction is drawn between the shower type of meteor which originates 
within the solar system and pursues an elliptic orbit about the sun and the so- 
called sporadic or cosmic meteor that arrives from outer space. It is now known 
that a large majority of the meteors incident upon the earth are of the sporadic 
type and that definite diurnal and seasonal cycles are associated with their rate 
of occurrence at a given point upon the earth. On a typical night the rate of 
occurrence of meteors increases to a maximum in the early hours of the morning 
and then the activity falls off until dawn. Similarly the nocturnal activity 
follows a smooth seasonal cycle with its maximum in the autumn and minimum 
in the spring. 

These variations of optical meteor activity have been explained in terms of 
the change of position of the observing station relative to the vector which defines 
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the earth’s instantaneous orbital velocity. This vector is defined by its point of 
intersection with the celestial sphere, a point which is termed the meteoric apex. 
The maximum number of observable meteor impacts will occur when the hemi- 
sphere containing the station is the leading face of the earth in its orbital motion 
and when the meteoric apex lies on the station’s meridian. This occurs at 
0600 hours; also at 1800 hours there will be a minimum of impacts since the 
station then lies on the rearward face of the earth and the meteoric apex is at 
its lowest point. Other configurations will yield activities between these extreme 
values. Similarly, the total daily activity will change with the elevation of the 
diurna] path of the meteoric apex over the celestial sphere, and thus show a 
seasonal variation with maximum in the autumn. 

It will be immediately noticed that the description we have given of the 
diurnal and seasonal changes of ionospheric burst activity closely parallels the 
above account of the time variation of sporadic meteor activity. A maximum 
of activity in the early morning and a minimum about 1800 hours are features 
of the diurnal cycle required by the above theory of meteors and these same features. 
have been found characteristic of the daily curve of burst activity (figure 1). 
The form of the seasonal curve of burst activity given in figure 3 shows similar 
agreement with the above description of the seasonal variation of meteor activity. 
The amplitude of the early morning peak also shows a seasonal change in con- 
formity with the changing angle of elevation of the meteoric apex, and in figure 1 
it is seen that the autumn and winter amplitudes exceed those for winter and. 
summer. 

The late morning peak of the burst cycle, however, cannot be explained in 
terms of known characteristics of meteors, so that some other approach is necessary. 
Appleton and Naismith (1947) have recently found a correlation between the 
noon level of burst activity and the state of ionization of the lower E region, 
and they suggest that the susceptibility of the ionosphere to impact ionization 
by a meteoric particle probably varies with the existing intensity of ionization. 
The ionized state of the E region is subject to solar control in a manner which is 
fully understood, so that, following Appleton and Naismith, it now appears 
reasonable to attribute the second maximum of the daily cycle of burst activity 
to the changing sensitivity of the lower E region. From this it would immediately 
follow that the amplitude of this peak should vary seasonally and be a maximum 
at the summer solstice; figure 1 shows that this is thegcase. 

It thus appears that the diurnal and seasonal cycles of burst activity may be 
explained in terms of two main factors, (a) the rate of incidence of meteors upon 
the area of the ionosphere visible from the station, (4) sensitivity of the burst- 
producing layer to impact ionization. Effect (a) is measured by the orientation 
of the earth’s orbital motion vector relative to the station, while effect (b) depends 
upon the state of ionization of the E region which is controlled by the elevation 
of the sun. Both these effects produce an afternoon fall in activity and so must 
result in the pronounced minimum which we have observed to be so consistently 
present on our daily curves. 

A further point of correspondence with the meteor theory is that of the 
heights of occurrence of the bursts. Olivier (1925) and Hoffmeister (1937) 
have collected the results of a number of observers on the mean heights of 
appearance and disappearance of meteors; the former value varies between 
160 and 90km., and disappearance occurs below 120km. and mostly at about 
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80km. In particular, Trowbridge (1907) has carefully studied meteor trains 
and finds that they tend to occur at a mean height of 87km. ‘These various 
figures may be compared with our result that the bursts appear to produce their 
scattered radiation from within a thin layer of the ionosphere at a height of 86 km. 

We therefore conclude that the evidence of the present investigation supports 
the theory that a large proportion of ionospheric bursts may be attributed to 
impact ionization produced by sporadic meteors. The absence of an eclipse 
effect, moreover, is not at variance with this view, for the eclipse of the solar 
radiation would cause only a temporary and slight reduction in the sensitivity 
of the critical layer at a time coincident with the normal fall in activity, so that 
the detection of such a slight change would be difficult by the burst-counting 
method. The daily variation in the rate of incidence of meteors would also 
obscure the small effect that might be expected. 


ACKNOWLEDGMENTS 
This work was carried out under R.A.F. auspices by the permission of Air 
Ministry and the A.O.C. of No. 60 Group; it was initially prompted by Sir 
Edward V. Appleton, F.R.S., to whom our thanks are due for helpful discussion 
and advice. We would also like to acknowledge the valuable assistance received 
from many R.A.F. personnel and particularly from F/O. N. R. Benzie. 


REFERENCES 

_APPLETON, E. V., and NaismiTH, R., 1935, Proc. Roy. Soc. A, 150, 685; 1940, Proc. Phys. 
Soc., 52, 402 ; 1946, Nature, Lond., 158, 936 ; 1947, Proc. Phys. Soc., 59, 461. 

APPLETON, E. V., NAISMITH, R., and INcRaM, L. J., 1937, Philos. Trans. A, 236, 191 

APPLETON, E. V., and PrppincTon, J. H., 1938, Proc. Roy. Soc. A, 164, 467. 

-EckerSLEY, T. L., 1932, 7. Instn. Elect. Engrs., 71, 405 ; 1937, Nature, Lond., 140, 846 : 
1940, f. Instn. Elect. Engrs., 86, 548. 

EcKersLey, T. L., and Farmer, F. T., 1945, Proc. Roy. Soc. A, 184, 1096. 

Hey, J. S., and Stewart, G. S., 1946, Nature, Lond., 158, 481. 

HoFrMEIsTER, C., 1937, Die Meteore (Berlin: Springer). 

Nacaoka, H., 1929, Proc. Imp. Acad. Tokyo, 6, 5. 

OuivierR, C. P., 1925, Meteors (Baltimore : Williams and Wilkins). 

Pierce, J. A., 1938, Proc. Inst. Radio Engrs., 26, 892. 

ScuaFer, J. P., and GoopaLi, W. M., 1932, Proc. Inst. Radio Engrs., 20, 1941. 

SKELLETT, A. M., 1932, Proc. Inst. Radio Engrs., 20, 1933 ; 1935, Ibid., 23, 132 ; 1938, 
Nature, Lond., 141, 472. 

‘TROWBRIDGE, C. C., 1907, Astrophys. F., 26, 95. 


Measurements of the Interaction of Radio Waves 
in the Ionosphere 


BY eG wi: HUXLEY, H. G. FOSTER AND-C, CaNEW TON 


Department of Electrical Engineering, University of Birmingham 
MS. received 28 September 1947 


ABSTRACT. ‘This paper describes measurements of radio wave interaction, in special 
test transmissions, between the B.B.C. transmitters at Droitwich and Lisnagarvey 
(N. Ireland). It is shown that the theoretical formula M=M,/[1+(n/Gv)2]3 accurately 
describes the dependence of the interaction modulation on the modulation frequency (2/27). 
‘The electronic collisional frequency v at the seat of interaction is deduced from Gy and 
the laboratory value of G. Investigations of the phase of the modulation are also described. 
The paper contains a sketch of the theory of wave interaction. There is also an adden- 
-dum dealing with the possible use of wave interaction as a tool in ionospheric research. 
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§1. INTRODUCTION 
HIS paper describes the results of a number of observations of the pheno- 
menon of non-resonant radio wave interaction in the ionosphere in special 
test transmissions. . "The equipment employed was similar to, but an 
improvement on, that used by us in earlier observations on gyro-interaction 
that have already been reported (Huxley, Foster and Newton 1947). The 
test transmissions which concern us in this paper were arranged by Mr. J. A. 
Ratcliffe of the Cavendish Laboratory in consultation with the Engineering 
Research Department of the British Broadcating Corporation, and through 
their courtesy we were informed of these transmissions. Mr. Ratcliffe also 
informed us of his method for roughly locating the region of interaction described 
in §5 below. 
A special efiort was made in our observations to obtain accurate measure- 
ments of the depth of interaction modulation as a function of audio-frequency. 
The transmitters used in the tests were Droitwich (170kw., 200kc/s., 80°, 
modulation) and Lisnagarvey, N. Ireland (60kw., 1050kc/s., plain carrier), 
although Droitwich was replaced by Ottringham, near Hull (170 kw., 167 kc/s.), 
on a few occasions. ‘lhe point of observation was in Selly Oak, Birmingham. 
Their relative positions are shown in figure 1. 
Before proceeding to a description of the experimental results, a brief account 
of the phenomenon of wave interaction and the elements of the theory are given in 
order that the significance of the observations may be the better appreciated. 


§2. DESCRIPTION OF THE PHENOMENON AND NOMENCLATURE 

Consider two radio transmitters; the one, which we shall call the interfering 
transmitter, radiates at high power a modulated carrier, whereas the carrier of the 
other, called the wanted transmitter, is unmodulated. An observer, so placed 
that he can receive the components of the transmission from the wanted station 
that have passed through the E region of the ionosphere in the vicinity of the 
interfering station, finds that, in general, the modulation of the interfering trans- 
mission is impressed upon the wanted transmission in the course of its passage 
through the E region. ‘This phenomenon is called Wave Interaction (previously 
the Luxembourg Effect). The phenomenon was first reported by Tellegen (1933) 
and the first satisfactory attempt at a theoretical explanation was made by Bailey 
and Martyn (1934). Later Bailey (1937) revised the original theory and gave 
a more rigorous and complete treatment. Systematic investigations of -the 
phenomenon were made by Grosskopf (1938). 


§3. SIMPLIFIED THEORY OF WAVE INTERACTION 
- The basic ideas of the theory given by Bailey and Martyn and later elaborated 
by Bailey are the following. The absorption coefficient ky of the ionosphere for 
the wanted wave whose frequency is p/27 at a place where the collisional frequency 
is v is known to be (neglecting the earth’s magnetic field) 


27e*?cN v 
ky= aie Spee SS (1) 
m Vp 
where e is the electronic charge in E.M.U., m the mass of the electron, c=3 x 101° 


and N is the electronic concentration incm~*. ‘The absorption coefficient is thus 
a function of the electronic collisional frequency and the frequency of the wanted 


wave. 
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The E region in the vicinity of the interfering station (frequency p,/27) is 
irradiated by the electric field of the latter, 

E,(1+D cosnt)icospjtaee. 0) eee (2) 
where D is the modulation depth and n/27 the modulation frequency. Suppose 
first that the modulation D is zero. The field £, cosp,t then supplies power W 
to an electron at a specified region in the ionosphere. When W is zero, the 
kinetic energy of agitation of an electron is the same as the thermal energy of 
agitation O, of the molecules of the gas. When, however, power W is supplied to 
an electron at a steady rate, a new condition of dynamical equilibrium results in 
which the mean kinetic energy Q of the electron exceeds that of the gas molecules, 
so that we may write O =kyQ;, where ky is Townsend’s energy factor. 

In this condition, let the collisional frequency of the electron be v and yO the 
mean energy lost per collision with gas molecules. ‘Then 


WO=W.. .e.8 oye (3) 
When the collisions between electrons and gas molecules resemble those between 
smooth rigid spheres, it is known that with a Maxwellian distribution of velocities 
y =2:66(m/M’')(1 —1/kp) =2:-66(m/M’)(1 —v?2/v?), 
where M’ is the mass of a molecule of the gas and 1, is the collisional frequency of 
the electron when in thermal equilibrium with the gas (W=0). 

Bailey has supposed that a similar formula is applicable to an electronic motion 
in a diatomic gas when 2:66m/M’ is replaced by a factor G which is characteristic 
of the gas. ‘Thus, 

| y= Gil-1h)=Gl—rie)2 = eee (4) 
In general, G greatly exceeds 2:66m/M’, as for instance in air, where 
2:66) Vs lO but. G-= 10 

When the coefficient of modulation D in equation (2) is zero, suppose FE, to 

supply power to an electron at the mean rate 

WARS ee ee (5) 
where A depends on the frequency p,/27, the collisional frequency v and e2/m 
(Appendix). 

When D is not zero the power supplied to an electron becomes time-dependent 
as follows: 

W(t)=W(1+ Dcosnt)?=W+2W(Dcosnt+4D? cos2nt), ...... (6) 
where W(t. D2) 0s & Gees eee eee (7 
The energy of agitation Q now becomes a cic of the time which we determine. 

Let Q be the mean value of Q when power Wis communicated to each el ectron, 
and Q, its value when W=0 (thermal equilibrium). It follows from (4) and (6) 
with kp= 0/O,= = (0/0) that 


dO Idi (Gr (OO) = Vt). eee eee (8) 
Similarly, (Gv)\(O —Q:) = WW D7/2) at ee ee (9) 
Hence dQ/dt + (Gv)(Q —Q)=2W[D cosnt+4D? cos2nt], ...... (10) 
2W {Dcos(nt--¢) D*cos(2nt —¢4’) 
hen ie SARIS 8 5 NEN | 
whence 9-8 [14 eGo Y ATOM [> OD 


where tan d =(n/Gv); tan di (2n/GP). we ee (12) 
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Thus the agitational energy Q fluctuates about its mean value O at a frequency 
(n/27). Since the collisional frequencies v and v corresponding to O and QO satisfy 


the relation (v/v)? =(Q/O), and since in practice the coefficient 2W, [O(Gr)] <1 
in (11), it follows that v is modulated as follows: 


qe E feels ia poset) Ge Cre amit | eet (13) 
Q(Gv) (1+ (m/Gv)’)? — 4(1 + (2n/Gr)?)3 

Thus the ratio of the amplitude of the second harmonic to that of the fundamental 

in the modulation of Q and » is 


L Din? + (Gy)2}/ [4n?-(Gr), 1 a (14) 
It is the custom to neglect the second harmonic but its existence in the correct 
proportion required by (14) has been established in measurements of wave 
interaction (Huxley, Foster and Newton 1947). 


‘In what follows we consider the fundamental term alone and write (13) 
v=v+y, with 


vWD_ cos(nt —¢) (15) 
el = eeeeae TPT STO Te bce ana eae : 
O(G7) [1 + (n/Gv)?]? 
It follows from (9) and (14) that the coefficient WD/O(G7v) in (13) may be written 
in the alternative forms 


WD = RS ele = pee) aS eed (16) 
OG) W+O(Gn [C+D 72)40(Gn/wy) 
and WDHO(Gr)=Db = Wid aD De. ae (17) 


When the wanted wave traverses a region in which the collisional frequency is 
modulated as in (15) by the field of the interfering station, the absorption coefficient 
(equation (1)) also becomes modulated as follows: 


k=ko +, dko/dv. 


Suppose the modulation of & to occur within a region which is traversed by the 
wanted wave alonga ray of lengths. "The emergent amplitude of the wanted wave 
in terms of its amplitude on entering the region is 


H= Eexp ( mal kds) =, exp ( — f. ie as) exp ( — I, v,,(dRo/dv) as) ; 
Fees (19) 


Since the total attenuation of the wave is in practice not excessive, the exponent 


8 s , : 
-| k, ds is of the order of magnitude of unity or less. ‘The increment in the 
0 


s 
absorption is still less and we may write | v,(dky/dv) ds <1; consequently (19) may 
0 


8 s 
par exp ( - i he as) E k: { v, (dk, |dv) as. cues (20) 
0 0 


Thus, the depth of the impressed modulation, expressed as a fraction of the 
received amplitude of the wanted wave is 


be written, 


ae | HiGid cee eT. (21) 
0 
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Suppose for simplicity that ky and y are constant along the path of integration 
through the region of interaction. This implies either that the path of integration 
lies at the top of the trajectory of the ray, or that, if the interaction region lies on the 
incident or reflected ray, it is effectively localized. The expression for M then 
simplifies (from (16)) to 

Dv(dko/dv)s cos (nt —¢) (22) 
[1+ D2/2+(Gv)O./W] [L+(n/Goy}e 

Which alternative is the more appropriate is to be decided from experience. 
In any case, it is supposed that where interaction occurs p?>v®, so that according 
to (1) the absorption coefficient is proportional to v. The total absorption of 


M= 


§ - . 
the wanted wave in nepers is | k, ds, which is, therefore, of the form Cv where C is 
0 


determined by medium and the path. The term s(dkp/dvo)v in (22) may be 
replaced by (Cvs), which can be estimated from the reflection coefficient of the 
E region for the wanted wave. 

The expression (22) for M shows that in this case the dependence of MV on the 
modulation frequency is relatively simple and of the form 


MMHG ee (23) 


where, from (22), M,y=DCsv/[(1 + D?/2) + O,(Gv/W)). 
Further, the interaction modulation lags in phase with respect to the inter- 
fering field £, at the region of modulation by the angle 


p=ian(n/Gr at ee) eee (24) 


In the measurements described below special attention was paid to the investi- 
gation of Masa function of m and it was found that formula (23) describes accurately 
what is observed when n/27 does not exceed 750c/s. From such measurements 
the quantity Gv may be deduced in each instance. 

In practice, the term (Gv)Q,/W>(1+D?/2) in (22) so that M, is effectively 
equal to WDCsv/O,(Gy). 


§4. MEASUREMENT OF INTERACTION MODULATION AS A 
FUNCTION OF MODULATION FREQUENCY 


(1) Calibration of equipment 


Because the depths of interaction modulation are relatively small in practice, 
it was necessary to use measuring equipment capable of giving reliable results 
down to modulation depths of about 0-1°% on a carrier received from a relatively 
distant (300-400 km.) medium wave station. It was decided therefore to use not 
an oscilloscope but a specially calibrated receiver and to obtain the carrier voltage | 
and the audio-frequency voltage at the detector from meter readings, from which 
the modulation depth could be obtained from calibration curves. The receiver 
was a modified commercial communications receiver which covered the required 
range of wavelengths. ‘The detector stage was modified to feed directly into a 
valve voltmeter circuit which was built into the receiver chassis, and which was 
calibrated to provide the carrier level at the detector in the form of a meter reading. 

The audio frequency voltages from the detector were fed, through suitable 
filter circuits, to a cathode follower, the output of which was connected directly toa 
Marconi wave analyser which was available in the laboratory. The advantage 
of a wave analyser lies in its narrow bandwidth of only 4c/s. which enables the 
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wanted modulation tone from the detector to be measured relatively free from 
background noise and interference. 

Through a laboratory calibration it was possible to use the direct readings of the 
wave analyser in combination with those of the carrier-level meter mentioned 
above to obtain the required modulation depth immediately. The wave analyser 
had, in fact, been originally used to measure the second harmonic of the interaction 
modulation shown in equation (13) as mentioned in the paragraph that follows it. 

The calibration of the system to read modulation depths direct was achieved 
by applying to the input, a radio frequency carrier whose frequency was near that 
of the wanted transmission. This carrier was modulated by a series of pure tones 
covering the range of frequencies of 50 to 2000c/s. The modulation depths, for 
the purposes of calibration were set in turn at 24°, 5°/ and 10°. Later this 
calibration was found to agree with a calibration tone of known depth specially 
imposed upon the carrier of the wanted transmission radiated by Lisnagarvey 
(and Westerglen on another occasion). 

The laboratory calibration unit comprised a crystal oscillator driving a power 
amplifier which was modulated by the appropriate tone. 

By extrapolation of the calibration curves which are straight lines, modulation 
depths down to 0:1°% could be confidently assessed. 
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Figure 1. Gnomonic projection centred on B, showing stations and point of observation. 


(ii) Experimental results 


The transmitters concerned are shown in figure 1, in which D, L and O 
represent respectively Droitwich, Lisnagarvey and Ottringham. ‘The point of 
observation was in Birmingham and is indicated by B. ‘The distance LB is 
340 km., whereas DB is 20km. and the angle LBD is approximately a right angle. 

The transmitter at D radiated 170 kw. in the form of a carrier on 200kc’s. 
modulated by a series of pure tones to a depth of 80°, while at the same time the 
transmitter L radiated 60 kw. as a plain carrier on a frequency of 1050kc’s. 


10-2 
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The times of observation were between the hours of 2230-0300 c.M.T. on the 
nights of June 15/16 and 16/17, July 3/4, 4/5, 5/6 and 6/7, and August 6/7, 7/8 and 
8/9, 1947. 

At intervals during the night, pulse transmissions were sent out by Lisnagarvey 
which showed that the ground wave received in Birmingham was relatively weak, 
and that the predominant wave received during the early night was a single 
reflection from the E region, often accompanied later by other forms of reflected 
waves from the E or F regions with relative strengths varying continually. 

In the test runs, in which M was determined as a function of (7/27) the modu- 
lation frequency, Droitwich was modulated in turn, for two minutes each, at the 
frequencies 50, 75, 100, 150, 200, 300, 400, 500, 750, 1000, 1250, 1750 and 2000 c/s. 

The modulation depths were found from simultaneous readings of the carrier 
level meter and the wave-analyser as described in the previous section, using 
A.v.c. to hold the carrier level constant at the detector. It was confirmed that the 
use of a.v.c. did not affect the measured modulation depth, except to improve the 
accuracy. 

Two examples of the results of such measurements are shown in the curves of 
figures 2(a) and (6). They refer to the times stated on the figures. 


0 200 400 600 800 1000 0 100 200 300 400 500 
Modulation Frequency Modulation Frequency 
Figure 2(a). Interaction modulation Figure 2(b). Interaction modulation 
depth vs modulation frequency : depth vs modulation frequency 
—Droitwich—Lisnagarvey 0230- —Droitwich—Lisnagarvey 2330— 
0300 G.m.T. 16/6/47. 2400 a.m.T. 3/7/47. 


If the theoretical formula (23), 

M= Mail (n/ Gre” Se oe eee (25) 
actually describes the dependence of M on n, then M and n should also be related 
as follows :— . 

n= Mi GuiM ye =(Gryr, os es (26) 
in which M, is the value of M as n+0. 

Thus, the curve which represents n? as a function of 1/M? is a straight line 
whose intercept on the n? axis is (Gv)?. In figures 3 (a) and (b) the experimental 
data exhibited in figures 2(a) and (6) are plotted with n? as a function of 1/M/?. 
It can be seen that for frequencies n/27 from 50 to 500 c/s. a linear relation subsists ; 
consequently, at the lower frequencies n/27, where the measurements are more 
accurate, equation (25), does in fact closely describe the experimental facts. 

Evidently the quantity Gy may be found from the intercept on the n? axis of the 
straight line through the experimental points. Further, M, is found from the 
intercept on the 1//M? axis: . 
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It is therefore necessary only to give the zero frequency modulation M, and 


i. value of Gy in order to specify completely the modulation M at any frequency 
n|2r. 


100 


80 


60 


(/aq? x 10° 
Figure 3. Points in figure 2 replotted to exhibit 1//? as a function of n?=(37f)?. 


In what follows the modulation is specified in this manner. 
Tables 1 and 2 summarize the experimental results on the dependence of 
Monn. The values of Gy ranged from 1100 to more than 2000 in these tests. 


Table 1. Observed depths on interaction modulation. Wanted station— 
Lisnagarvey (1050 kce/s.) 
Interfering transmitter—Droitwich (200 ke/s.; 170 kw.) 80°%. 


Date Time cm.t. 1, (%) Gv Date Time c.M.T.. My (%) Gv 
June 16 0230-0300 0-6 1400 July 5 2230-2300 1-63 2320 
ee lor 92300=2330 ive7/5) 1480 nape 2330-2400 2S — 
pet (an 0230-0300 0-9. irregular ee 2230-2300 Ahoy 2860 
July 3 2330-2400 Ilo7S 1950 ay 0 8) 2330-2400 1-4 1480 
a 4 2230-2300 175 1870 ~» 2 Olson Del 1480 


4 2330-2400 1°85 1520 
Interfering transmitter—Ottringham (167 kc/s.; 200 kw.) 80%. 
Date Timec.m.T. My (%) Gp 
June 16 0200-0230 0-58 1550 
» 17 0200-0230 0-45 1480 


Table 2. Observed depths of interaction modulation. Wanted station— 


Lisnagarvey (1050 kc/s.) 


Interfering transmitter—Droitwich (200 kc/s.; 170 kw.). ; 
Date Time c.m.T. My (%) Gr h(km.)| Date Time c.m.7. My) (%) Gr h(km.) 


Aug. 6 2300-2328 2°6 1550 75 |Aug. 7 2300-2328 2°6 1340 95 
» 6 2340-2348 2°6 1100 = 68 ed e223 45 2-1 — 86 
», 7 0010-0018 f-552° 1050 ~ 70 ,, 8 0010-0018 1:9 1720 = 86 
», 7 0040-0048 Dee 1570 = 88 », 8 0040-0048 225° 1870 86 
EEO TI0-0155 1-9 1450 = 88 i non OlSO=015s Deh 1207, 
5, 7 0210-0218 —_ — 83 », 8 0210-0218 ISS 1790 = 86 
> 4 0240-0248 2:0 1950 = 83 », 8 0240-0248 iles — 86 

9 0100-0128 0°77 1140 82 


h=height of interaction. 
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It was found in an isolated observation of wave interaction between Westerglen 
(767 ke/s.) near Falkirk and Ottringham (167 kc/s.-170 kw.) on 8th May (0030— 
0100) that M,=4°% and Gy=1800. Thus, Gy is of the order, in these experi- 
ments, of 1-5 x 10%. If it be assumed that the temperature and composition of the 
air at the region of interaction are approximately the same as those at the ground, 
then, as Bailey has done, we may suppose that G=2:6 x 10-3, deduced from the 
experiments of Townsend and Tizard on the diffusion of slow electrons in air. 
Then from the experimental value of Gv we find that v is approximately 6 x 10°, 
and exceeds 1 by a few per cent. 

Recent experiments to determine G by Mr. Zaazou and one of us in this depart- 
ment have shown that its value for air is 1-3 x 10-3 and not 2-6 x 10% as supposed 
hitherto. ‘Thus » becomes 1:2 10°. Further, if oxygen is present as atomic 
oxygen at the interaction region, then G should be somewhat smaller than the 
laboratory value. 


§5. THE PHASE OF THE INTERACTION MODULATION 

It was pointed out by Ratcliffe that an examination of the phase of the inter- 
action modulation should provide a further test of the theory. According to 
equations (12), the phase of the modulation of the collision frequency lags on that 
of the interfering field E, at the place of interaction by an angle ¢=tan-!(n/Gy). 
It follows therefore, from (22) that the phase of the interaction modulation, with 
that of interfering field as the standard of reference as above, lags by 
@=tan!(n/Gvy) or tan+(n/Gvo)+180° according as (dk,/dv) in (21) has a 
negative or a positive sign. Since, from (1), kj is a maximum when p=», and 
(dk/dv) is negative when v>p, but positive when v<g, it follows that the very 
low frequency modulation M,(n—> 0) should vibrate in phase with the interfering 
modulation if interaction occurs below the height where v=/p, but should vibrate 
in anti-phase if injected in a region at a greater height. 

It was found that the zero frequency modulation was anti-phased with respect to 
that of the interfering field. It may therefore be concluded that interaction takes 
place above the level at which v = p(where p/27 is the frequency of the wanted wave). 

In figure 4, D represents 
the interfering transmitter 
(Droitwich), B the point of —oL 
observation (Selly Oak, Bir-  (Usnagarvey 
mingham) and I the centre of 
a region where interaction 
occurs, ~~ hus), 1B: is 7a uray 
from the wanted - station 
(Lisnagarvey) that reaches B. 
Let DB=/, (20km.) and put 
l,=D1+IB. Take as refer- 
ence phase for the modulation ‘Ds 
tones that of the interfering. Figure 4 me ee ; ; 
modulation at D, its source. gure 4. Relative positions of the wanted and interfering 

, stations and of the reception point. . 
The phase lag #5, of this modu- 
lation when received at B on the carrier of the interfering station is 


ib, = 300/E ie degreccsemnaetey 9 eee (27) 


where the modulation frequency f=/27c/s. and v=3 x 10° km/sec. 


Region of 
Interaction 
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Similarly, the phase of the interaction modulation generated at I, when 
received at B on the carrier of the wanted wave, lags on the reference phase at D by 


bo = 360fl,/v + tan (27f/Gv) +180 degrees. —....... (28) 


Thus, if the audio outputs of identical receivers are compared on a double-beam 
oscilloscope, the relative phase lag at B of the interaction modulation with respect 
to the interfering modulation, which is 


(t. — 1) = 3 60f(/, —1,)/v + tan (2nf/Gv) +180 degrees, ...... (29) 


may be measured. If the receivers are not identical the phase shift introduced by 
each at any frequency f is easily determined and the difference used to correct the 
value of (¥, —y,) given by the oscilloscope. 

Let #% be the frequency-dependent part of (29), that is 


p = 360f(, —4)/v + tan! (27f/Gy,). yee SU) 
It follows that s=0 when f=0, and that 
p— 360f(2, —1,)/v+90°, as f> o. 


Thus the curve representing % as a function of f passes through the origin and 
possesses an asymptote whose slope is 
povl@yid] \=300(=L)/v, ee ees (31) 
and whose intercept with the y-axis is 90°. 
When the measured values of # are plotted against f, the distance /, may be 
obtained from the slope of the asymptote as follows: 


ee a, Ose, ee eee ee (32) 

It follows that if the elevation of the downcoming ray IB and the angle LBD are 
known, then the point I may be determined. ‘This method of fixing I from (30) 
and (32) was pointed out to us by J. A. Ratcliffe. 

Figure 5 shows that a 
_ typical experimental curve 
representing ¢ as a function 
of f conforms closely to 
what is to be expected from ai 
(30). On the assumption 
that the region of inter- 
action I lies on a ray 
reflected at a height of 
100km., it was found that 
in most cases the apparent 
interaction point obtained 
from (32) was located at a 
height of (86+5)km. The 299 
method therefore provides 
a rough estimate only of 


1000 


E LAYER 


Phase Delay (degrees) 


the height. 0 

The most stable read- ‘ or Medulstion Faauaner (o/s) ee 
ings of amplitude and phase Bicure.st 
are obtained when one of Note. Inset: Point 1 obtained from asymptote 1. 


the waves reflected from the Point 2 obtained from asymptote 2. 


144 L. G. H. Huxley, H. G. Foster and C. C. Newton 


ionosphere is predominant. In these tests the predominant wave was 
the wave undergoing a single reflection at the E region as was shown 
by pulse transmissions from Lisnagarvey. However, except during the early 
evening, other waves reflected fromthe F or E regions appeared in varying strengths, 
and the interaction modulation on these components was able, on occasions, to 
shift the phase in an irregular fashion and to alter the amplitude of the interaction 
as measured at the point of observation. 

It was noted that when a.v.c. was employed in measurements of modulation 
depth, the audio-frequency amplitude showed sudden departures from a steady 
value, although the carrier was held fixed. The steady level was assumed to 
correspond to the single E wave and the fluctuations to the presence of other sky 
waves. 


§6. REMARK ON THE MAGNITUDE OF THE INTERACTION 


The maximum depths of interaction modulation M) measured in these tests 
were small, as can be seen from tables 1 and 2, and we believe that this relatively 
small value is to be attributed to the relatively high frequency of the wanted 
station (1050 kc/s.). 

It is of interest to compare the values of M, listed in tables 1 and 2 with those 
obtained in other transmissions. . 

Thus, van der Pol and van der Mark found that the interaction M, impressed 
on Beromunster (556kc/s.) by Luxembourg (150 kw.) was about 8°, and as 
mentioned below table 1, we observed that Ottringham (100 kw.) impressed upon 
Westerglen (767 kc/s.) a modulation M, of about 4%. 

In these tests the modulations M, impressed upon Lisnagarvey (1050 kc/s.) 
by Droitwich (170 kw.) were less than 2%. ‘Thus, M, appears to diminish as p 
increases. 

ACKNOWLEDGMENT 

The authors wish to express their appreciation of the co-operation and interest 

shown by Mr. H. L. Kirke and Mr. R. A. Rowden of the Engineering Research 


Section of the British Broadcasting Corporation in providing them with schedules 
of the transmissions. . 


ADDENDUM 
Suggested measurements of wave interaction 
By L. G. H.sHUXERY. 


It is explained, in the above paper, how a determination of the amplitude and 
phase of the impressed modulation M as a function of modulation frequency 
n/27 leads to value of the product Gy. It is natural to enquire whether the pheno- 
menon of interaction of radio waves is able to provide estimates of the value of other 
molecular and electronic parameters of the ionosphere. To this end, consider the 
general formula (22) for the impressed modulation. According to this formula, 
and what immediately follows it, the amplitude of the zero frequency modulation 
M, is 

DCvs 
- (Danone. ce) 


The mean carrier power W supplied by the interfering field E, to an electron 
is proportional to the radiated power P of the interfering transmitter, so that in 
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general M, is not proportional to the radiated power P, but becomes so in effect 
when the term (Gv)O,/W>1. 

We therefore first consider the order of magnitude of (Gv)Q,/W in the practical 
case where P=200kw. and the frequency p,/27=2x 105. We may consider, 
with good approximation to the truth, that the aerial with its image radiates as a 
Hertz doublet. 

The amplitude £, of the field produced at a point at height # in the E layer with 
an angle of elevation « with respect to the transmitter is given by 

Ey =45 x 10-*(P sin 20)/h? (voltsiem)?, aves: (A 2) 
where P is in kw., # in km. and F, in v/cm. 

It follows that the amplitude (peak) of the field EF, when h=100km., and 
P=200kw. is 2:124/2 x 10-4 sin 2« v/cm. 

To estimate W, assume a collisional frequency »=10®. We may suppose 
therefore that pj=(47 x 10°)?>v?, and proceed to estimate W=AE?, with 
A= {te?v/m(pj + v?)} cos? 0 
We find W=1-8 x 10" cos? @ sin 2x erg/sec. If the temperature of the region 
where interaction occurs lies between 200° and 300°k., then Q, is approximately 

le cig ince Gre 15.10", we tind 
(Gv)Q,/W = 4-2/cos? @ sin 2«. 

The maximum modulation is obtained when the point of observation is so 
chosen, if possible, that 0, the angle between £, and the earth’s magnetic field, is 
zero and «is 45°. The term (Gv)Q,/W in (A 1) is then approximately (1 + D?/2) 
and the impressed modulation M, would be less at P=200kw. than would be 
predicted from its value when P is small, on the assumption that M) is proportional 
tOP: ; 

If, however, the maximum power P were notably less than 200 kw. or the angles 
@ and « appreciably different from zero and 45° respectively, then (Gv)Q,/W at 
maximum power would certainly exceed (1 + D?/2) by a factor of eight or more. 
An observer, unless his measurements of modulation M,) were very accurate, 
would judge M, to be proportional to P up to the maximum value of P. 

If, on the other hand, a very powerful transmitter capable of radiating 300 kw. 
or more, or a strong gyro-transmitter, were employed to produce the interfering 
field, then measurements of M, by the sensitive method described in §4.1 should 
reveal the lack of proportionality between M, and P under optimum geometrical 
conditions. 

Such measurements should also provide additional information about the 
physical properties of the E layer, as may be appreciated from the following 
considerations. 

For equation (A 2) write E7=KP, then, W=AE{=AKP=D5P, where 


e? v cos? 6 sin 2% 
Ob ee iy oat a oie (A 3) 
Write (A 1) in the form 
M,=a/{1 + D2/2+(GYQgbP], eee (A 4) 
where, from (1), a DCis eee ay (A 5) 


Expression (A 4) is equivalent to 


1/M,=(1+D?/2)/a+(Gr)QyabP. we (A 6) 
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If, therefore, the experimental values of 1/M, were plotted against 1/P. the 
points should fall near a line whose slope is (Gv)Q,/ab and whose intercept on the 
1/M, axis is (1 +.D?/2)/a A : 

Thus a and (Gv)Q;/5 could be found individually. Since Gy is obtained from 
measurements of M, as a function of modulation frequency 7/27, it follows that 
we obtain, from (A 6), a and Q,/b, and the absorption coefficient Csv froma. Hf, 
however, proportionality between M, and P persists at very large values of P, then 
it is possible to set lower limits only to Csv, the total absorption in nepers along 
the path s. 

Appendix 

An important quantity that appears in the theory of wave interaction as 
outlined in §3 is W, the mean power given to a free electron in a gas by a high 
frequency alternating electric field E, cos p,t, which for simplicity was written in 
equation (5) as W= AE}. Expressions for W were derived by Bailey (1937), 
and independently at about the same time, slightly more accurate expressions were 
found by one of us (Huxley 1937). (Note that Bailey’s w is 1/v times the W used 
in §3. It is the work communicated per collision.) It is found that the high 
frequency conductivity of a gas containing NV electrons per cm? is 


52 ae Pi ) A 
Claes (=a) (1+, EMU CO Oe aceoe (Ay) 
and " aia 1 E (p, —o)? 
2 ST “ye (Pia + (p, —w)? 
“PE top [t+ (Pie en a 


in which oy and o, are the real parts of the conductivities parallel and perp 
to the earth’s magnetic field H and w=He/m (E.M.U.). 

The work done by the components of Ey parallel and perpendicular to H is. 
W, = (40,2? cos? He and W, =(40,Eysin?6)/N, where @ is the angie between EF, 
and H. When p}<v?<w? as is the case in these tests, then 


2 Ne? v ( fi =); 
all 5 eae ii ee A9 
os pate ee pi + a 
and WeAR*. ©) 4 a5 Gene (A 10) 
: 26 e Vv pe 
: _ corte ‘at Semeeeaiia 
es pam a Gap) (1s PtP) eae 
The following simplifications of (A 11) are applicable with the restrictions. 
specified : (a) p? >v2: AZ 2e%v cos? 6/3mp?: 


(6) pi~v?: AZ ev cos? 6/2m(v? + p?); 
(c) pi<v: Az e®cos?4/3mp. 
In these special cases W is easily calculated from (A 10); the appropriate 
value is then used in (22). 
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The Diffraction Structure of the Elementary Coma Image 


By R. KINGSLAKE 
Eastman Kodak Company, Rochester, N.Y. * 


MS. received 25 October 1947 


ABSTRACT. The integrations reported in this paper demonstrate conclusively that the 
elementary coma image is of the size and shape determined by ordinary geometrical optics, 
and that the only effect of diffraction is to break up the image into an elaborate fine structure 
of dots and lines of light. Practically all the light in the image is confined within the tri- 
angular space between the principal ray at the tip of the figure and the sagittal focus. 
Between the sagittal and tangential foci there is a small amount of light broken up into a 
series of approximately concentric curved bands of darkness and light, centred about the 
brightest part of the image. Even when the amount of coma is very small, comparable to 
the Rayleigh Limit, the characteristic shape of a typical coma-image is already making its 
appearance. ‘The accuracy of the theoretical predictions is fully and completely borne out 
by the actual photographs of a comatic image. 


V1 ENE RODWIC TlON 
y the aberration ‘‘coma’”’ is meant that third-order or Seidel aberration 
which leads to a wavefront differing from the ideal sphere by an amount 
P, given by 


Eze Oy Vat degen get “Ate oea: (1) 


where y, z represent the coordinates of a point in the exit pupil of the lens, and 
a the amount of coma present; a varies linearly with the image height above 
the lens axis. The centre of curvature of the ideal or reference sphere is assumed 
to be the ordinary Gaussian extra-axial image point. In the absence of distortion, 
this is the point where the principal ray of the oblique pencil pierces the paraxial 
focal plane. The minus sign in (1) indicates that with positive coma, the wave 
is lagging for points above the middle of the lens. 

Since we are discussing one of the Seidel aberrations, we must assume that 
both the aperture and the obliquity are very small. ‘The longitudinal departure, 
P, of the wavefront in the exit pupil, relative to the reference sphere, amounts 
at most to a few wavelengths and is thus infinitesimal compared with y and sz. 


§2. GEOMETRICAL IMAGE 


The y and z components of the point where a given ray pierces the paraxial 
focal plane, relative to the Gaussian image-point, are given by : 


y = —l—; ef = —1— a) sie eee. (2) 


where J is the longitudinal distance from the exit pupil to the image plane. 


Thus 
wy" = al(3y? +2°); Cet ces | i ne (3) 


* The major part of the computations involved in this paper were made by the author in 1936-37 
when he was on the faculty of the Institute of Applied Optics, University of Rochester. 
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These quantities are best interpreted by expressing the aperture-point 
position in polar coordinates. Writing 
y=" cose: sz=rsinGe. 5 +) EAeerie (4) 
the components of the ray-point displacement become 
y’ =alr*(2 + cos 20) ; eat] gal (361W'40) ee cok oS (5) 
The values of y’ and 3’, for a single zone of the aperture, give a double circle 
of radius alr? in the image plane as shown in figure 1, the centre of the circle 


‘ 
Y 
° j 
o 
225)/45 
8 AS 
270 ; >7 
3alr* 
160 
7! 
GAUSSIAN IMAGE 
Figure 1. Ray-point in exit pupil and in Figure 2. The 
image plane. geometrical coma image. 


being at a height 2a/r? above the Gaussian image-point. Each elementary zone 
of the lens forms such a double circle, and the complete family of double circles 
for the circular lens aperture fit between two common envelope lines, at 60° 
to each other, which intersect at the Gaussian image or principal-ray point 
(figure 2). : ; 

The top and bottom (“‘ tangential’’) rays of any zone intersect at the top-most 
point of the double circle in figure 1, and the front and rear (‘“‘sagittal’’) rays 
of a zone meet at the lowest point of the double circle. If the height of the 
tangential image above the principal ray is comay, and the height of the sagittal 
image comag, then 


cOmay=3'comay, 5). yee (6) 


$32 DH PIAS BaR BEAD TONS HES sO LC @IvMEN 
The phase difference between light arriving at the Gaussian image-point 


from the point (y, 2) and from the centre of the exit pupil is given by equation (1), 
since phase=27P/A. If, however, 


we shift our reference-point to 
some point vertically above the 
Gaussian image, we add a further 
phase difference which can be used 
to simplify computation of the \ 

light distribution. From figure 3 \ 

it is seen that the phase difference Tata) aCe cane 
at D between light travelling along ee dea shir of the eee : : 
AD and along PD is given by : 


2m 2n[y : 
=| pp -ap |==]# : on radians. pee) 
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Hence from (1) and (7), the phase-difference between light arriving at any 
reference-point on the centre-line of the image from (y, z) and from the middle 
of the aperture is (277/A)[ —ay(y? +22) +(y/L)8h]. 

Since for any single zone in the exit-pupil aperture x is negligible compared 
with y or zg, this may be written 


(ZajA)\y/D| alr? -+8h|: © ee 6s... (8) 
The phase-difference will be zero for one entire zone if 
ONS Ores Ng ap = 6 4 AT IRN gees (9) 


and, since alr? is equal to the sagittal coma of the zone, we reach the very important 
conclusion that each separate zone of the lens forms a perfect equiphase focus at 
the sagittal image of the zone. Hence to determine the light distribution in the 
focal plane, for a complete lens, we divide the lens aperture into a number of 
circular zones of equal area, and integrate the amplitude contributions of all 
the zones at each of a number of discrete points in the focal plane. Contours 
of light intensity can then be drawn through the various points so computed. 


S45 TDHESAMPLITUDE DISTRIBUTION IN THE IMAGE OF A POINT 
FORMED BY A SINGLE NARROW CIRCULAR ZONE OF THE LENS 


This problem is classical, and a good treatment has been given by Martin 
(1930), who shows that the amplitude at a radial distance ¢ from the centre of 
the image of a single zone of radius 7 1s given by 


(7.4). eae (10) 


For convenience in computation and in plotting the results, we shall assume: 
that the radius of the lens aperture is unity. We can then introduce a symbol 
W =2rt/Al to represent the radial distance of our integration point R from the: 
centre of a zonal image. By comparing (10) with (7) it may be seen that the 
physical meaning of W is the phase difference between light reaching the point 
R from the edge and centre of the (unit) marginal zone of the lens aperture. 
In the case of coma, there is no path difference for any zone at its sagittal focus, 
and here W is zero. Thus the amplitude contribution of a zone of radius 7, 
at a radial distance W (radians) from the centre of the zonal image, is given by 


Upto! athlle Be Bee (11) 


§5. RADIAN MEASURE OF THE AMOUNT OF COMA IN THE LENS 


The amount of coma present in the lens may be expressed by a type of radian 
measure similar to that used in the W notation introduced in §4 by writing 


B2Or Nl) (Gl) 22naj autor ee (12) 


By comparing (1) with (12), it is clear that B represents the phase difference, 
in radians, between light reaching the Gaussian image-point from the top and 


centre of the lens aperture. 
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These two quantities, B and W, may therefore be used to replace comag 
and ¢ respectively, with complete assurance that the units will be properly 


compatible. 


§6.: THE RAYLEIGH LIMIT FOR COMA 
Conrady (1929) has stated that the Rayleigh Limit for Seidel coma is reached 
when 
cOmay = 05 A/sin Uap 2 ee ee eee (13) 


which becomes A//2 for a lens of unit radius. Converting this into the B notation, 
we find that : 


B-tolerance =(27/Al). (Al/2) =7 radians. 


Now Rayleigh based his limit on the fundamental idea that the aberration 
in a lens will be acceptable if it is just sufficient to cause a diminution of the 
peak intensity in the image to 80%, of its value for a perfect lens. It is therefore 
of interest to plot the variation of intensity along the centre-line of the image 
when B=7 to see how closely Conrady’s statement is borne out. 

The amplitude at a point on the centre-line of the image, distant QO (radian 
measure) from the Gaussian image-point, can be found by integrating the zonal 
contributions given by equation (11). The integration must be taken over 
a set of infinitesimal zones of equal area, i.e. of equalincrementin7?. In figure 4, 


Z Z 
shoe ees 
2 Se 
Br Sw \ 
eZ R Nes 
—— 
Guns icc 
{5 pans \ ) 
— Qs 
SK VA 
aE? 
Figure 4. Procedure for integrating the light Figure 7. The distance W from a zonal 
distribution along the centre-line of the focus to the point of integration, R. 


coma figure. 


the heavy dots represent the sagittal foci of such a series of zones, equally spaced 
along the mid-line since the 7? of the zones are equally spaced; G is the Gaussian 
image-point, S the sagittal focus of the marginal zone for which r=1, and R 
the point at which the net amplitude is to be determined by integration. Since 
the distance GS represents B, and GR represents Q, the distance of R from the 
equiphase focus of any zone r is W=(Q—Br*) radians. The total amplitude 
is thus, by equation (11), 


2 A 
| Fol Wr)d(r*)= | Jolr(Q — BF) dQ eu, 5) eee (14) 
‘This equation has been given by Steward (1926a, b) and Buxton (1926). 

A series of values of (14) were calculated for a succession of Q values, for 
B=7 and by plotting graphs of the J, values against 7? and integrating with a 
planimeter, the central section of the coma image could be plotted (figure 5). 
It is clear that unit intensity can be obtained only when B=0 and O=0F that 
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is, at the centre of a perfect Airy disc. It will be seen from figure 5 that Conrady’s 
tolerance is not in accordance with Rayleigh’s rule, for the maximum intensity 
is about 89%, of the central Airy disc intensity when B=z. 

Steward’s (1926a) data for a central section of intensity of the coma image 
for the case B=4 are shown in figure 6. This case does actually represent the 
Rayleigh limit, for the peak intensity is now 80%. Conrady’s formula (13) 
above should therefore be emended to read 


Rayleigh limit for comag = 0-64 A/sin Uy. 
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Figure 5. ‘The centre-line section Figure 6. The centre-line section of the coma 

of the coma image, for B=7z. image, for B=4 (from Steward). 


§7. LIGHT CONTOURS 
Light contours for the case B=4 


To determine the light distribution over the whole image, we assume that 
the lens aperture has unit radius, and divide the aperture into a number of zones 
equispaced in 7?. ‘The procedure follows that suggested by Conrady (1919), 
except that he used Simpson’s rule for the mechanical integration instead of 
a planimeter. A point R is taken in the focal plane (figure 7), and the oblique 
distances W from R to each of the equiphase foci of the various zones are tabulated. 
J,(Wr) is then plotted against 7? and the area of the graph determined in the 
usual way with a planimeter. ‘This area is the net amplitude at the point R 
due to the superposition of the contributions of all the zones in the entire lens. 
By taking a succession of points R in the image plane, a number of longitudinal 
and cross sections of the light distribution (in amplitude) can be plotted. The 
points where these sections cross the desired contour lines of intensity can then 
be joined with smooth curves. 

For the case of B=4, the lens was divided into eight zones having radii r 
equal to \/(1/8), »/(2/8), »/(3/8),. . . 1-0 respectively. ‘The amplitude-contribu- 
tion graphs were plotted for each of about 150 points in the image plane, having 
Y values lying between —10 and +14 (radian measure) and Z values lying 
between 0 and 8. As the pattern is symmetrical about the Y axis, only one side 
need be computed. Each curve was then integrated with a planimeter giving 
directly the net amplitude at that particular point. The results of the integrations 
are summarized in table 1. 
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Table 1. Computed amplitudes for B=4 


Yh, 0) +1 +2 Fass) aes +5 +6 +8 

Pe 

—10 0 

— 9 —0-008 

— 8 —0-002 —0-006 

407 0-002 

= 6 0-020 0-024 — 0-004 —0-002 —0-:012 

=) 5 0-044 0-038 0-028 0-014 0-0 

a= ot 0-064 0-060 0-070 0-048 0-040 0-016 

— 3 0-078 0-078 0-070 0-052 0-048 0-052 0-036 

— 2 0-104 0-088 0:066 0-038 0-036 0-040 0-052 

— 1 0-218 0-186 0-124 0-052 0-004 0-008 0-034 
0 0-426 0-368 0-218 0:078 —0-016 —0-046 0-014 
1 0-648 0-580 0-376 0-134 —0-028 —0-086 —0-032 
2, 0-848 0-734 0-474 0-168 —0-048 —0-118 —0-070 0-060: 
3 0-888 0-772 0-472 0-158 —0-078 —0-156 —0-110 0-064 
4 0-710 0-584 0-342 0-052 —0-142 —0-186 —0-110 0-066 
5 0-366 0-294 0-106 —0-106 —Q-212 —0-190 —0-090 0-090 
6 0-040 —0-018 —0-128 —0-218 —(0-236 —0-158 —0-028 0-100: 
7 —0-208 —0:230 —0-246 —0-:246 —0-194 —0-:078 0-032 0-110: 
8 —0-264 —0-258 —0-240 —0-176 —0-088 0-024 0-102 0-088 
9 —0-152 —0-148 —0:096 —0:026 0-046 0:106 0-136 0-046 
10 0-020 0-012 0-054 0-106 0-130 0-144 0-130 —0-008. 
11 0-134 0-118 0-152 0-160 0-136 0-100 0-058 —0-052 
12 0-130 0-128 0-128 0-116 0-070 0-038 —0-018 
13 0-076 0-050 0:052 0-020 —0-020 
14 —0-034 —0-044 —0-:052 


The centre-line amplitudes are a recomputation of those given by Steward and shown in figure 6. 


A series of longitudinal and transverse sections were then plotted from the 
above figures, and contour points were taken from the section graphs at amplitude 
values ,// corresponding to intensities /=0, 2, 4, 6, 8, 10, 20, 30, 40, 50, 60, 70, 
80%. 

The resulting contour map of the light distribution is seen in figure 8. * 
In this figure, the zero-intensity loci are shown dotted, and the outline of the 
geometrical image is included. 

Light contours for the case B=3 
It has recently been brought to my attention that the intensity contours 


for the case of B=3 have been very fully worked out by Nijboer (1946). This 
plot is reproduced in figure 9 as it may not be readily available. 


Light contours for the case B= 6:14 radians (352°) 

This case has been integrated by Martin (1922 a), and the central part of the 
pattern plotted. It was later verified photographically by the use of a specially 
designed microscope objective having the desired amount of coma (1922b). 

Light contours for the case B= 20 

For this case, 16 zones were taken of radii y/(1/16), \/(2/16), .. . 1, respec- 
tively, the equiphase-points of the successive zones being, therefore, 1-25 (radian) 
units apart. ‘The amplitude-contribution graphs were plotted for some 190 
points in the focal plane, and each graph was integrated as before. The resulting 
net amplitudes are tabulated in table 2. The contour pattern (figure 10) was 
plotted as before by drawing a series of longitudinal and transverse sections 
and marking the points for the amplitudes \/J corresponding to J=0-5, 1, 2, 4 
6, 8, 10%. 

* This pattern differs somewhat from that given by Steward (1926 a, p. 156) ; there are also 


small differences between the amplitude distribution along the centre-line of the pattern as given in 
table 1 and as given by Steward. 
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3° fob}? 4° 


Figure 13. Photographs of a comatic star image showing progressively increasing amounts of coma. 
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Figure 14. 


Photographs of a comatic star image showin 
increasing departures from the focal plane. 


g the effects of progressively 


Contours of intensity for B=4 


Y 
y=0 


Figure 8. 


yer 
Figure 9. Contours of intensity for B=3 (from Nijboer). 
Reproduced by kind permission of the Elsevier Publishing Co., Inc. 


Contours of intensity for B=20. 


Figure 10. 
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The diffraction structure of the elementary coma image I! 


As in the previous cases, the 
loci of zero intensity are shown 
dotted, and the outline of the 
geometric coma pattern is included 
for reference. The pattern was 
computed to Y =50, which is just 
beyond the centre of the circle in 
the geometrical image belonging 


to the marginal zone. mS = 225 
8 Zee 


Light contours for the case B=40_ gaara 


For this case, the number of 
zones was increased to 32, of 
radii 4/(1/32), 4/(2/32) etc. the 
equiphase-points of successive 
zones being 1-25 (radians) apart 
as for the B=20 case. This was a 
convenience in computation since 
the oblique distances QO from all 
the equiphase foci to the integra- 
tion points had already been 
evaluated. The labour of com- 
puting the net amplitude at each 
integration point was considerable, 
as it involved computation of 32 
four-digit products of ~ and O, Figure 11. Some typical amplitude-contribution 
determination of the J, curve for EE I I A 
each product and integration by 
planimeter. ‘Tl’o show how rapidly the amplitude contributions alternate from 
one zone to another, a few examples are included in figure 11. 

In all, the net amplitudes at some 256 points in the image plane were computed 
in this way (table 3). The longitudinal and cross section graphs were then 
plotted, and contours were drawn (figure 12) at amplitude levels corresponding 
Eominrensities Woy lt) 2.927553 3-354 9/5. 
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S38 PHOTOGRAPHIC CONFIRMATION OF THE 
COMPUTED PATTERNS 

The problem here is to obtain a lens in which there is a considerable amount 
of coma and negligible spherical aberration and astigmatism. Such a lens 
might be, for example, the rear half of a symmetrical system designed to give 
excellent anastigmatic definition at unit magnification, with parallel light entering 
through the central stop. 

Fortunately a suitable lens of this type was available. ‘The rear half had a 
focal length of about 8 inches and an aperture of f/5-6, the field out to 5° from the 
axis being flat and highly corrected for astigmatism. ‘The spherical aberration 
was so small that even at f/5-6 the axial image in monochromatic light was a 
perfect Airy disc. 

With this lens set at 4° obliquity, an attempt was made by adjusting the iris 
diaphragm to reproduce as closely as possible the computed pattern for the 
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The diffraction structure of the elementary coma image E57 


case B=40. ‘Then, since coma is proportional to image height, it was expected 
that each degree of obliquity would represent 10 units of B. However, this 
operation proved to be surprisingly difficult, and in the actual series of enlarged 
photographs of the star image at every half degree of obliquity, reproduced in 
figure 13, it will be seen that the computed patterns for B=4, 20, and 40 very 
Closely resemble the photographs at }°, 2°, and 34° respectively. The resem- 
blance between the computations and the photographs is remarkably close, 
even to the shape of the individual little spots and bands of light in the images. 
The original star-image photographs were made through an 8mm. apo- 
chromatic microscope objective on ordinary 16mm. motion-picture panchromatic 
film, in sodium D light, the negatives being then further enlarged in making 
the prints. 
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Figure 12.. Contours of intensity for B=40. 


§9. OUT-OF-FOCUS COMA IMAGES 

The problem of computing the light distribution in an out-of-focus coma 
image is very difficult, as the convenient fact that each zone forms a sharp equiphase 
focus at its sagittal image can no longer be utilized. Consequently, use was 
made of the lens system mentioned in the last paragraph to provide this informa- 
tion. The lens was tilted to 44° obliquity (corresponding approximately to 
B=50), and a succession of star-image photographs were made for a series of 
focal settings at successive intervals of 0-1 mm. out of focus. It was found to 
be immaterial whether the image plane was chosen within or beyond the best 
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focus plane, thus again confirming the very high degree of correction of all the 
aberrations other than coma. These out-of-focus star images are shown in 
figure 14. 
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Diffraction and Optical Image Formation 
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The Thomas Young Oration, delivered 24 September 1947 


very slow development of the fundamental concepts and methods. ‘The wave 

theory of light seemed not only well established, but even nearly completed, 
by about 1820, after the work of Young and Fresnel. However, the first applica- 
tion to the resolving power of optical instruments was made by Airy in 1835, the 
second by Helmholtz and by Abbe forty years later. The next step, the extension 
to the case of lens errors, was not made until about 1900 by Strehl, and 1920 by 
Conrady and by Richter. A somewhat different branch, that of coherence, 
begun hesitatingly by Verdet in 1860 and developed by Michelson in 1890, 
found its practical application about 1920 at Mt. Wilson, and its further theoretical 
foundation by van Cittert (1934). You will understand that I have often been 
. astonished to find that problems of such old standing still showed themselves 
open to further development. 

I shall treat of three different, though interconnected, subjects: the coherent 
background, the degree of coherence and the diffraction theory of aberrations. 
As an introduction, I would demonstrate a few diffraction experiments. By 
the aid of a small arc lamp, condenser and vertical slit I throw the shadow of a 
thin vertical needle on the screen. In order to make the details visible through 
the whole room, I enlarge the shadow, in the horizontal direction only, by a short 
focus cylindrical lens. ‘The external fringes are clearly seen, but also the internal 
fringes, which Thomas Young explained by interference between the two beams 
diffracted at the edges (figure 1(a)). With a wedge-shaped needle (figure 1(4)) 
the different behaviour of the two kinds with increasing thickness of the needle 
is seen ata glance. The internal fringes soon become too dark to be seen. Yet 
with another thicker needle they reappear very clearly, a dark fringe at the centre 
(figure 1(c)). I have used a simple trick here: the needle is double, through 


if one reviews the historical development of this subject, one is struck by the 
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(a) (b) (¢) 
Figure 1. 


Figure 2. Figure 3. Diffraction figure of a slit, above without background, 
middle with coherent background, below background alone. 


Figure 4. F. G. Pease at the eyepiece of his 15-metre stellar interferometer 


Figure 5. Astigmatism patterns for B=8, 15, 26 and 45 respectively. 
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Figure 6. Astigmatism for B=1, left without background, middle same with 
background, vight background alone. | 


Figure 7. Coma patterns with B=0-7, 2:5, 10, 20 and 50 respectively. 


Figure 8. Astigmatism with B=12, without background, with background in 
phase and with background one quarter behind. 


(a) (b) () 
Figure 9. Astigmatism, 8 =17, (a) complete pattern, (b) edge effect alone, (c) pattern at focal line. 
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the narrow slit between some light passes and by diffraction is spread over the 
dark shadow. Evidently, the background thus created is coherent with the 
faint light of the fringes. Therefore the latter are effective here with their 
amplitude. A simple calculation shows that the fringes are much more easily 
seen in this way. ; 

In order to obtain a similar effect with the light bent into the shadow at a 
single edge, I use a different artifice: the shadow-throwing screen is made slightly 
transparent. A number of fringes is now seen in the shadow (figure 2). In 
general it is found advantageous to combine both methods. Thus the Fraunhofer 
diffraction image of a slit was thrown on a coherent background in the following 
way: ‘The slit is covered by a plane-parallel glass covered with a transparent 
layer of a strongly absorbing metal (silver or aluminium). In this layer a narrow 
scratch is made exactly in the centre of the slit. The extra light through this 
scratch gives rise to a much broader, but not too much weaker, diffraction image, 
the central fringe of which may cover the whole diffraction image of the slit 
(figure 3). 

In these cases one can in different ways also change the phase of the back- 
ground at will. ‘The result is therefore that the amplitudes and phases at various 
points of a diffraction image can be observed experimentally by aid of the coherent 
background. 

In this method there is further a clear parallelism with theory. In order 
to define clearly the meaning of the phase difference between points of a diffraction 
image, the theorist must be sure to introduce a surface of reference, which may be 
plane, cylindrical, spherical, etc. in a more or less arbitrary way. Indeed, 
controversies have sometimes arisen through neglect of this. In the same way 
the experimenter must introduce an auxiliary coherent wave, which may radiate 
from a more or less arbitrary point or line. 

The question of partial coherence started with Verdet, who asked at what 
distance apart two points on a screen illuminated by the sun would still have 
coherent vibrations. He found that this was determined by the apparent diameter 
of the sun, the actual distance being less than 1/20 millimetre. Let us imagine 
an experimenter who wants to verify this. He will take a piece of tinfoil and 
prick very small holes in it, in pairs of various distances. ‘Through each hole 
a cone of light will pass and adjacent cones will overlap and show interference 
fringes. But if they do not, shall we call the adjacent cones zmcoherent? | agree, 
but I must warn you that this is quite a daring step, undertaken only recently. 
Indeed, opticians have been very cautious, calling vibrations incoherent only 
when they came from different sources, thus making sure that their haphazardly 
changing phases would be statistically independent. Well, as soon as we take 
the more daring point of view, a new opportunity presents itself, namely to call 
vibrations partially coherent when they give fringes of lower visibility and to 
define their degree of coherence y to be equal to the number between 0 and 1 which 


expresses their visibility : 


Perce I eee 
visibility = —“*“— = degree of coherence y. 
Imax + Lin 
This new concept of degree of coherence leads to various remarkable results. 
For our sunlit screen y vanishes at a distance of 0:07 mm., but at larger distances 
it rises again and goes up and down many times. More exactly, its course 1s 
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equal to that of the amplitude of the diffraction image of a star formed by a telescope 
objective of angular aperture equal to the sun’s diameter (Airy disc with rings). 
It would carry us too far from our subject to discuss the application of the same 
theorem (discovered by van Cittert 1934) to the role of a microscope condenser 
or of a condenser in front of a spectroscope slit. 

Another general theorem, on the contrary, is directly connected with image 
formation. It is the theorem about the propagation of y, which states that a 
knowledge of the distribution of intensity, of degree of coherence and of relative 
phases in any surface intersecting a beam of light enables us to calculate the same 
quantities at a following or at a preceding surface. For instance, the image 
that will be formed in a photographic camera—i.e. the distribution of intensity 
on the sensitive layer—is present in an invisible, mysterious way in the aperture 
of the lens, where the intensity is equal at all points, namely in the distribution of 
y in this aperture. And if you ask for the mathematical connection between 
the two: one is the Fourier transform of the other. ‘This was probably known, 
in an incomplete form, to Michelson in 1890, but for lack of the requisite term 
and even of the requisite concept, he could not adequately express it. In such 
cases one feels the truth of E. Mach’s statement that science serves to economize 
thinking. Indeed a single term may stand for a whole theory, may convey its 
‘deas, theorems and concepts. 

Michelson (1890) had suggested that it is possible to obtain the apparent 
diameters of stars from the visibilities of interference fringes, or, in our terminology, 
to find them from a determination of the degree of coherence as a function of 
distance apart. In the practical execution of this idea at the Mt. Wilson Obser- 
vatory, our minute holes in the tinfoil were represented by two eight-inch mirrors 
under 45° which were movable along a six-metre steel beam, mounted across 
the opening of the 100-inch telescope. ‘The mirrors reflect the light from a 
star towards the centre of the beam. ‘Two fixed mirrors mounted there throw 
the light into the telescope. The observer first sees two star discs which.he 
brings into coincidence by adjusting the mirrors. He then estimates the visi- 
bility of the fringes that appear and repeats this for various distances of the 
first mirrors. F. G. Pease, the designer of the 100-inch telescope, is to be 
credited for most of this work. It may be said that he devoted much of his time 
during the last fifteen years of his life to measuring degrees of coherence—without 
knowing it. All the same he measured various star-diameters and “‘resolved”’ 
Mizar, before only known as a spectroscopic double star (distance 0-011”). 
Figure 4 shows Dr. Pease at the eyepiece of a specially constructed 15-metre 
instrument. Perhaps future astronomers will build larger instruments of this 
kind, say of 50 or 100 metres, which will indirectly show details down to one 
thousandth of a second of arc. 

This indirect method of studying celestial objects may well be compared 
with the study of crystal structure by x-ray diffraction. There also the synthesis 
from the Fourier transform to the image cannot be obtained by the direct optical 
method, that is, the image cannot be seen but must be calculated. 

. I shall dwell somewhat longer on my last subject, the réle of diffraction in 
image formation in the presence of lens errors. The prevailing attitude among 
opticians was even recently the following. Diffraction causes a certain unavoidable 
deterioration of the ideal point image, which may be expressed by the radius a 
of the Airy disc. Any aberration present will also give a certain diffusion, 
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expressible by a ‘“‘radius”’ 6 of the corresponding geometrical pattern... Both 
together may then be assumed to give a diffusion a+b. There is one thing in 
favour of this crude estimate: the tolerances for aberrations deduced from it 
are very much on the safe side. 

My own work on the change of the Airy pattern caused by small aberrations 
was started in 1934 and completed by Nijboer in his thesis of 1942. A short 
survey of the results must suffice here. To begin with, the aberration will of 
course be expressed, not in terms of light rays and their intersection with the 
receiving plane, but as deviations of the wave surface in the exit pupil from the 
ideal spherical form. Let these deviations be represented by V(y, x) as a function 
of the plane rectangular coordinates y, x in the circular opening, or as V(r, ¢) 
in polar coordinates, such that r=1 at the edge. — Secondly, this characteristic 
function V is developed in a series of orthogonal polynomials, 


Vig, $) =2 Damien 1) cos mo, “SI psct (1) ; 
which were specially constructed for the purpose. In the receiving plane with 
polar coordinates p, ys the resulting diffraction image has then, to a first approxi- 
mation, the amplitude 


A(p, b) = 2p {J4(p) a eros fear F +1(p) COS np}. eS One (2) 


Nijboer also gives the general form of the terms up to the fourth powers of the 
coefficients f. ; 

At the centre of the pattern, the point of maximum amplitude for small 
errors, the result is 
Ay=1—- B;,,/4(n +1). ene) 
The intensity at the centre, A?, is a good measure of the quality of the image 
(Strehl’s “ definition-brightness”’). A diminution of 10°% may well be tolerated, 
A,=0-95. Suppose, for instance that there is only ordinary spherical aberration, 
B19 must then be less than 1, whereas considered geometrically, this value would 
give a circle of least confusion of radius 6, or 3-1 times the radius of the Airy 
‘disc. Another remarkable advance lies in the fact that there are no mixed 
terms in (3), so that a higher aberration cannot be improved by small amounts 
of alower one. In other words, the balancing of aberrations has been completely 
attained by our introduction of orthogonal polynomials. 

The method described was thus fully successful for small errors.. Dr. 
Nijboer found it increasingly difficult for larger aberrations. I can show a 
slide with the pattern for astigmatism with B,,.=4. The formula for this case 
fills a whole page. As it appeared hopeless to get any further theoretically, 
especially to get an insight into the gradual transition to the geometrical pattern 
for increasing errors, we turned to experiment. Let me give a few details of 
the way we obtained pure third order errors. 


Astigmatism was obtained with a symmetrical feos: lens of 1m. focus 
and magnification one, thus excluding coma, turned into an oblique position 
through measurable angles up to 15°. The circular diaphragm remained fixed, 
perpendicular to the beam. The astigmatism is proportional to the square 
-of the angle. 

Coma was obtained from an ordinary achromatic telescope objective by 
shifting its components laterally in opposite directions. ‘The coma is propor- 
tional to the shift and to the third power of the aperture. 
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Spherical aberration was obtained from a single meniscus lens in a reversed 
position, the amount being changed by changing the aperture. It is proportional 
to the fourth power of the aperture. 

In all cases green mercury light was used and the diffraction pattern enlarged 
5-10 times by an auxiliary lens. The largest amounts of aberration used could 
always be determined experimentally from geometrical optics (axial distance 
of focal lines etc.). Let me show a few examples of astigmatism and of coma 
(figures 5, 6 and 7). Mr. Nienhuis, who made these experiments and whose 
thesis is to appear in a few months, now proceeded to the experimental investiga- 
tion of the formation of these patterns. 

As an example I take the largest coma pattern. The wave retardation, 
expressed in radians, is in this case 


Vix, 6) =38r c0sh=3BVy" 27), 2 eee (4) 
there being no need here for the polynomial. Let us start from the geometric 
optical pattern. As is well known, each zone of the lens aperture gives rise 
to a circle with displaced centre and which is described twice. In fact, the optical 
path to the point (7, ¢) in the receiving plane becomes 


Vy, 2) —99n 326) ee (5) 


and the point of intersection of the ray is found by equating the derivatives of 
(5) to zero: 


n=0V ay, 6=0V/ds, 


which in our case reduces to 7 =38(27?+7r?cos2¢), €=3fr?sin2¢. Therefore 
two rays from diametrically opposite points of a zone intersect in the same point 
7, € and must show interference. Substituting (6) into (5), the path difference 
with the principal ray becomes generally V —y dV /¢y —z OV /0z=V —3V = —2P, 
or twice this amount between the two interfering rays. Nienhuis finds that 
the observed appearance of the interference fringes is quantitatively explained 
in this way. In the lower part of the coma pattern, near the tip, the lower parts. 
of large circles will overlap with the upper parts of much smaller ones. The 
two crossed systems of fringes give rise to the diamond pattern observed. 

One further detail can better be illustrated on the astigmatism pattern. In 
this case the geometrical pattern is to be deduced in the same way from 
V(y, 2) = Br’ cos24¢ = B(y? —2") giving y=2By, €= —2Bz, that is, all zones. 
give concentric and proportional circles which are described in opposite directions. 
‘Therefore there is no overlapping here and we should expect an evenly illuminated. 
circle instead of the observed four-pointed star. In the same way as above,. 
we further obtain for the path difference in the receiving plane, 


V -2V= —V= —B(y? =2?). 


This means that our pattern, on interfering with an auxiliary spherical wave,. 
should show the same hyperbolic fringes, with reversed sign, that would be 
found in the lens aperture in a Twyman interferometer. The method of the: 
coherent background realizes this; figure 8 shows the result. 

The four-pointed star may be explained by the diffraction at the diaphragm-. 
edge. As is well known, the Fresnel diffraction at an edge of any form may 
be ascribed to rays emerging from the edge and spreading from the undiffracted 
ray only in directions perpendicular to the edge. In our case we must therefore- 
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expect diffraction streaks issuing from each point on the circle of least confusion. 
As this circle is described in the opposite direction, the streaks will turn against 
the radius vector. The geometrical problem is easily solved, the streaks will 
envelop an asteroid, the equation to which will be yt+¢#=(4)#. This was 
tested experimentally by inserting a metal disc in the diaphragm opening, leaving 
free only a narrow annulus. ‘The asteroid pattern then appeared, unobstructed 
by the intense geometrical pattern (figure 9(b)). The reverse is also possible: 
by throwing the shadow of a small circular opening on the lens aperture, the 
latter may be illuminated with an intensity decreasing towards the edge and 
vanishing at the edge itself. The edge effect was indeed absent in this case, 
only a simple circular pattern remaining. Such experimental tricks are not 
even necessary. It is found geometrically that the asteroid edge-pattern must 
remain unchanged, at least in form, when the focus is changed. It therefore 
comes out much clearer when the receiving plane is placed at one of the focal 
lines (figure 9 (c)). 

After these experiments the question arose: must we rest content with this. 
solution of the problem? Evidently the problem is in essence a mathematical 
one; we have no reason to doubt the validity of the relatively simple diffraction 
integral, which was also the starting point for the development in case of small 
errors. ‘The real difficulty was to find an asymptotic expansion for large wave 
numbers. Only a few months ago my assistant N. G. van Kampen attacked 
this problem again and found the solution. ‘The ‘‘ method of stationary phase”’ 
used goes back to Stokes and Kelvin. Its mathematical elaboration was well 
known to us, as it is due to my colleague J. G. van der Corput (1936). However, 
van Kampen had to extend it to two variables. Let me give a brief summary 
of the mathematical formulation. ‘The problem is to develop an integral of 
the general form 


tkf (y,2) S 
| e(o ae” dy de 


into an asymptotic series for large values of the parameter k, the integral being 
extended over the domain D, the boundary of which consists of a finite number 
of analytical curves. 

It is found that the ever increasing rapidity of fluctuation of the exponential 
causes the series to depend on the behaviour at a limited number of decisive 
points, which are of three kinds: (a) internal points at which the argument of 
the integrand is stationary, i.e. 0f/¢y=df/dx=0, (6) boundary points at which 
the argument is stationary along the boundary, 1.e. df/ds=0, (c) corner points, 
or boundary points where two analytical curves join. 

In the neighbourhood of any internal decisive point (yo, 29), the exponential 
may be partially developed into a power series 


exp (tkfy) exp {0R(a1.y? + 2oyo VF +92") HI +... thay, y? + tke y*} 
and this is integrated term by term between limits and —oo + «©. The resulting 
integrations are easily performed, they give a series with principal term 


ikf, |%11 %12 | 


TB € 


with additional details about signs into which I shall not enter. ‘The following 
terms are of order k-1, k-2 etc. Ina similar way, each decisive boundary point 
gives a series beginning with k*, k-? etc. and each corner point terms of order 
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k, k-2 etc. The contributions of all decisive points must finally be added. 
Of course it is no easy matter to prove the mathematical validity of the whole 
procedure, but Professor van der Corput has just now succeeded in it. 

Now, applying this to our diffraction integral, we may ordinarily put the 
amplitude g equal to one and f= V(y, z) —yn —2€, and consider a circular boundary. 
The internal decisive points are then to be found from (6) and the principal 
term (7) is that of geometrical optics, with phases and interference taken into 
account. Even the intensity agrees exactly, the square of the last factor of (7) 
corresponding to the concentration of light rays by the curvature of the wave 
surface. 

The second term comes from the boundary points, it corresponds to the 
beams diffracted by the edge. But here the theory gives more than the experi- 
mental treatment, which was not able to predict the intensity. Of course the 
new development also gives more than these two terms and these further terms 
could not be found in another way. We are now endeavouring to fill the gap 
between small and large errors by calculating some intermediate case, say astig- 
matism with 8 = 10, from both sides. 

It is especially the general insight, however, which gains very much by the 
discovery of the asymptotic development. It shows that physical intuition 
combined with experimental ability may go far towards elucidating the main 
characteristics of phenomena, but that only an adequate mathematical treatment 
-can give a satisfactory final solution. 
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ABSTRACT. Thestructure of PbS deposits condensed from the vapour in vacuo on to 
{001}, {110}, {111} and {443} rocksalt faces has been investigated by electron diffraction. 
The results suggest that the deposit atoms take up positions of least potential energy 
relative to the substrate, as far as is permitted by the disturbing effects of collisions of 
incident atoms with the initial deposit crystal nuclei, and by the limited surface mobility 
of the deposited atoms over the substrate. "This view is also supported by the nature of 


the changes in crystal orientation which occur when initially random deposits are heated 
in vacuo. 
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§1. INTRODUCTION 
BAD sulphide is a semiconductor which has important applications owing 
[o: its photoconductive properties and especially its electron-emission 
sensitivity to infra-red radiation. Outlines of the properties of sensitive 
deposits, and of their production by sublimation in vacuo or low pressure oxygen 
have been published by Lee and Parker (1946), Starkiewicz e¢ al. (1946), and 
Sosnowski et al. (1947 a, b). 

A recent investigation (Wilman 1948) of such deposits by electron diffraction 
established that lanarkite, PbO.PbSO,, was formed as an oxidation product 
in the sensitization of the PbS by sublimation and baking in oxygen. With 
a view to studying further the form of the oxidation product and its connection 
with photoconductivity, single-crystal lead-sulphide deposits have now been 
prepared by condensation on rocksalt substrates. Rocksalt was chosen as 
likely to initiate a parallel orientation in the PbS because its axial length (5-639 a.) 
is within 6% of that of PbS (5-929 a.). The structure of these PbS deposits 
as shown by electron diffraction is described below. The results afford a picture 
of the mode of growth of a semiconductor condensed from its vapour on to an 
ionic substrate and provide a useful comparison with the growth of similarly 
prepared deposits of metals (Kirchner 1932, Lassen 1934, Lassen and Briick 
1935, Rudiger 1937, Wilman 1939). 


§2. EXPERIMENTAL 

Fresh rocksalt {001} faces about 8x5mm. were prepared immediately 
before use by cleavage from an almost perfect natural crystal. The specimen 
was handled only with tweezers, and any stray fragments were lightly brushed 
off the face with a camel hair brush. Other faces used were ground on large 
cleavage blocks of rocksalt and smoothed on 0000 emery paper before lightly 
etching during a second or two in a fast stream of tap-water, followed by plunging 
immediately into absolute alcohol to stop further etching and to remove water 
and solution. Diffraction patterns of Kikuchi lines and spots elongated per- 
pendicular to the shadow edge from such surfaces showed that the resulting 
smooth surfaces were free from random recrystallized NaCl. 

The lead sulphide was prepared by precipitation from concentrated lead- 
nitrate solution by H,S gas, followed by repeated washing with distilled water 
and drying in air at about 80°c. After driving off excess sulphur by heating 
in vacuo, a few milligrams were transferred to another pyrex tube and sublimed 
in a vacuum of the order of 10-?mm. Hg on to rocksalt substrates. 

The deposits were normally allowed to cool for at least five minutes before 
letting in air and transferring to a Finch-type electron-diffraction camera (Finch 
and Wilman 1937). A camera length of about 47cm. and 50-65kv. electrons 
were used. For transmission examination the rocksalt substrate was dissolved 
away in distilled water and the floating PbS film was picked up on a nickel gauze 
and dried zm vacuo in the diffraction camera. 

The rate of deposition and thickness of deposits were controlled roughly 
by eye to standard conditions, namely a rate of the order of 40 a./sec. and thick- 
nesses about 50, 250 or 1000. (corresponding respectively to light yellow, 
translucent brown and opaque, or nearly so, by transmitted light). The thick- 
nesses were estimated for a few cases by subliming known weights of PbS on 
to a measured area and visual comparisons with these were made ap: the remaining 
specimens. 
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$3. THE STRUCTURE. ORS GHE DEPOSITS 

The results are outlined in tables 1 and 2. The deposits prepared were 
mostly of medium thickness in view of the purpose for which they were later 
intended. 

Deposits on apparently good ground and etched surfaces sometimes consisted 
of randomly disposed PbS crystals, even at the higher temperature range where 
strong orientation of the PbS was usually observed, all other conditions being 
kept constant. These results are not recorded in table 2 because they were 
attributed to the accidental presence of submicroscopic randomly disposed 
NaCl crystals on the main rocksalt substrate. 


(i) Orientations of the PbS 


In addition to the orientations listed in the tables, small proportions of 
randomly disposed PbS crystals were occasionally indicated by faint PbS rings 
in the diffraction pattern as in figure 2, but the orientation was usually very 
strong as shown by well-defined diffraction spot patterns, figures 1 and 6. . In 
the case of deposits on {110} and {111} NaCl faces below about 150° c. the only 
orientation observed was with {111} planes parallel to the substrate with other- 
wise random azimuthal distribution. 

(ii) Twinning 

About one-third of the deposits on the {001} NaCl faces were found to possess 
strongly twinned structures on either {111} or {332} planes as shown by figure 2 
and figures 3 and 4 respectively. The {332} twinning appears not to have been 
observed hitherto in the mineral galena. Deposits condensed from the vapour 
must of course grow from many nuclei formed initially on the substrate, not 
necessarily all in the same orientation; thus it is not clear whether these PbS 
crystals first grew in an orientation parallel to that of the NaCl and then twinned, 
or whether they were separate crystals in twin orientation. In the above cases, 
the crystals in the twin orientations to those which are parallel to the rocksalt 
are believed to arise actually by twinning from initially formed PbS crystals 
whose axes are parallel to those of the NaCl, because these latter crystals have 
in these cases developed faces (see (iii) below) which are parallel to the twinning 
plane. ‘The planes in the {111} and {332} twin crystals, which would be junction 
planes with the NaCl {001} substrate, have lattice-point distributions which are 
only similar in spacing to the NaCl along one direction, i.e. the <110> rows 
of the NaCl, and growth in these orientations on the NaCl is, therefore, less 
likely than in the {001} orientation, though possible. 

Diffraction spots due to {hhl} types of twinning are most prominently observed 
in the <110> azimuth, their positions being then obtained by the rotation 
of the main pattern (the stronger spots) about the central undeflected spot, 
through 2 tan~* //hy/2, namely 70°31’ for {111} twinning (figure 2), 50° 28’ for 
{332} twinning (figure 3). None of the seven deposits on {111} NaCl nor that 
on a {443} NaCl face showed any evidence of twinning, and only one on {110} 
NaCl showed a few extra spots not yet explained. 


(1) Development of external crystal faces 


In most of the PbS deposits the crystals had tended to develop definite external 
faces to a greater or lesser extent as shown by elongation of the spots, especially 
those near to the shadow edge. In all cases the faces parallel to the substrate 
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Table 2. The structure of PbS deposits on {110}, {111} and {443} NaCl faces 
| 
No. of specimens with following characteristics 

Approx. | Approx. Order of 

temp. of | deposit crystal diameters (A.) 

substrate | thickness | Total Orientations 

(xcs) (A.) number observed |L-+ TERS}: 

20-50 | 50-200 | 20-50 | 50-200 


Deposits on {110} NaCl face 


Slight {111} 


ay ot ! | L, to substrate : | ee | J - 
100-200 50 1 Random = 1 = 1 
200-300 50 1 ier Ah 1 - 1 = 
250 2 Ne 2 1 il i 1 

Deposits on {111} NaCl face 
100-200 | 50 1 Seley 1 {7 oe 

| L substrate 

1000 1 Random = 1 = 1 
200-300 50 4 ||L* 4 4 - 4 ~ 
250 2 dee 2 - 2 ~ 

Deposit on {443} NaCl face 
200-300 250 1 ees 1 | 1 | = | di | - 


* \|L, denotes PbS cubic axes parallel to those of the NaCl. 
+ ||L. and LR are reckoned with respect to the substrate surface. 


were the most prominent and were often very extensive and atomically smooth 
as shown by the degeneration of the spot pattern into almost continuous parallel 
lines perpendicular to the shadow edge. ‘This is shown in figure 1, characteristic 
of those deposits which were kept at the deposition temperature (200-300° c.) 
for a further 30sec. after deposition. The larger size of the crystals in these 
deposits was shown by the appearance of diffuse Kikuchi bands, as well as by 
the sharpness of the vertical streaks. 

The patterns from one of the {332}-twinned deposits on {001} NaCl in the 
<100> and the <110> azimuths showed that the crystals in the main orientation 
parallel to the NaCl had developed {100}, {110}, {120}, and {335} faces, the streaks 
in the pattern corresponding to {335} faces making an angle of 40° with the plane 
of incidence in figure 3. ‘The elongation of the twin spots perpendicular to the 
shadow edge showed that the twin lattices also had developed boundary faces 
parallel to the substrate and therefore of {335} type with respect to the cubic 
axes of the twin lattices. In another twinned PbS deposit having a smaller 
proportion of {332} twinning there were strongly developed {001} faces, and 
{120} faces were prominent though {110}! and {335} were not observed. 

In PbS deposits on {110} and {111} NaCl faces there was slight elongation 
of spots near and normal to the shadow edge, showing a tendency of the crystals’ 
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Figure |. PbS on {001} NaCl; 
{00I} orientation , <110)Az. 


Figure 3. PbS on {00I} NaCl; 
{001}, <100>Az., {332} twinning 


Figure 5. PbS ool NaCl. 
{001} ead {Il} eee es 


Figure7, Random PbS onf{III} NaCl. 
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aie 2. PbS on {001} 
{001} ,<1T0> Az., {III} twinning. 


NaCl ; 


Figure 4. PbS on {001} NaCl; 
{001} ,<110) Az., {332} twinning. 


Figure6. PbS on {I10} NaCl, 
parallel orientation. 


Figure 8. Deposit of Figure7 
heated in vacuum to ~ 250°C, 
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in the deposit to form faces parallel to the substrate. In two specimens on 
{110} NaCl faces there was pronounced development of {001} boundary faces, 
inclined at 45° to the substrate, and the large extent of the spot pattern is also 
consistent with the development of such small projections on the deposit surface. 

The thin PbS deposit prepared on a {443} NaCl face had all its crystals very 
strongly oriented with their cubic axes parallel to those of the NaCl. The spot 
patterns were as free from arcing as those of figures 3 and 4, though without 
any elongations representing definite boundary faces. The pattern previously 
obtained from the {443} NaCl face showed Kikuchi bands and relatively few but 
strong spots elongated normal to the shadow edge, showing almost perfect 
{443} boundary faces. 

(iv) Lattice constants 

In two cases of PbS deposits at 200-300° c. on {001} NaCl faces, the lattice 
constants of the PbS were measured from the patterns obtained by transmission 
after dissolving away the NaCl. In one case the reference material was graphite 
(dy) = 1-230.) giving appg=5-931a., and in the other case NaCl crystallized 
from solution (dy9)=2:820 4.) giving appg=5-927 a. These values are close 
to the average of a=5-929 a. previously found for PbS layers formed by passing 
H,S gas over lead-nitrate or lead-acetate solution or sublimation of such material 
in vacuo in pyrex cells (Wilman 1948). 

{4b Eee ERE Ch OF HEATING PbS) DEPOSITS 
ON NaCl SUBSTRATES IN VACUO 

It has been mentioned above that the continued heating of PbS deposits 
in vacuo at 200-300° c. after deposition on the {001} NaCl faces, only resulted 
in elimination of twinned structure, increase of crystal diameter and strong 
development of faces parallel to the substrate. No change of the main orientation 
was observed. It seemed possible, however, that (1) deposits consisting of 
random crystals as a result of low mobility of the deposit atoms on the substrate 
at the temperature of deposition might develop a two-degree orientation relative 
to the substrate lattice when heated im vacuo; and (ii) that the rate and kind of 
orientation developed should depend on the temperature at which the heat 
treatment is carried out and on the distribution of atoms or ions in the substrate 
surface. 

To test this a 10004. thick PbS layer was deposited at 100-150°c. on a 
{111} NaCl face and was found to consist of random crystals about 120. in 
diameter. When this was heated for six minutes at 200-300°c. im vacuo it 
developed strong single-crystal structure (figure 8) with a {111} plane parallel 
to the {111} substrate surface and the cubic axes parallel to those of the NaCl. 
A further confirmation was obtained in the case of a 50a. thick random layer 
deposited at room temperature on a {001} NaCl cleavage face. ‘Though initially 
composed of random crystals of only 10-20. mean diameter, on heating to 
200-300° c. in vacuo for four minutes it developed strong two-degree {001} 
and {111} orientation similar to that in layers deposited at this temperature on. 
{001} NaCl. In both cases there was a large and rapid increase in crystal size 
at this temperature of heating, to a mean diameter of at least 300 a. though less 
than 500 a. 

§5. ORIENTATION OF NaCl ON PbS 

Several cases were observed, e.g. figure 6, where traces of NaCl had evidently 

been present on the walls of the pyrex tube and has sublimed on to the already 
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deposited and oriented PbS. ‘These NaCl deposits were found in all cases 
to have the same orientation as the underlying PbS deposit, whether this was 
oriented or random. 


§6. DISCUSSION 
(i) The orientations observed 

The orientations of the PbS crystals relative to the {001}, {110}, {111} and 
‘443? NaCl substrates agree with the conditions generally observed, that in such 
cases of strongly oriented overgrowths there is a similarity of periodicity in at 
least one set of parallel lattice rows in the deposit and substrate crystal. The 
development of {111} orientations on {001} NaCl at temperatures above about 
250° c. as well as the first-developed {001} orientations is analogous to the case 
of silver deposits formed on {001} NaCl faces over a similar range of temperatures 
(Briick 1936). On {111} and {443} NaCl faces, the consistent growth of PbS 
with cubic axes parallel to those of the NaCl, and not in more than one orientation 
geometrically equivalent relative to the surface plane, is analogous to the continued 
growth of NaCl on NaCl, building up a continuous pattern of positive and 
negative ions. Thus it would seem that not only the surface atoms govern the 
orientation in the growth of such ionic crystals but also those in the underlying 
net-plane. 

(11) Twinning 

The proportion of specimens having a {111} or {332} twinned structure is 
only about 30° compared with almost 100°{ for those of silver. This may 
perhaps be associated with the smaller difference of lattice constants of the 
substrate and PbS deposit (cf. Wilman 1939). On the other hand, it seems 
likely that here the PbS crystals grow initially in an orientation parallel to the 
NaCl and twin during their further growth, whereas silver is strongly {111} 
twinned in very thin layers; but the crystals in the four twin orientations grow 
together (Menzer 1938), so that surface regions of thick silver layers are formed 
of crystals whose axes are parallel to those of the NaCl. 


(11) Lattice constants 
The present estimations from transmission patterns show more accurately 
than was previously possible from reflection patterns (Wilman 1948) the striking 


constancy of the lattice axial length, 5-929 a., of PbS prepared both by chemical 
deposition and by sublimation in vacuo. 


(iv) Interpretation of deposit structure in terms of a hypothesis of lowest potential 
energy state modified bv surface mobility and rate of deposition of atoms 
from the vapour stream 

It is now clear that oriented deposits must arise through the initially deposited 
atoms taking up stable positions where they have lowest potential energy on 
substrate crystal surfaces which consist of a periodic distribution of atoms or 

ions. Similar views have been expressed for example by Stranski (1928). 

In like manner further atoms take up packing positions having least potential 

energy, and thus consistent with the normal deposit crystal structure and atomic 

distribution, i.e. orientation, in the first layer. The attainment of such a state 
is, however, dependent especially upon the rate of deposition and on the thermal 
movements and the strength of forces between the deposit and substrate ions 
compared with those of the deposit and substrate alone. Increase in temperature 
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of the substrate, by increasing the amplitude of thermal vibrations of the substrate 
surface atoms, should assist the growth of strongly oriented deposits, and this 
has indeed been confirmed by many experiments including the present results 
on PbS. 

The effect of rate of deposition on the number and size of the deposit crystals 
at any stage of the deposition on single crystal or amorphous substrates has not 
been quantitatively explored in the present experiments on PbS, but the general 
observations so far made may be explained as follows: ‘The atoms first falling on 
the substrate have a certain degree of mobility over the surface before losing most 
of the kinetic energy which they possessed on arrival; thus, sooner or later, they 
meet other atoms and aggregate to form small nuclei. The greater the number 
of atoms incident on unit surface area per unit time, the larger will be the number 
of these nuclei formed per cm?/sec. and the smaller their lateral extent, and also 
the smaller the distances between these nuclei on the still uncovered substrate. 
The higher the rate of arrival of the atoms per cm2/sec., the greater is the number 
of atoms bombarding each nucleus per sec., and these will therefore hinder the 
development of the orientation which would be imposed by the substrate and 
result in a correspondingly weaker orientation or even a completely random 
crystal deposition. This picture is analogous to the Brownian movement of 
small particles in a liquid due to bombardment by the atoms of the liquid surround- 
ing them. 

The above experiments on PbS show that, even with a relatively high rate 
of deposition of about 40 a. thickness per second, strongly oriented PbS deposits 
on NaCl were obtained at substrate temperatures of about 150°c. As in the 
case of silver and other metals under similar conditions, random or weakly 
oriented PbS deposits tend to be formed at or near room temperature; and no 
strong development of orientation with increase of thickness was observed, such 
as occurs when the crystals develop large plane faces perpendicular to the incident 
beam direction. (Kirchner 1932, Burgers and Dippel 1934, Beeching 1936, 
Burgers and van Amstel 1936). ‘These results, together with the observation 
that random-crystalline PbS deposits do not develop orientation unless heated 
well above room temperature, suggest that the PbS molecules have low mobility 
over NaCl surfaces at room temperature, but that the mobility rapidly increases 
with temperature. ; 

The PbS deposits prepared at room temperature tended to be orientated 
apparently because the substrate surface temperature became more or less 
raised by absorption of radiant energy from the hot test-tube wall from which 
the PbS was being sublimed. The temperature of the substrate at its surface 
is difficult either to determine or control. It is a result of equilibrium between 
energy falling on to the surface by radiation from the oven—or filament—source 
of atoms, and energy leaving the surface by re-radiation and by conduction to 
the lower layers of the substrate. It was found that PbS deposits about 50 a. 
thick on an NaCl-crystal substrate 2 to 3mm. thick were often weakly oriented, 
but when 10mm. thick NaCl crystals of about the same area were used, the 
50. thick PbS deposits obtained were of quite randomly disposed crystals 
only 10-20. in diameter. This demonstrates the need for using substrates 
of high thermal capacity, if necessary backed by a metal mass of high conductivity, 
to ensure that the substrate surface temperature is not much higher than that 
of the rest of the substrate. 
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(v) Recrystallization of initially random PbS deposits 


According to the above view of the dependence of crystal growth on the 
deposit conditions, we would expect that even if the crystals of the deposit are 
initially randomly disposed, subsequent heat treatment should enable rearrange- 
ment of the atoms to take place in such a way that the state of least potential 
energy with the substrate, consistent with the surface mobility, is attained. 
This state should be the same as that which would have been obtained if the 
layer had been deposited at this temperature. ‘The results of the experiments 
in §4 do, in fact, lend further support to the hypothesis that the lowest potential 
energy state tends to be attained as far as the limited surface mobility of the 
deposit atoms permits. In this connection it is noteworthy that both {001} 
and {111} orientations of the PbS were developed together on the {001} NaCl 
face when heated zn vacuo to 200-300" c. 

he fact that random crystalline (or amorphous) deposits on heating take 
up an orientation characteristic of the substrate and its surface temperature, 
which under these conditions becomes that of the heating enclosure, enables 
the relation to be studied between the orientation of thin deposits and the nature 
and surface temperature of the substrate, without the uncertainty of measurement 
of the actual surface temperature during the deposition process. This should 
be especially helpful in the case of deposit materials which show a variety of 
orientations within a relatively narrow temperature range. 


§7. SUMMARY 

The structure of PbS deposits condensed from the vapour im vacuo on to 
{001}, {110}, {111} and {443} rocksalt faces has been investigated by electron 
diffraction. ‘The results obtained agree with the theory that the deposit atoms 
take up positions of least potential energy relative to the substrate, as far as is 
permitted by the disturbing effects of collisions of incident atoms with the initial 
deposit crystal nuclei and by the limited surface mobility of the deposited atoms 
over the substrate. Further confirmation of our theory is afforded by changes 
of orientation of initially random deposits when heated zm vacuo, and lead to 
the conclusion that such experiments should enable a more precise investigation 
to be made of the relation between the deposit orientation and the temperature 
of the substrate surface. 
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SUN O DiCiGas TON 


N a recent article (Lattes 1947), we gave an account of a determination, by 
] the method of grain-counting, of the ratio of the masses of the z- and ju-mesons 
observed in photographic emulsions. We have found more than thirty events 

in which the secondary p-meson produced by a z-meson at rest comes to the end 
of its range in our emulsions. ‘This new material, of which the main features are 
summarized in table 1, allows us to increase the accuracy of the determination, 


Table 1. Ranges of z- and u-mesons in individual examples of the p-decay 


Ranges of Ranges of 


ee m-meson p-meson Ey sie m-meson p-meson Ey 

o (microns) (microns) CINE 4 (microns) (microns) (ive) 
1 134 613 4-08 i 477 585 397, 
2 1287 615 4-09 18 12 685 4-36 
g 124 591 4-00 19 52 644 4-21 
4 106 633 4-16 20 62 594 4-00 
5 453 615 4-09 2] 230 570 3-91 
6 105 567 3°91 22, IANS 570 8.91 
7 906 606 4-05 23 488 680 4-34 
8 248 627 4-13 24 96 667 4-29 
9 231 648 4-22 25 409 626 4-13 
10 81 590 eo) 26 91 623 4-12 
11 57 609 4-06 pH | 990 602 4-04 
12 634 567 3-90 28 92 619 4-10 
is 48 611 4-06 29 15 598 4-01 
14 1297 603 4-04 30 91 639 4-18 
15 i! 591 4-00 Si 2290 660 4-27 
16 57 590 3:99 


Mean values: R,=614 4; E,=4-08 Mev. 


and in this article we describe in more detail the method of making the measure- 
“ments and the results which follow from them. We find that the mass ratio 
m,,m, =1:65 +0-11, the limits corresponding to the ‘probable error”’ arising 


* An account of the main features of a number of lectures given in Copenhagen in December 
947, in Zurich in January 1948, and in a number of colloquia in English Universities. 
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from statistical fluctuations alone. It is reasonable to assume that the u-mesons 
appear as a result of a spontaneous decay of the heavier 7-mesons, and that the 
former are identical in type with the ordinary cosmic-ray mesons of mass 200 m,. 
It then follows from the present results that the momentum balance in the decay 
is provided by a neutral particle, a neutral meson with a rest-mass equal to 
115 +30 m,, where m, is the mass of the electron. The observations also allow 
us to draw conclusions about the decay constants of the different particles, and 
we have thus determined the following limits for the half-value periods : 
T_ > 0-4 x 10> "sec.,. 7,52 10- sec, andy 2 67a, acce? where 7, is the 
half-value period for the postulated neutral meson. 


$2°sPRINGCIRLE OF LAE MEDPEO., 

In determining the quantity m,,/m,, by grain-counts, we make the usual assump- 
tion that a particle with a charge |e|, energy EMev., and range in the emulsion R 
microns, loses energy at a rate which depends only on its velocity, v, and is 
independent of its mass, m: OE/OR=f,(v). It is then easily shown that 
OF /OR=f,(R/m) and E=md¢d(R/m). 

We assume, secondly, that if N(R) is the total number of grains in the track of a 
particle of range R microns, the number of grains per micron, 0N/0R, in an element 
of length of the trajectory, is a function only of the rate of loss of energy, 0E/0R, 
of the particle. It follows that 

N-Sml (Rin). See ee eee (1) 
Equation (1) allows us, in principle, to determine the ratio m,/m,, from the observed 
numbers of grains in the track of 7- and u-mesons. 

In making grain-counts, we employ a x 95 ‘‘ Achromatic”’ objective with x 10 
eyepieces, and count the numbers of grains in successive intervals along the 
trajectory of length 50. ‘These intervals are measured along the projection of 
the track on the plane of the emulsion. With this procedure we shall tend to 
over-estimate the grain-density in a track which dips relative to the plane of the 
emulsion, through ‘‘fore-shortening”’. We have therefore confined our measure- 
ments to those events in which the angle of dip, 8, at every point in the trajectories 
of the primary and secondary particles is never greater than 15°. The errors due 
to ‘‘fore-shortening”’ are then small. 

A second source of error exists owing to the variation of the degree of develop- 
ment of an emulsion with depth which occurs when a plate is processed by con- 
ventional methods. ‘This has the consequence that the grain density may be 
different in the tracks of two particles produced contemporaneously, in regions of 
equal ionization, if they are situated at different depths in the emulsion. 

Thirdly, as we pointed out in a previous communication, there is not a 
completely random distribution of grains in the gelatine in the emulsions at present 
available. ‘‘Islands”’ appear to exist in which the concentration of silver-halide 
grains is significantly higher, or significantly lower, than we should expect to find 
as a result of random fluctuations alone. This feature seems to be particularly 
noticeable in boron-loaded plates. ‘The measurements show, however, that the 
influence of this factor on the total numbers of grains in the relatively long tracks on 
which our observations are made is not more serious than the purely statistical 
fluctuations associated with the finite number of grains in a track. 

All the above sources of error will effect the measurements on 7- and 4-mesons 
indiscriminately; and the errors they produce in the determination of the mean 
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Nn 


value of the mass-ratio will be small if the measurements are made on a sufficiently 
large number of events. They contribute, however, to the variation among the 
individual values of the mass-ratio. ; 

The following difficulties cannot at present be avoided. First, at the end of 
the range of a particle the ionization is sufficiently intense to produce an almost 
continuous succession of grains which cannot be individually resolved. We 
refer to this effect as ‘‘clogging’’. No accurate counts can be made in this region 
of a track, and we therefore adopt the following convention. We define the 
number of grains in a segment of the trajectory of length 7, in which there is a 
contiguous succession of grains, as equal to r/a, where a is the mean “ diameter”’ 
of the grains. In meson tracks “clogging”’ is serious only in the last 501 from 
the end of the range, and for such short tracks no reliable mass determinations are 
possible. 

Secondly, when the rate of loss of energy of the particle is very small, the grain 
density becomes so low that the counts become unreliable because of the general 
“background” of grains. ‘This sometimes makes it difficult to decide whether a 
particular grain was produced by the passage of the particle or is due to “fog”’. 
In our emulsion this difficulty begins to become appreciable at 400 and 6504, 
respectively, from the ends of the ranges of x- and z-mesons. 


§3. CONSISTENCY OF GRAIN-COUNTS 


Even in those parts of the trajectory of a particle in which ‘‘clogging”’ is not 
serious, the individual grains cannot always be completely resolved. Itis therefore 
necessary to adopt a convention that a particular size of grain, for example, shall be 
regarded as two, that very small grains shall be neglected, etc. We have found 
in practice that if different observers establish a common convention, by reference 
to suitable photo-micrographs of characteristic tracks, it is possible to obtain the 
following degree of consistency between different measurements. An experienced 
observer can determine the total number of grains in the last 400 » of the track of a 
particular meson with a consistency of 3°, and there is a similar degree of con- 
sistency between the observations of two different observers. In practice we have 
made a succession of five determinations of the number of grains, N, and N., in 

_ each of the tracks of the 7- and p-mesons of a particular event, making counts on 
the two tracks alternately. The ratio of the numbers, N,/N,, on which the mass 
ratio depends, is then found to be consistent to about 3%. Because of the finite 
number of grains in a track, the statistical fluctuation in the numbers N,, or N,, to 
be expected in going from one event to another is of the order of 6%, and the 
difference in the average number of grains in equal lengths of the tracks of 7- and 
-mesons about 13%. 


§4> RESULTS 


In table 2 we show the results of grain counts on those events which satisfy the 
conditions that at no point in either of the trajectories is the angle of dip greater 
than 15°. The table shows the number of grains, N(R), in a residual range R 
microns—for intervals 50 long—up to a maximum of 400 for 4-mesons and 
of 650, for 7-mesons. ‘The numbers in each vertical column correspond to the 
measurements on a particular event, so that the results for the two types of 
mesons can be compared. 

It will be seen that the corresponding results for different events are subject 
to considerably greater fluctuations than those corresponding to statistical effects 


Table 2. Grain counts on the track of 7- and pz-mesons 
Event No. 2 14 7 WZ 5 17 8 9 1 4 3 
Range No. of grains in tracks of 
(microns) e-mesons 

50 66 55 54 Te 61 60 60 66 45 33 69 
100 109 97 90 127 1055) 104 MO0Se eto 80 91 118 
150 144° 1:25) 1 1G GR Sent OO 97 125162 
200 179 157. 132 2045 180" 462° 18355 205s 5 1105 sot ee 
250 211 184 157 235° 216 “T98" “210m 9242) 142 Sion a 
300 236 = 201 184-266 —239 233. 2438°99276° MG7" 196 ezG7, 
350 265s Ope? Oil 287-2005 Ol DTA BOS See? 24ers 
400 289: 233: . 220 » 9312) «293. 285 “301-524 1ha20 456 2S eee 

m-mesons 

50 75 57 66 72 68 65 71 80 41 58 WS 
100 129 96> 114 0p 121 103 3 12ore 142 81 101 1 
150 171 hs Geto NEO a7 lee AS 194 «111 140 
200 213 172) 180) 19322 Oe ee li) 
250 26005 9) 2008 20S 23455 Oe eS 
300 20 Seon 28 ee. OMNES OIE OS 
350 330. 257 243 3000) 30725286 
400 SOs SP ESD SD) Bil 
450 391 S005 292) 36:7 eso, 


100 


100 200 300 500 700 1000 100 200 300 500 700 1000 
Range (microns) Range (microns) 

Figure 1. Logarithmic plots of the number of grains in the tracks of 7- and y-mesons for 

four examples of the y-decay. (d) shows the results for the event which displays the 


least difference in the number of grains in the tracks of the associated 7- and «-mesons 
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alone. We attribute these variations to “ fading” and to the other factors discussed 
above. The results for four typical events are represented in figure 1, in which 
log N is plotted against logR. Figure 1(d) corresponds to the event which 
displays the least difference between the results for the 7- and jz-mesons. 

In order to compare the results for the different events, we normalize the values 
for the different observations. We first determine a normalizing factor for each 
event, k,, defined as the ratio of the average number of grains in the last 400 py 
of the different -mesons, Nj” (400), tothe number of grains, N’ (400), in the same 
length of the particular ~-meson under consideration : k,;= N” (400)/N; (400). 
Table 3 shows the values ./(R) obtained by multiplying the numbers N(R) for 


Table 3. Values of (Rk), for 7- and p-mesons, after normalizing to 
correspond to a constant value of 4”, (400) 


No. 2D 14 a, 12 5 7 8 gy 1 as 3 N™R) 


k; 0-956 1-197 1-268 0-894 0-952 0-979 0-927 0-818 1-368 1-112 0-823 


Range 
(microns) fe-mesons 
50 64 66 69 64 58 59 56 54 62 59 57 | 60+2 


100 ings N14 I NOX) TO 97 Ys Oil 9% | 104=-2 
150 NED WSR IO) ISO ISIS IS IBID GNSS I ISD? 
200 7 Smee Onn Oe S27) Ol OR OS ae OS mmm / nl OO | male/slee es 
250 Vy ZS SY DAO) Oey IN PRION Ry TSR) AND) |) OWES 
300 ys ia) PTB S) UB} NS I} PY PNG ty PIR KG) || PEK) 
350 BIO AY DD ADI PIB SOD Ss ON DT PS) KGS}. || DOSES: 
400 ZI) 2D BR BY BID BIS VP DD BE DIE) BD || Vi 228 


7-mesons 

50 AZ, 68 84 64 65 64 66 65 56 64 62 66+2 
100 We ATS) IAS MO SS OE altho. || Wilko se 
150 iOS ies al = aks 7 lil a ISD ales 160+3 
200 PQ) Aer PRS ily Ailik Ney AI ARS 198+3 
250 eee oo e204 LOD 23S ae 217i 239 232-4 
300 23 Sn nS ee SO OE TO 263-4 
350 SIS 308 BOs Aes BWP AeKo) 290+5 
400 B50) 9338-5 340-297" 324 9306 320-5 
450 SI SH SHO SMS sh) 348 +6 
500 Ail SB Saxe) Bio 0 
550 435 425 421 386 404+-7 
600 461 450 435. 417 426+7 
650 490 469 455 449 8 


The values of 4” (R) for z-mesons have been obtained by determining the average number of 
grains, AV "(R), in successive intervals of length 50 », and summing over the total range R. 


the different events given in table 2, by the appropriate normalizing factor k;, 
i.e. W(R)=k,N(R). It may be emphasized that the mass ratio m,/m, depends 
only on N,/N,, and the values deduced from the observations on the individual 
events will not be changed by the use of the normalizing factor. Table 3 also 
contains the mean values for the number of grains, 4” ”"(R), in a range R microns. 
The results are represented diagrammatically in figure 2, in which values of 
log V"(R) are plotted against log R, for both m-mesons (upper line) and 
p-mesons (lower line). 


17 rg 21 23 Logk 25 27 2-9 


Figure 2. Plots of log N(R) and log Sy (R) against log R. The upper and lower curves. 
are drawn through the results for 7- and u-mesons respectively. 


§5. DETERMINATION OF m,/m 


1 

The determination of the quantity m,,/m,, can be carried out by any one of the 
following methods, the first of which was suggested to us by D. H. Perkins. It 
may be emphasized that they represent merely alternative modes of analysis, and. 
not independent determinations. 


(a) Perkins’ method 


It follows from the equation N = mF(R/m) that if we establish a graph showing 
the relationship between .’"(R) and R for u-mesons, we can deduce the corre- 
sponding curve for 7-mesons for any assumed value of the mass-ratio Mp, [1,. 
For each point (;(R), R) in the graph for u-mesons we construct the correspond- 
ing point W7"(R)m,/m,, Rm,/m,). The locus of these points then corresponds to 
the relation between 7(R) and R for z-mesons. Alternatively, if we plot the 
relation between log.’"(R) and log R, the corresponding values of log A’”"(R) 
and log Rk are given by a simple translation of magnitude log (m,,/m,) in both 
coordinates. 

Conversely, if the relations between (log. /’"(R), log R) and (log. V"(R), log R) 
have been established experimentally, values of m,/m,, can be deduced by drawing 
a succession of lines at 45° to the axes and measuring the translations, AB in figure 2, 
to which they correspond. The values so obtained are shown in table 4. 
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Table 4. Values of m,/m, deduced from normalized grain-counts 
Range (microns) 50 100 150 200 250 300 350 400 
Mz/ My 1-70 1-67 73 1-63 1:58 Ikoe ys) 1-60 1-64 


** Probable error’ (--) 0:18 0-14 0-13 0-13 0-13 On? Oat? Git 
The various values are not independent, and that with the greatest statistical weight is 
1-64+0-11. 
The values of m,/m, for the different values of R, are not independent, and 
that with greatest statistical weight is m,/m, = 164 +0-12, the limits corresponding 
to the magnitude of the “ probable error”. 


(b) Determination of the form of the function N=mF(R/m) 


It will be seen from figure 2 that the relationships between log.”"(R) and 
log R and between log ."(R) and log R are linear for values of R greater than 
100 » to within the standard deviations of the measurements. We can therefore 


write 
Nea REI Iki) Ga vOrr hae Ri eR a (2) 


We have dete1mined the value of m from the best straight line drawn through 
the points representing the measurements on z- and yx-mesons in figure 3, and have 
obtained the values 7,=0-701+0-014, ,=0-721+0-014, and the mean value 
70 aise OLS “Equation: (2) can itherefore be written 4 =kmo22 Ro 
From the position of the straight lines in figure 2 we find that km)?" =3 95 + 0-05 
Bnei = 4759 0-0/4 S0) thatey == 4-53 RO land 3-95 RO a re 
values of 1%” and 4)" calculated from this formula for different values of R 
show no significant difference from those based on the grain-counts given in 
table 3, except for values of R less than 100 ,. 

Alternatively we can plot the quantity a=./'/R” from the observed values of 


A™: a should be a constant for the results for a given type of particle 
C= i Rhyl i 


The values obtained are shown in table 5 and in figure 3 for both z- and u-mesons.. 


Table 5 Table 6 
Range ere) soe we CR) Individual values of the ratio m,/m,, 
(microns) ““— pom =“ ~ “porn 
50 .3°7340-10 4:10-+0-10 Event m|My Error Normalized 
100  3:94-40-08  4:39-+0-08 No. (%)  “‘ weight” 
150 3:924+0:07 4:55+0-08 2 2:06 +0-36 17-4 1-24 
200 3:95+0:06 4:57+0-08 14 1:52+0:27 17-6 (bowel 
250 3:98+0:06 4:58+0-08 8 2:04-40-45 YT 1:19 
300 3:95+0-05 4:58+0-08 12 1-07=-0°19 18-0 (Lois, 
350 3:97+0:05 4:50+0-08 5 Ie SSEOsIS 18-7 1-07 
400 3:9440:05 4:52+0-08 17 1:-45+0:28 19-3 1-00 
450 4-52+0-08 8 1-84+0-38 20-5 0-89 
500 4-53+0-08 9 1:65-+0-36 21°8 0:79 
550 4-55 +0-08 1 1237220300 22-0 0:77 
600 4-51+0-08 4 1:82+0°-42 23:3 0-69 
650 4:49 0-08 Weighted mean: 1:66+0-11. 
Weighted means: a,=3-95+0°05 ; 
Qn = 4°53 40-08. 


In this case, also, the different values are not independent, so that the “‘ probable 
error” is not reduced by taking the mean. 
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From equation (3) we can write 
a, a, = [m.,/m,]°?89 = 1-145 + 0-024; m,,m, = 1-62 + 0-12, 


We can calculate values of .””, corresponding to the observed values OLN is 


from the formula N = km289R°71, assuming that m,/m,=2-1, 1-65, 1-32 in turn. 
The corresponding values of the quantity ./”"/R" are represented by the dotted 
lines in figure 3 : (a) corresponds to the case in which the 7-meson disintegrates to 
produce a neutral meson, equal in mass to the p-meson, while (c) corresponds to 
the assumption that the momentum balance in the p-decay is provided by an 
emitted photon. 


§6. INDIVIDUAL VALUES OF THE RATIO m,/m, 


By using the equation N=km°?8°R°71 we can determine a value of m,/m, 
from the observations on the individual examples of the y-decay, and the values 
thus obtained are summarized in table 6. The results are also represented in 
figure 4. In this diagram the size of the oblong, corresponding to a particular 
determination, has an area proportional to its statistical weight. It will be seen 
that of the ten individual determinations, only one has a magnitude less than 1:32, 
and not one greater than 2-1. 


§7. THE EXISTENCE OF A NEUTRAL MESON 


The three methods which have been employed in the analysis of our obser- 
vations are consistent in that they show no significant differences from the mean 
value 1:64 +0-11. Although the fluctuations of the individual values are so large 
that we cannot regard our results as final, they strongly suggest that the momentum 
balance in the process of u-decay is provided by a neutral particle and not by a 
photon, and that this particle is less massive than the y-meson. .The application 
of the conservation laws of energy and momentum allows us to calculate the mass, 
m,, of this “‘neutretto’’, and we thus obtain the value m,=115 + 30 m.. 

Recent photographs by Anderson, Adams, Lloyd and Rau (1947) show that 
the cosmic-ray mesons brought to rest in the gas of an expansion chamber decay 
with the emission of an electron of energy 25 Mev. Assuming the mass of these 
mesons to be 200 m,, this result is consistent with the assumption that a neutral 
particle with a mass of about 130m, is ejected in the opposite direction to the 
electron in the decay process. ‘The agreement, within the limits of experimental 
error of the masses of the neutral particles indicated by the two experiments, is 
suggestive. We have therefore considered the consequences of assuming that 
the two particles are of the same type, and that the ordinary cosmic-ray mesons are 
identical with our z-mesons. ‘The decay sequence that follows from these assump- 
tions is represented in figure 5. Since the electron has half-integral spin we can 
conclude from the process 1>e+v that ,- and v-mesons differ by half-integral 
spin. It then follows, from the process 7 > p + v, that 7-mesons have half-integral 
spin. 

§8. THE DECAY CONSTANTS OF MESONS 
(a) m-mesons 


The distribution in range of the u-mesons, of which details are given in table 1, 
are represented in figure 6(a). Following Livingston and Bethe, we can define 
a. straggling coefficient, s, by the relation s=/(7/2)./{2,(R; —R,,)*/n}, where 
R; is the range of the 7th p-meson, the observed mean value of all the 1 j-mesons 
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being R,,. The measurements shown in table 1 are thus found to correspond 
to a value of s given by s/R,, =6-24%. 

Only extrapolated values of the range-energy relation are available, for 
particles of high velocity and charge |e|, in the ‘‘ Nuclear Research” emulsions. 
These are given with sufficient accuracy by the equation E =0-262M25R0575, 
where EF is the energy in Mev. and R the range in microns of a particle of M mass- 
units (proton=1). It follows from this equation that 5£/E =0-5755R/R. 
Writing 6R/R=s/R,,=0-0624 we find, that SE/E=0-0359. This value does not 
differ significantly from that corresponding to the width at half-maxima of the 
distribution in energy of the y-mesons shown in figure 6(b), which has been 
calculated from the observed range distribution on the assumption that the mass 
of the u-mesons equals 200m,. The new observations therefore confirm the 
previous conclusions that the u-mesons are emitted with values of the velocity 
which are constant within narrow limits, for the observed straggling coefficient 
is of the same order of magnitude as that suggested by the classical approach to the 
problem. 

Now the directions of ejection of the «-mesons are oriented at random with 
respect to the line of motion of the j.-mesons at the extreme end of their range. 
Let the difference in the direction of motion of the two particles, at the point of 
juncture of the two tracks, be 6. If the 7-mesons come to “rest” before decaying, 
this angle is without significance, for the particles will suffer frequent changes 
in direction, due to scattering, when moving at low velocities at the end of the 
range. If, however, a z-meson decays in flight, we should expect the velocity of 
ejection of the x-meson, relative to the emulsion, to depend on the value of 6. 
In these circumstances, the range of the ejected mesons would be very sensitive 
to the direction of emission relative to the line of motion of the parent 
m-meson. 

We have calculated the magnitude of the effect to be anticipated, making the 
simplifying assumption that the 7-mesons always decay when 5 from the end 
of their range. It is easily shown that in these circumstances the distribution 
in range of the p-mesons would extend from about 200, to 1200. Further, we 
have examined the distribution in range of the p-mesons in these events for which 
6 <z/2, and, separately, those for which 6>7/2. It is found that there is no 
significant difference between the two distributions. We therefore conclude that 
the z-mesons, in the observed events, were at or very near the end of their range 
when they decayed. 

We cannot conclude, however, that the observed events represent all the 
7-mesons which, having decayed in the emulsion, led to the production of a 
y-meson which stopped in the emulsion. If, for example, a 7-meson decayed 
when its residual kinetic energy exceeded 2 Mev., it would be difficult to distinguish 
the event from one in which a meson is scattered by a nuclear collision. We 
believe, however, that if the residual range of the 7-meson were less than 1004 
at the instant of decay, we should recognize the event because of the sudden 
decrease of the grain density in the track immediately after the point of “ scatter- 
ing”. We therefore conclude that in twenty cases a 7m-meson traversed a distance 
of 100 at the end of its range without decaying, and that in no case did it decay in 
this interval. The time of flight of such a meson is 3 x 10~™sec., and we can 
therefore write 20 x3x10-!2xk,<1, where k, is the decay constant of the 
m-mesons. It follows that k,<1:6 x 10'° and 7,>0-4 x 10-1°sec. 
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In every case in which the jz-meson reaches the end of its range in the emulsion 
it stops without leading to a disintegration with the emission of heavy particles. 
"This suggests that the - and p.-particles are positively charged —for the Coulomb 
repulsion would then prevent their close approach to a nucleus when moving with 
low velocities, and their interaction with it—and we shall assume this to be the 
ease:* 

We have now observed 150 events in which a z-meson, at the end of its range, 
produces a p-meson. In every case we can distinguish no displacement of the 
line of motion of the ejected p-meson from the last grain in the track of the 7-meson 
(see figure 7). We regard this as evidence that in no case did the 7-meson, at the 
end of its range, diffuse a distance greater than 2. before decaying. If it is 
assumed that a positively charged positive meson diffuses through the solid 
material of the emulsions as through a heavy gas, an estimate can be made of an 
upper limit to the half-period value r,. The result so obtained is +, >5 x 107“ sec. 
Dr. Frohlich has pointed out to us, however, that the problem of the motion of a 
positively charged meson in a solid substance is one of great complexity. It is 
likely that a positive meson will be “trapped”? wh:n brought to rest in a solid 
and that its diffusion will thus be prevented. We can therefore attach little 
weight to the limit obtained by this method. 


(b) p-mesons 


We have seen that in thirty cases the »-mesons reach the end of their range in 
the emulsion. In no case have we observed the track of a secondary p-meson to 
terminate in the emulsion without having the characteristics associated with that 
of a particle at the end of its range. We can therefore write for the decay constant 
k,, of these mesons 30k, x 107-1, 1.€.°k, <3:3 10? ye 2 10 Sec ire 
lower limits to the values of 7, and 7, derived from these observations have the 
following important consequence. ‘They indicate that a meson, of either type, 
with an energy of the order of 4Mev., will be unlikely to decay in free flight, zm 
vacuo or in a gas at normal pressures over a length of path p of at least several 
millimetres. Thus for .-mesons we have p, >8:7 x 10° x 2 x 10-!cm. =1-7cm., 
and for z-mesons p, >6:7 x 10° x 0-6 x 10-!°cm. =4mm._ It is therefore reason- 
able to undertake magnetic deflection experiments which demand for their 
success that the particles should not suffer spontaneous decay in a length of path 
of the order of 1 or 2cm. (Powell and Rosenblum 1948). 


(c) The postulated neutral meson 


The observed degree of homogeneity of the range distribution of the z-mesons 
produced by the decay of 7-mesons at rest indicates that the kinetic energy of the 
ejected particles £), is subject to a variation of less than 0-15 Mev., or 2-3 x 10-“ ergs. 
Further, the lower limits to the lifetime of the z- and y-mesons show that the 
energies to which their rest-masses are equivalent must be constant to within 


* These observations do not, of course, provide conclusive evidence for the sign of the charge 
for the “-mesons may have a very weak interaction with nuclei. Further, although evidence 2 
accumulating in support of our original suggestion that they are identical in type with the mesons 
‘commonly observed in experiments with expansion chambers and counters, there is no decisive 
‘evidence available to show that negatively charged mesons of this type, when brought to rest in 
materials of high atomic number, produce disintegrations with the emission of heavy particles. All 


that is known definitely, at the present time, is that they disappear without giving rise to a delayed 
coincidence. : 
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10“ ergs. Any variation, 5£,, in the energy of ejection of the n-meson, of the 
order of 10~’ ergs, would be associated, therefore, with a variation in the mass ‘of 
the postulated neutral meson. It follows that the observed degree of constancy 
of the values of E,, allows us to set a lower limit to the half-value period of the 
neutral meson 7,. We write 5E.8t~h, where A is Planck’s constant. Since 
SE <2:3 x-10-’ ergs, 7,>2-6 x 10-7! sec. 
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ABSTRACT. A study has been made of scattering of 40 protons and 160 mesons in photo- 
graphic emulsions in order to determine the mass of the individual particles. The spread 
in the values so obtained is large, but the evidence suggests that the majority of mesons 
recorded by the emulsion can be identified both with the mesons, of mass 200 me, commonly 
observed in cloud-chamber experiments and with counters, and with the jz-mesons observed 
in the photographic plates. The slow particles producing nuclear disintegrations, c-mesons, 
appear to contain a large proportion of particles with a mass equal to that of the 7-mesons. 


SNE ROD WiC ETON 


REVIOUS work (Bose and Chowdhry 1941, Perkins 1947, Occhialini and 
Powell 1947) has shown that some of the characteristics of the particles 
observed in photographic emulsions can be determined by an inspection of 

the “grain-density”’ and multiple Coulomb scattering associated with the tracks 
caused by the passage of these charged particles. Such methods can be ex- 
tended to give quantitative information about the masses and charges of the 
particles. 

Grain counting is especially favourable for the determination of the ratio of the 
masses of particles involved in contemporaneous events, such as the 7-1 meson 
decay process (Lattes, Occhialini and Powell 1947); but it is not well suited for 
absolute mass determinations, for reasons discussed by the authors. 

The method based on observations of multiple Coulomb scattering is indepen- 
dent of the fading of the latent image and gives absolute values. Although 
inherently less accurate than grain counting, it is therefore more suited for work 
with particles occurring in events of which the relationship in time is unknown, 
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The theory of multiple scattering has been developed by Williams (1939, 1940), 
and several experiments have been made to verify the formulae which he estab- 
lished. Probably the most accurate and complete is that of Kulchitsky and 
Latyshev (1942) on the scattering of fast electrons by light and heavy elements. 
These authors found agreement with Williams’ formula for light elements, but 
noted systematic deviations for the heavier elements. These deviations can be 
explained by the more accurate calculations of Moliére (1947) on the small angle 
single scattering by a screened ‘Thomas-Fermi atom. 

It would appear, therefore, that the calculations of Williams and of Moliere 
forma reliable basis for determining masses by observations on multiple scattering. 


§2) THEORY OF METHOD 

A charged particle traversing a medium suffers frequent small-angle deflections 
in elastic collisions with the nuclei of the atoms composing the medium. As a 
result of large numbers of these small deflections, the track has a curved appear- 
ance. ‘This scattering is purely statistical in nature, and Williams has evaluated 
the probability of a change of direction, «, occurring in a given medium, for a 
path length ¢, from the distribution of the values of the angles of scattering occur- 
ring in single collisions. Williams’ theory shows that the probability distribution 
for small values of « is very nearly Gaussian; for large values the distribution 
approaches that given by the Rutherford formula for single scattering. The 
deflections which he considers are the projection of the true deflections on a two- 
dimensional plane. In our experiments, also, we define « as the two-dimensional 
projection of the true deflection. 

The arithmetic mean value of the deviations considered, <«> yy, is given by a 
formula of the form <a> sy={f(8, t, N, Z)}/E, where B=ratio of velocity of 
particle to velocity of light, t=thickness of scattering medium traversed by 
particle, N=number of atoms per cm® of the scattering medium, Z=charge 
number on nucleus of atom of scattering medium, E=mean kinetic energy of 
particle in element of trajectory under consideration. 

a varies with the energy and accordingly with the residual range of the particle. 
It is convenient to define a quantity, e, which is independent of range and which is 
given by the relation «=«E/f(,t, N, Z). The values of ¢ will have a probability 
distribution similar to that for «, and <e>,y=1. 

Consider the trajectory of a charged particle, of unknown mass, moving in a 
given medium and undergoing successive deflections «,, «;, 73; etc. If we assume 
the particle to have a mass M, we can, from the range-energy curve for protons, 
calculate the velocity of the particle for any value of the residual range. It is then 
possible to calculate the quantity <e> ,y for the trajectories of particles of this 
type; and similar calculations can be made for any assumed value of the 
mass. ‘The value of M, which enables us to satisfy the relation <> ,y =1, for the 
trajectory of a particular particle, then gives a measure of its mass. 

M can be determined analytically in the following manner. We evaluate 
<€,.> av for an approximately correct value m of the mass. 

Now 

2 Ee ipaeea 
- En (Pu, t, N, Z) 
H(B,t, N, Z) is relatively insensitive to a change in mass, in comparison with E. 
We can therefore write ey, —~e,,.Ey,/E,,. The relationship between mass 
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and energy may be calculated from the range-energy data for photographic 
emulsions (Camerini 1947). This gives E,,/E,, ~(M/m)4, since the range 
exponent is 0:56 for Tees of charge |e]. Therefore ¢,,-~(M/m)°: Ae 
and <ey>sy=(M/m)**<e,,> ay- But by definition <<y,>,y=1, -so that 
M=m\ <em> av} *™. 

The values of ¢ are distributed according to a probability function and the 
mean, <e> y, is a Statistical quantity. It is therefore important to decide 
whether its value is best determined from the arithmetic mean, mean square or 
median of the individual values. An analysis shows that the arithmetic mean of 
all values of « below a certain cut-off point should be determined in order to give 
the least probable error. ‘This mean, which we shall denote as <egg> yy, is the 
arithmetic mean of all values of ¢ less than 4<egg>,y. The determination of 
<eco> av depends upon an approximate knowledge of the answer; in practice the 
problem is rapidly solved by a relaxation method. The necessity for ignoring 
terms greater than 4<eco> ay arises from the fact that very large deviations due 
to single collisions may occur with finite probability ; and these would, if included 
in the statistics, give rise to very large fluctuations in mass values. 

The distribution of the individual values of « is very nearly Gaussian. 
From an examination of this distribution, for « up to the cut-off point, the 
probable MESS error in the determination of <eog>,y was found to be 

+ (0:55/4/n) <eco> ay, where n is the number of independent observations made 
on the deflections. When x has large values, this leads to a probable error in the 
mass determinations of + 123/,/n %. 


$3, EXPERIMENTAL 


After attempts had been made to measure the scattering along a track in a 
photographic emulsion by means of a normal high magnification microscope, it 
was decided to use the mechanized projection microscope (Occhialini, Payne and 
Powell 1948) for this purpose. Large-scale drawings were taken with this instru- 
ment. The advantage thus gained over direct measurements were threefold : 
firstly, the spurious scattering on the track was reduced (see below) ; secondly, the 
residual spurious scattering could be accurately evaluated ; and lastly, a permanent 
facsimile of the track at high magnification was available for inspection and analysis. 

In order to facilitate the actual drawing, a magnification of x 4000 was used. 
This magnification was obviously too large relative to the resolving power, but it 
allowed us to obtain drawings on a convenient scale for analysis. Only the centre 
of the field of the projection microscope was used in making the drawings. 

The drawing of the track was then divided up into a series of segments, the 
lengths of which were chosen to fulfil the conditions indicated below. ‘The changes 
in direction between the middle points of alternate segments were then measured. 
The distances between these middle points are the cell lengths to which 
reference is made below. 

There were two possible methods of measurement: of angles either between 
successive tangents, or between successive chords (‘‘sagitta””’ measurements are a 
special case of chords). Since the track is defined by a small number of grains, it 
is obviously impossible to draw accurate tangents. ‘The most accurate procedure 
is to introduce deliberately a small systematic error into the measurements by 
measuring angles between mean lines drawn through grains over certain segments 
of a track. The length over which this mean line is drawn is chosen as half the 
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cell length (figure 1). The mean line will generally be intermediate between a 
chord and tangent. The theory of a multiply scattered track has been given by 
Rossi (1941) and leads to an easy calculation of the systematic error, as shown in 
Appendix 1. It is found that measured angles must be increased by 442%, to 
correspond with those in Williams’ theory (Appendix 2). 

We have found it possible to increase the effective number of measurements on 
a single track by the following procedure, illustrated in figure 2. ‘The track is 


A. Method of measurement of angular deviations. 


; Pp Meson F/13B 


B. Reproduction of section of drawing made on typical meson track, showing 
lines drawn through grains in successive segments. 


Figure 1. 


Figure 2. Overlapping half cells drawn on track. 


divided into overlapping half cells, 1, 2, 3, etc. In each of these a straight line is 
drawn to pass, as closely as possible, through the centres of gravity of the grains. 
The angular deviations between the lines in 1 and 3, 2 and 4,3 and 5, etc., are then 
measured. It can be shown that for a large number of observations of this type 
two-thirds of the total number of measurements are statistically independent. 

The cell length used must be such as to keep the experimental errors low com- 
pared with the expected scattering. The experimental errors arise from the 
‘finite size of the grains, the distortion of the microscopic field, the errors in drawing 
and the errors in measuring. These effects give a “spurious scattering’ 


” super- 
imposed on the Coulomb multiple scattering. 


Protons: Number of half cell 
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The spurious scattering has been evaluated by means of observations on the 
tracks of very fast particles (protons of 40Mev., one «-particle of 120 ev.). 
These tracks have a negligible multiple scattering, and, in the case of protons, 
the same grain density at corresponding points in the tracks as mesons of the same 
velocity. Almost all the changes in direction observed in measurements on these 
tracks are due to the sources of error indicated above ; they correspond to spurious 
scattering. ‘The results of these measurements are given in Appendix 3. They 
show, as would be expected, that the spurious scattering decreases with longer 
cells and high grain density. On the contrary, the Coulomb multiple scattering 
increases with longer cells. ‘The length of the successive cells can therefore be 
chosen to keep the spurious scattering less than one-quarter of the true multiple 
Coulomb scattering. 

The values chosen are: 


1 
Length of half cell (jx) (ey I 
Lengths of successive | 
overlapping cells (2) ip 


Mesons : Number of half cell 1 2 3 + 5 6 Wt 8 10) 11 iy 


Length of half cell (2) LO NO alo key IS gy I aI aS DO) 


Lengths of successive 20 20 22:5 275 30. 30 30 30 32-5 37:5 40 
overlapping cells (x) 


§4. RESULTS 


As a check on the method, measurements were made on two series of protons. 
‘Twenty proton tracks (each 300 microns in length) were selected at random from 
a plate exposed to the disintegration particles arising from the reaction B1°(d, p)B™. 
‘The mean value obtained for the experimental mass was 1-10 +0 10 proton masses. 
This value depends upon a knowledge of the atomic composition of the emulsion, 
and we have used the values given by Messrs. Ilford Ltd. ‘The limits of error 
indicated correspond to statistical fluctuations alone. There is an additional 
uncertainty, due to the inherent errors of the experiment, equal to about +0-07 
proton masses. 

_ Similar observations on twenty proton tracks observed in cosmic-ray plates 
gave a value 1-20 +010 proton masses. ‘There is no significant deviation from 
the expected value, particularly as we cannot exclude the possibility that tracks due 
to deuterons or tritons may have been included among those selected for measure- 
ment. Histograms of the distribution of the individual values of the mass are 
shown in figure 3 by full lines. ‘The dotted lines show the corresponding distri- 
bution calculated for a gamma distribution and normalized to correspond to the 
same number of tracks. A comparison of the two distributions indicates a 
satisfactory internal consistency between the observations. ‘The mean <€co> av 
was cross-checked for constancy with range and varying cell size by taking the 
values of « at the same residual range for all the particles of the same class, and then 
determining <€c9>,y for these terms. ‘The values of the means thus obtained 
show no significant variations. 

The tracks of particles on the cosmic-ray plates (Ilford, Boron-loaded C.2 
exposed at 2800m. and 5500m.) with grain density less than those of protons 
were divided into four groups. The 7- and yz-mesons were taken from the 7-1 
decay process. Some 7-mesons were taken in which the jz-meson resulting from 


ee 
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the decay did not end in the emulsion. In these cases, however, only those 
m-mesons were accepted for measurement of which the track of the associated 
y-meson was of sufficient length to allow a definite identification. Similarly, a 
y-meson track was accepted for measurement only if it was accompanied by the 
track of the primary 7-meson, and if the track of the y-meson of length of the order 
of 600 microns terminated in the emulsion. Following the previous definition 
adopted in this laboratory (Lattes, Occhialini and Powell 1947), we distinguish a 


Scattering Units Scattering Units 
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N°: of Protons 
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75000 5000 -2000-+4000 70000 ~~ 19000 5000 2000 (1000 700 500 
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A. 20 protons from H.T. set. B. 20 protons from cosmic rays. 


Figure 3. Histograms of the distribution of values for proton masses, in scattering units and 


mass units. The broken line represents the expected distribution ; full line represents 
the observed distribution. 


third class of mesons, o-mesons, which suffer capture and lead to observable 
nuclear disintegrations in the emulsion. We can also define a fourth group, 
p-mesons, which stop in the emulsions without giving rise to any observed charged. 
particles. The origin of these mesons is uncertain, but most of them can probably 
be identified as the mesons commonly observed in experiments with Wilson. 
chambers and counters. 

Combination of the scattering results with previous grain-density measure- 
ments permits an unequivocal assignment of charge |e| to all measured particles. 
If the charge were greater, the results would lead to a gross inconsistency between 
the masses as measured by scattering and by grain counting. 

The results of the mass determinations on 207-, 20 u.-mesons, 60o- and 60 p- 
mesons are shown below together with the “ probable” statistical error. 

Mass of z-mesons: 260 +30m,, 
Mass of p-mesons: 205 +20m,, 
Mass of o-mesons: 275 +15m,, 
Mass of p-mesons: 200 +10m,. 


N° of Mesons 


N° of Mesons 
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The systematic experimental error is again of the order of +7°4. Histograms 


similar to those for the protons are given for these four classes in figure 4. 


hte histogram for the p-mesons shows a slight tailing at the low-mass end, two 
mesons lying outside the expected histogram. The statistical errors are, however, 
too great to allow us to regard these observations as evidence for the existence of 
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Figure 4. Expected (broken line) and observed (full line) histograms for mesons. 


mesons of mass approximately 80m,. ‘The distribution in «€ is not “normal” 
as the method of taking overlapping cells gives rise to fluctuations in the number 
of independent readings, and the imposition of a “cut-off” will also give rise to 
occasional large values of <egg>,y- The theoretically expected histogram was 
calculated by taking the normal frequency distribution with the same probable 
error. In practice for the above-mentioned readings a slight tailing of such a 
distribution should occur. ‘The results otherwise appear consistent with the 


expected probable error. 
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§5. CONCLUSION 


It has previously been suggested (Lattes, Occhialini and Powell 1947, Marshak 
and Bethe 1947) that the constitution of the z-, u-, o- and p-mesons is as 
follows : . 

The 7-mesons are heavy mesons with masses of 300-420 m, and strong nuclear 
interaction. ‘The p-mesons are identical with those hitherto observed in cloud 
chamber and delayed coincidence experiments with approximate mass 200m, 
(Fretter 1946, Hughes 1947, Valley 1947). The p-mesons were supposed to be 
predominantly j.-mesons with a slight admixture of z-mesons for which the decay 
had escaped observation, and the o-mesons to contain a high proportion of 
negative 7-mesons. 

Our results give support to these assumptions. The mean mass of the 
o-mesons is significantly higher than that of the p-mesons. ‘The masses of the 
y- and p-mesons are in agreement, within the probable error, with the results of 
Fretter’s mass determinations at sea level (mass 202 +5 m,). The best value of the 
mass ratio of 7- to u-mesons from grain counting is 1-66 + 0-11 (Lattes, Occhialini 
and Powell 1948). There is no significant disagreement with the mass ratio 
obtained from our experiments of 1-3 + 0-2. 


APPENDIX! 
Correction for systematic errors 


The number of useful independent measurements which can be made on the 
track of a particle is controlled by two conflicting factors. ‘To minimize the 
statistical errors, it is desirable to make the largest possible number of independent 
observations on the scattering, and therefore to choose the smallest possible cell 
length. In practice, however, the effects of spurious scattering impose a lower 
limit to the length of the cell which can be employed in a given region of the track. 
The error introduced in the evaluation of <eco> ay is the result of a compromise 
between the errors due to statistical fluctuations and to measurement. 

In conformity with these considerations, cell lengths along the track are chosen. 
so that at each point the spurious scattering does not exceed a known small fraction 
of the multiple scattering. Within each cell the deviation is measured with the 
least possible error by drawing the mean line through the grains contained in 
successive segments of the track and measuring the angular change between the 
directions of the track in these segments. A limit is set to the size of the segment 
by the magnitude of the systematic errors introduced. These systematic errors 
arise in the following way. It can be proved quite generally that the arithmetic 
mean of the angular deviations between successive tangents along the multiply 
scattered track is larger than the arithmetic mean of the deviations between 
successive chords. The mean line drawn through consecutive grains will in 
some cases correspond more closely to a tangent and in othercasestoachord. The 
arithmetic mean of such a series of measurements, <x» > ay, Will have upper and. — 
lower limits set by the values for the means of tangents and chords. The calcu- 
lation of the limits is most easily carried out by the methods given by Rossi and 
Greisen (1941). 

Consider figure 5. Chords AB and CD of length kt and tangents A’B and CD’ 
are drawn. The cell length is t; kt is a fraction of this length. We require to 
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find the mean value <a> ,y of the angle BA—CD. Now 
| %exp = BA —CD =(BA —B’A) +(B’A —CD’) -(CD —CD’), 
and, therefore, 
(Oexp> av=<(BA —B’A)®» av + <(B’A — CD’) ay + <(CD - CD’)? > av- 


The cross-terms vanish 
in the average as they con- 
tain positive and negative 
contributions of equal mag- 
nitude. Rossi and Greisen 
(1941) showed that 


<(ED -CD'®> ay 
=((BA—BA’)®> ,y = dake, 
where a) is the mean square of angles between successive tangents at intervals of 
one unit of length. Similarly <(B’A—CD’)?>.,=«a2t(1 —k). Substituting 
these values in equation (1), <«?> yy =aét(1 — 4h). 

As was explained above, <> ,y has upper and lower limits corresponding 


to the mean deviations between successive chords and tangents. The mean 
square deviation for the tangent is by definition «5, and, therefore, 


V (2/7) toV/t > exp» av > V (2/m)atovV/t(1 — BR), 
since the transformation from mean square to arithmetic mean involves a constant 
which is the same throughout the terms. 
The most suitable value for k from the experimental data available was taken as 
k=}. For this value «/(2/m)a9v/t><Gexp> av > V(2/7)0-92a94/t, and, therefore, 
<Gexp> av = 0:964/(2/m)av/t + 0:02 \/(2/m7)a9+/t. 


Similar  calcula- 
tions were made for 
sagitta with smooth- 
ing as shown in figure 
6. This treatment 
leads to difficulties in 
calculation of the 


scattering distribution 
function; the connection between tangents and sagitta was only derived by Rossi 


by effectively ignoring the single scattering “tail”. | At the same time, the accuracy 
of the measurement is not improved. The method described above was therefore 


adopted in preference to sagitta. 


Figure 5. Correction for systematic error in measurement. 


Figure 6. Sagitta measurement on smoothed track. Vertical line 
indicates sagitta distance. 


APPENDIX 2 
Application of Williams’ theory to mixed media 


The extension of Williams’ theory to media consisting of a mixture of atoms is 
comparatively simple. Williams gives the following summary of his method. 
The units used are units where, 8 = 2Zze?(Nt)*/(M?c?£); Zis the atomic number 

-of scattering nuclei; ze, M and fc are charge, mass and velocity respectively of 
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the scattered particle, €=(1 —f?)!; N is number of atoms per cm; and f is 
thickness of scattering material. 
The arithmetic mean scattering in units of 5 is given by Williams as 


<a> py =0-80< a, > + 1-45, 
Co Ee 2=2In (Prac! Pann: 


Williams defines ¢,,,, a8 equal to 4/(7/2) 6-units and ¢,,;, aS a characteristic angle 
for the atoms of the media, given (by Williams) for low and high velocities of the 
incident particles. Moliére (1947) has calculated 4,,;, for all velocities. 

Before Williams’ theory can be applied to mixtures, the following changes must 
be made, 6 being defined as 


PB wapt NJ #\ 4 
Ss sy tix Nit 2 
; M?B4¢4e2 ? 
where the subscript 7 refers to the different molecular species. <«,,>? now 
becomes 2 
Q 2,4 uN,Z; In Karey) (dmin), | 
OSC t 
Me pita 82 


As before, @ma.=4/(7/2) 6-units 


, 
and (¢ in), refers to the characteristic i = 
angle for each molecular species 7 in 33 ae 
the medium. | 

For the purpose of finding mc 
<aeo> ay With a large angle cut-off, @ 31 = 
the distribution inw is wellrepresented 5 
for large values of « by ine 

(z7/a°)(1 a 3m <a> |”). ai | pee | | | | 

In figure 7, curve A gives the value of eaees Ve | ia 
<G ays) and curve Bagg as both k fo df | ee | 
measured in terms of 8-units, as a 2 ei re 60700 
function of cell length in the Ilford C.2 Cell Lene imer ee! 
Boron-loaded emulsion. In curve B Figure 7. The variation for different cell lengths 
all angles greater than 128 are cut off. of (A) the arithmetic mean scattering <%>av 
The value for 8 was taken as 1:564/t/E, in 6-units, (B) <«co>avy in similar units with 


iy i 5 angles greater than 12 6 cut-off. 
where ¢ is in units of 10 microns. nce 
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Determination of spurious scattering 


The main limitations on the accuracy with which the scattering of a track can 
be measured are due (a) to the finite size of the grains making up a track, and (0) 
to the introduction of experimental inaccuracies. ‘These effects produce spurious 
scattering, whose distribution is approximately a Gaussian error function. The real 
multiple scattering of the track also has an approximately Gaussian distribution. 
The experimentally measured scattering is therefore formed by the Gaussian 
addition of the real and spurious scattering <a>.» = <a> rent “t (O" > goa gels 
therefore, the quantity required is <a> sa =<@"> exp —<®"> spur, Which may 
be written as <a> real— <a> expt as 2 <ae> ate! a eens for <a> exp > <a> spur? 
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In order to determine <x, spur the following series of experiments was 
carried out. 

Drawings were made, with the aid of the projection microscope, of the tracks of 
12 fast protons (c. 40 Mev.) and one fast a-particle (c. 120Mev.). Sections of 
tracks selected for measurement had grain densities varying from 3 to 10 grains 
per 10 microns, and were traced out under a wide range of optical conditions and 
with wide variations of the density of the background grains. For these tracks 
> spur Was much greater than <x >,.,, and the measurements were effectively 
a determination of the spurious scattering. 

To determine spurious scattering two types of measurements were made. 
Those shown in figure 8(a) are the same as those previously described in the 


; Sa 
| Pepto te fr et Ooo pe SBS RR tp OT eal os eee te 


(a) Mean lines through grains of overlapping cells for fast particle. 


U | 

| yl] 0) | 
| ~ +t | 
los cee, 


(6) “ 'Triple-point ’’ determinations. 


Figure 8. Determination of spurious scattering. 


determination of the multiple scattering along the track. Those shown in 
figure 8 (b) consisted of so-called triple-point determinations in which the displace- 
ments, X, of the centres of gravity of the grains were measured, from a line 
running roughly parallel to the direction of the track, and 


<X> ay = (1/(m —2)) ~ | (Ges ae ep Ce A) | 


was measured. <X> qy is, therefore, a factor which indicates the mean displace- 
ment of the individual grains about the trajectory of the particle, and provides a 
useful indication of the spurious scattering on all types of track, even those with 
large multiple scattering. 

The results of the measurements show the following relationship between 
<X> sy in microns for a grain density of N per 10 microns, cell-length ¢ in 10 
microns units (as in figure 1), and (a> 7 =39<X> ay/V/N . t!?, measured in 
degrees. <X> y for individual tracks had values ranging from 0-125 micron 
to 0-25 micron. Measurements on 10 meson tracks showed the same limits. 
The cell length for measurements on mesons and protons were chosen so that 
PORRC Xa = O22 uMICLON <B> wpmel <> pos 

The correction for spurious scattering over a whole track was thus maintained 
at less than 3%. 
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A Modification of Benoit’s Method of Exact Fractions 


BY Es G:. MICHAEEIS 
Physics Department, Birkbeck College, London 


Communicated by R. Fiirth ; MS. received 3 February 1948, and in amended 
form 12 April 1948 


ABSTRACT. A process is outlined by which the computations arising in the application 
of Benoit’s Method of Exact Fractions to interferometric measurements can be shortened 
and simplified. A graphical illustration of the method is given, and it is compared with other 
modifications of Benoit’s original technique. 


§1. INTRODUCTION 


ORK carried out in this department has shown that Benoit’s ‘‘ Method 

of Exact Fractions’ as used in interferometry is capable of a slight 

modification, which has been found to shorten and facilitate the 
computations involved and which may also serve to extend the usefulness of the 
technique. 

Benoit’s method is commonly employed to find the order of interference of 
light between the plates of an interferometer from the approximately known 
distance between the plates and the measured fractional part of the order at the 
centre of the fringe-system. It was developed for use with the Michelson 
interferometer; nowadays it is widely used in work with the Fabry-Perot etalon 
and with such instruments as Késters’ interference comparator. Full accounts. 


have been given by Benoit (1898), Childs (1926) and Williams (1930). 
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§2. THE METHOD OF EXACT FRACTIONS 


For the purpose of the method the integral part p and the fractional part f of the 
order of interference m at the centre of the fringe system are determined separately. 
An approximate value p, of p is found from the separation of the plates as deter- 
mined by mechanical means and, in the case of the Fabry-Perot interferometer, f is 
obtained from the diameters of the interference rings. The true value of n, fora 
wavelength A, is then given by 


Pic le jhe ip) nee eo. (1) 
where m is a positive or negative integer to be determined. 
Now m has to be so chosen that the fractional parts of the central orders n 
corresponding to wavelengths A calculated by 


[REL A Oo i eer (2) 


agree with those found experimentally. Usually m is given successive integral 
values until this desired agreement is achieved simultaneously for all wavelengths 
used. This method, though simple enough in theory, may nevertheless involve 
a large number of seven- or eight-figure computations when more than four 
wavelengths are used for the calibration, and the assumed p, differs appreciably 
from the actual value. 

Methods for shortening this work have been developed by Kosters (1926) and 
by Pérard and Maudet (1927). Kosters describes a slide rule on which the frac- 
tional parts for various spectral lines are plotted against the separation, which may 
be read off when the required combination of fractional parts has been found. 
Pérard and Maudet show that, by a method to which the following js related, the 
plate-separation may be found in a series of approximations when the mere order 
of its magnitude is known. 


§3. BASIS OF THE MODIFICATION 

The method now suggested is based on the following consideration :— 

For m=0 we have, by (1) and (2),n=(p,+/4)Ai1/A=p + 4, say, where p and ¢. 
are once more the integral and fractional parts respectively. For another value 
of m, 

n=(pPytfrtm)r,/A=pt+d+mA/A=(p+m)+h+ mA, —A)/A. 

If, therefore, the fractional orders f for the wavelengths A are to agree with 
observation, then ¢ + m(A, —A)/A may only differ from f by an integer, 1. e. 

m(A, —A)/A=f—¢d+r, Fave l Dee lo terece (3) 
or m=(f—b+r)Af(Ay—A). nee es (4) 
The integer m to be added to p, must of course be the same for all values of A. 
Hence it is only necessary to solve (4) for a range of values of r for each wavelength 
until an m is found which is common to all. This will in general occur with a 
different r for each wavelength unless these lie fairly close together. Since m 
is a small integer, only a low degree of accuracy is required, and all auxiliary 
calculations may be carried out with a slide rule. 


$4. GRAPHICAL INTERPRETATION 


The method also permits of a graphical interpretation. By plotting m 
against r a straight line is obtained for each A. ‘The intersections of each line with 
r=0, +1,2,..., represent values of m yielding the correct f for the corresponding 


196 E. G. Michaelis 


wavelength. The m to be selected for a general coincidence is then charac- 

terized by such points of intersection which for all ’s lie on the same ordinate. 
The diagram shows the straight lines thus obtained in the measurement 

of the plate-separation of a 1-cm. etalon, using neon light. The following wave- 


x Partial near-coincidences. 


O General coincidence at m= 32. 


lengths were employed: A, =5852-4878a., A, =5881-8950 4 eserue 
Ay = 6029-9971 A., As = 6074-3377 a. : tg O77 342A 


Near-coincidences for A,, Ay and A; occur at m= —103, and for A,, d 
m= +99, but the correct valu - 
parts is clearly given by m= +32. 
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§5. CONCLUSIONS 


The method of calculating the correct value of m has been used by several 
students of this department and has been found to yield quick and reliable results. 
Using four wavelengths only eight seven-figure calculations are required, and it has 
been possible to carry out the entire work without the use of a calculating machine. 
Since the method is clearly capable of distinguishing the correct value of m from 
a range lying between + 100, the approximate separation of the plates need only 
be known to about 1/10mm. This fact might facilitate work where large plate- 
separations have to be used, e.g. in the determination of the metre. 

The process outlined would appear to represent a final step in the method of 
Pérard and Maudet; it also seems the only one necessary unless the conditions are 
as stringent as those assumed by these authors. The method of finding m cannot, 
of course, compare in rapidity with that made possible by the use of Kosters’ 
slide rule, but it will be recognized that the acquisition or the making of such an 
instrument will not be warranted unless a large number of routine measurements. 
with a standard source are to be made. 
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An Isotopic Abundance Rule 


F. C. Frank has recently (1948) drawn attention to the exceptions to the rule requiring 
even-odd nuclei to have at least one even-even neighbour of greater abundance. ‘To the 
exceptions which he quotes (1”°Xe, 17Sm, “Sm, 1*°Pt, 78°U) may be added °*Ru, which is 
slightly more abundant than 1°°Ru (Ewald 1943).. We should like to comment on these 
exceptions. 

"Ru, 2°Xe. The abundances of these isotopes have almost certainly been enhanced by 
the decay, during geological time, of })Tc (half-life 4 10® years) and of TANT respectively. 
It seems likely (see Katcoff 1947) that !°I has a sufficiently long half-life to have survived 
beyond the solidification of the earth’s crust. Assuming that the excess #?°Xe content of 
the atmosphere is derived from ™°I originally in rocks which have by now suffered erosion, 
one can make a rough calculation of the half-life of ?°I. The value we obtain is about 
4 108 years, in approximate agreement with the value suggested by Katcoff. Experiments 
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are in progress here to check this hypothesis by examining geologically-old iodine-containing 
minerals for xenon by a micro-gas analytical technique. These two exceptions are in fact 
attributable to the unexpected f-instability of the parent isotopes, rather than to any 
irregularity on the part of the daughter-substances themselves. 

M7Sm, 19Sm, 5U. These exceptions can readily be understood in terms of the 
a-instability of “8Sm and *4U. In these cases it is the decrease in the abundances of 
the neighbouring even-even isotopes, during geological time, which causes the anomaly. 

199Pt, No such explanation is available in this instance. It may be pointed out, 
however, that !°Ir has not been reported ; it may well be that it has a half-life too short to 
have survived to a significant extent, but still too long to give an observable activity in such 
bombardment experiments as have so far been carried out. 

K. F. CHACKETT. 

Londonderry Laboratory for Radiochemistry, G. R. Martin. 

University of Durham. 
5th May 1948. 


Ewa pb, 1943, Private communication to Fliigge and Mattauch, Ber. dtsch. Chem. Ges., 76, 1. 
Frank, F. C., 1948, Proc. Phys. Soc., 60, 211. 
Kartcorr, S., 1947, Phys. Rev., 71, 826. 


Distribution Coefficients for the Calculation of Colours in the 
C.L.E. Trichromatic System for a Total Radiator at 2450°x. 


Smith (1934) and Harding and Sisson (1947) have published tables of the distribution 
coefficients and the relative energy distributions for total radiators from 1500° to 3500° xk. 
in steps of 250° and for the equal energy stimulus and standard illuminants A, B,C. These 
did not include data for 2450° k. required by us in connection with the estimation of 
chromaticity coordinates and integral transmission of colour filters. 

The relative energy distribution and distribution coefficients for a total radiator at 
2450° have therefore been computed in the same manner and for the same wavelength 
intervals as those prepared by Harding and Sisson. ‘The use of ‘“‘ condensed ’”’ tables 
of this type having entries at every 0:01 micron reduces the time required in the calculation 
of chromaticity coordinates compared with that required for the use of the original C.1.E. 
tables which have twice the number of entries (Smith and Guild 1931-2). The previous 
tables have been calculated using a value of 14350 for the constant C, in the Planck formula, 
and this value was retained in the preparation of the tables here. 

The general method followed is that given in a recent paper by Harding and 
Sisson (1947). Relative energy values were calculated from the Planck formula 
Eyg= C,A~*/{exp (C2/A#)—1} to seven decimal places, the constant C, being eliminated 
by calculating all energies relative to that at a wavelength of 0:56 micron. From these 
energy values the distribution coefficients were calculated using the distribution 
coefficients for an equal energy stimulus published by Smith (1934). The values were 
then adjusted so that the sum of the (£jy,;) entries was 100 and the entries rounded off 
to four decimal places. 

The relative energy distribution and the distribution coefficients are given in tables 
and 2. 

The accuracy of the tables was checked by calculating the chromaticity coordinates 
of a total radiator at 2450° kK. and comparing ‘the colour equation with published data 
(Harding 1944). 

The figures were found to be in agreement, as is shown in the table below. 


Origin Colour equation (C2= 14350) 
Calculated 7-figure tables 0-4809261x+0-4141740y+0-1048998~ 
Table 2 0-480926x +0-414174y +0:104900z 


Published figures 0-480926x +0-414174y +0-1049002 


Table 1. 


Wavelength 
in microns 


0:38 
0-39 
0-40 
0-41 
0-42 
0-43 
0-44 
0-45 
0-46 
0-47 
0-48 
0-49 
0:50 
0-51 
Ory) 
0:53 
0:54 
0-55 
0:56 
OF5i7 
Table 2. 
W avelength E,x, 
‘m microns 
0:38 0-0005 
0-39 0:0024 
0-40 0-0104 
0-41 0-0392 
0-42 0-1591 
0:43 0-4046 
0-44 0-6040 
0-45 0-7012 
0-46 0-7214 
0:47 0-5740 
0-48 0:3229 
0-49 0:1264 
0:50 0:0219 
0:51 0-0482 
0:52 0:3686 
0-53 1-0834 
0-54 2:1292 
0-55 85307] 
0-56 5-3503 
0:57 7:5446 
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Relative energy values for a total radiator at 2450° k. 
calculated from Planck’s equation : C.=14350 


Relative 
energy 
4901883 
6°391445 
8:197502 

10-35542 
12-89910 
15-86034 
19-26807 
23:14781 
DVISSAN EY 
32-40599 
37:°81488 
43-75630 
50-23396 
57:-24681 
64:78936 
7228 5077 
81:42009 
90:47657 
100-0000 

109-9661 


Wavelength 
in microns 


0.58 
0-59 
0-60 
0-61 
0-62 
0-63 
0-64 
0-65 
0-66 


Relative 
energy 
120-3480 
131-1161 
142-2398 
153-6855 
165-4198 
177-4074 
189-6137 
202-0026 
214-5389 
227-1876 
239-9144 
252-6852 
265:°4676 
278-2301 
290-9424 
3035756 
316°1022 
328-4958 
340-7331 
352-7905 


Distribution coefficients for total.radiator at 2450° K. 


Research Laboratories of 
the General Electric Company Ltd., 
Wembley, Middlesex. 
26th April 1948. 


Wavelength 
in microns 


0-76 
0:77 


Ex; 


DOD Se) 
12-1188 
13-6079 
13-8791 
12-7058 
10-2816 

7:6341 

5:1653 

3-1789 

1:7946 

1:0010 

0:5198 

0-2689 

0-1465 

0:0746 
00383 
0-0199 
0-0089 
0-0046 
0:0016 


Totals : 116-1169 
Colour: 0:48093x-+-0:41417y+0-10490z 


Harvine, H. G. W., 1944, Proc. Phys. Soc., 56, 305. 
Harpine, H. G. W., and Sisson, R. B., 1947, Proc. Phys. Soc., 59, 814. 
SmitH, T., 1934, Proc. Phys. Soc., 46, 372. 
SmitH, T., and GuILp, J., 1931-2, Trans. Opt. Soc., 33, 73. 


Ey, 


A232 
O59873 
80789 
6:9630 
5-6661 
4:2409 
2:9828 
1:9484 
1:1761 
0:6575 
0:3639 
0:1877 
0:0968 
0:0526 
0:0262 
0:0137 
0-0085 
0:0029 
0-0015 
0:0000 


100-0000 


2o9 


Ez, 


0:0173 
0-0136 
0-0102 
0:0048 
0:0022 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
0-0000 
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SytviaA A. E. Hirp. 
H. F. STEPHENSON. 
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REVIEWS OF BOOKS 


Modern Physics, by G. E. M. Jauncey. Pp. xiv+561. Third edition. (New 
York: D. van Nostrand Co., Inc., 1948.) 33s. 


The fact that a book reaches a third edition in a comparatively short time shows that it 
has found favour amongst those to whom it was addressed. The task of a reviewer is 
therefore somewhat simplified, for he has only to point out possible improvements and note 
some of the small errors which escape the eye of even the most viligant proof-reader. 

It is noted that photographs of eminent physicists appear throughout the book but, 
unfortunately, the balance of selection appears to have been made in favour of American 
contemporaries : one wonders why Sir J. J. Thomson, Lord Rutherford, Mme Curie and 
Prof. Planck have been omitted ? If it is a question of space, one must ask, what is the 
value of giving a list of Nobel prize-winners or of including Chapter II—‘“‘ Some Useful 
Mathematics’ ? This chapter could easily be omitted from any future edition of the book, 
for there will undoubtedly be a cry for more space, and if a reader is not acquainted with the 
elementary calculus how can he be expected to appreciate, from the help given here, 
discussions on wave mechanics ? 

Now while it is the privilege of any author to use what symbols he likes to denote various 
physical quantities, provided, of course, that each is satisfactorily defined, yet in a book 
which is undoubtedly modern in outlook one is surprised to find that G is still used to denote 
the gravitational constant: why not y, and let G be the strength of the field? It would 
also appear to be a retrograde step to use R as the symbol for the strength of an electric field. 
Usually E is preferred, but if this symbol is used for the charge ona positively ionized particle 
then X [cf. J. J. Thomson’s Elements] could be adopted as the symbol: the use of R is 
particularly unfortunate as in the same book it is used to denote, among other things, the 
universal gas constant and an electrical resistance. "The reviewer would like to see the 
symbols (e/m)_ and (e/m), used to denote the specific charge of an electron and a positively 
charged ion respectively when it is necessary to differentiate between them. 

Again, should not the idea of “‘ lines of force ’’ be omitted from a book for senior students ? 
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If B is the magnetic induction, then | B.n dS, where fis a unit vector, etc., is the flux of B 
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across S, and this is all that the phrase “ lines of force’ would seem to imply. We also 
find that the “‘ gauss ”’ is used as the unit of magnetic field strength and that the units for the 
universal gas constant are given wrongly : they should be erg.mole-tdeg!. In connection 
with units, a lack of consistency in style is noted : one would have thought that an American 
author would have followed the recommendations of Birge and not have written c=3 x 101% 
cm/sec. and J =4-18 joules/cal. [p. 539], but c=3 x 101° cm. sec-! and J =4-18 joule.cal—. 
Also specific heat should have for its unit cal.gm~ deg- c.; the definition of specific heat, 
on p. 165, is wrong, dimensionally. 

The argument on p. 107 is very confusing. ‘The author is endeavouring to calculate the 
capacitance per unit length of a cylindrical condenser so that since dV is always an increase in 
potential this must occur when an element of positive charge is brought nearer to a positive 


a 


charge : in this instance dV = —(2q/r) dr and v=] (—Qair) dr) =2q log (b/a). 


There is sull no discussion of the Zeeman effect : this omission was pointed out by a 
reviewer of the second edition, and several errors and misprints occur, the most amusing of 
which is: “‘ In 1803 Rumford . . . . married the widow of Lavoisier .... Lavoisier had been 
guillotined in 1894....” [p. 24]. One also finds that on p. 171 a cross-reference exists to 
matter which has been removed from the present edition. Finally, why retain the term 
i Baas gas ’’ when an “ ideal gas’”’ is meant, and why print van der Waals with a capital 
BOM eat 

These are only small defects in what will still prove to be an exceedingly useful book to 
serious students of physics, but such errors are so common among writers and teachers that 
I have ventured to draw attention to them. C. J. SMITH. 
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